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Foreword 


It is almost 200 years, since Georg Ohm published his research on the galvanic 
response of a metal to an applied voltage, amid a storm of controversy. He gave us, 
one of the handfuls of physics equations that everybody learns, V — IR. Now, we 
prefer to relate the linear galvanic response to the electric stimulus, via the intrinsic 
electrical conductivity, o, of the metal, j = cE. Some 20 years later, William 
Thomson discovered both the quadratic decrease of conductivity of a normal metal 
in a magnetic field — the normal magnetoresistance — and its variation with the angle 
between the electric current and the magnetization M, in a ferromagnetic metal. 
Near the end of the 19th century, Edmund Hall found that in a perpendicular 
magnetic field, a current excites a transverse voltage proportional to the current in a 
normal metal, and proportional to magnetization in a ferromagnet — the normal and 
anomalous Hall effects. 

Shortly afterwards, the electron with its tiny mass and unchanging quantum of 
negative charge —e was identified as the mobile carrier of electric current, but its 
ability to transport the angular momentum associated with its intrinsic spin of h/2 
discovered in the 1920s, attracted little attention. This was because, unlike the 
charge, the angular momentum of the electron could be flipped in a collision with 
another electron. Spin diffusion lengths, measured in nanometres, are a small 
multiple of the electrons’ mean free path. Anisotropic magnetoresistance and 
anomalous Hall effect became familiar effects but remained unexplained for much of 
the 20th century. Only when it became possible, in the 1970s to prepare high-quality 
magnetic thin films and heterostructures, thinner than the spin diffusion length, did 
‘spin electronics’ become a practical possibility. 

Since then, there has been an avalanche of discovery of new magneto-electric 
phenomena, and spintronics has found important applications in contactless 
sensing, scalable non-volatile memory, and fast electronic switching. A bewildering 
array of new ideas and phenomena has emerged, many associated with spin-orbit 
interaction. Topology in direct and reciprocal space is an important consideration. 
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IV Foreword 


A page is needed to list the acronyms, let alone explain the physical effects. There is 
much for a newcomer to the field to master, both conceptually and practically. 
Vincent Baltz's new book is a welcome guide to the first aspect, for both newcomers 
and practitioners of the art of spin-dependent electron transport. 

Here is a concise, meticulously-illustrated account of the subject with necessary, 
but not excessive mathematical detail, which will allow the reader to grasp the basic 
concepts and learn how to use them in the ten extended exercises, that form an 
integral part of the text. A common notation and a single system of practical SI 
units is adopted throughout, which helps to reduce the confusion of conventions and 
units, found in the literature. Familiarity with the numerical values of the quantities 
involved will allow the reader to develop a critical, physical feel for the subject. 
Baltz's book is based on a series of lectures, given to students in Grenoble, in recent 
years. It will serve as a Baedeker for travellers in the rugged and testing terrain of 
contemporary spintronics. 


Michael COEY 
Dublin, November 2022 


Preface 


'This book is intended for readers, who wish to acquire a solid basis for under- 
standing electron transport in spintronics and the fundamental principles of some 
associated applications. It provides the reader with sufficient knowledge, to be able 
to invest further in this field of research. To this end, some of the prominent key 
notions widely used in the field were selected with care, with the aim of providing a 
simple, concise, and efficient framework. These selected notions are explained, using 
simple examples and analytical calculations. The technical terms and specialist 
jargon are explained, and the way subtleties complicate the phenomena, without 
altering their physical basis, is addressed. 

I am grateful to Mike Coey of Trinity College Dublin, and Daria Gusakova, 
Gilles Gaudin, and Matthieu Jamet of SPINTEC Grenoble, who took the time to 
read this book and gave me to benefit from their opinions, advice, and to correct 
some of the errors. I am very thankful to many students, as well as to colleagues 
from SPINTEC, and within the scientific community, for numerous stimulating 
discussions and for motivating my curiosity and interest in this topic. Finally, I 
would like to thank Héléne Béa and Olivier Fruchart, for allowing me to teach 
spintronics to complete their courses on nanomagnetism, at the University of 
Grenoble Alpes. 

The text has grown out of lectures, given to MSc and PhD students. Therefore, it 
is designed for advanced graduates and postgraduates, as well as researchers and 
engineers, with a background in condensed matter physics and magnetism. Inter- 
ested readers are encouraged to complement their knowledge by consulting, for 
example, the following books: N. W. Ashcroft and N. D. Mermin, Solid State Phy- 
sics, Saunders College Philadelphia (1976); C. Kittel, Introduction to Solid State 
Physics, John Wiley & Sons (2004); J. M. D. Coey, Magnetism and Magnetic 
Materials, Cambridge University Press (2010); E. du Tremolet de Lacheisserie et al. 
(eds) Magnetism I & II, EDP Sciences (2002). 
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VI Preface 


The text is structured in seven parts. Each part feeds on the previous one. 
A careful effort has been made, to indicate how the parts and the physical 
phenomena, they describe, are related to each other. While the text focuses on the 
fundamental aspects, the implementation of the fundamental physical effects in 
typical applications is described in detail. Orders of magnitude of the various 
parameters and key phenomena are provided in each section, based on experimental 
data. The references are concise, in order to provide a reasonable initial framework 
for the interested reader, to explore further. Physical effects, that are more 
peripheral to current concerns are mentioned wherever relevant, with references. An 
Index of key concepts allows the reader to navigate from one key area to another. 
Throughout this textbook, we use SI units. A List of symbols and units, as well as 
the corresponding formulas, is provided so that the reader can refer to them 
throughout the reading. When there are several possible definitions in the literature, 
the definition used in this book is made explicit, and the differences in prefactor 
depending on the articles and books consulted, are also made explicit. 

Chapter 1 gives a concise overview of spintronics and what this discipline 
contributes to society. 

Chapter 2 lays the foundation of spin-dependent electron transport. Based on 
the two-current and sd scattering models, we first show why the electron scattering 
probability is spin-dependent, and how this effect impacts electron transport, using 
the example of the current-in-plane (CIP) giant magnetoresistance (GMR) effect. 
We then show how the addition of spin-orbit interactions mixes the spin states and 
reshuffles the spin-dependent scattering probabilities, using the anisotropic 
magnetoresistance (AMR) effect as an example. 

In chapter 3, we describe in detail, the effect of spin accumulation, encountered 
whenever electrons flow across an interface, due to the distinct partial current den- 
sities in materials of different types. Here, we detail how the electron dynamics are 
described in this case, and what kind of spin-flip relaxation mechanisms prevail. We 
elaborate on the conditions of a key parameter known as the spin-coupled interface 
resistance, necessary to maintain the spin polarization across an interface, and pre- 
sent how spin accumulation is at the heart of the current perpendicular-to-plane 
(CPP) GMR effect. How intrinsic interface effects, such as spin memory loss, as well 
as how non-collinearity and non-uniformity alter spin accumulation, are presented. 

Chapter 4 focuses on the process by which spin angular momentum can be 
transferred from current (spin angular momentum flow) to magnetization, and the 
type of torque (STT), that this transfer generates. The key parameters involved, 
including the spin mixing conductance, are discussed in detail How STT alters 
electron transport, by providing an additional relaxation channel for the spins, is 
tackled, as well as how STT can trigger oscillations and magnetization reversal. 
The spin pumping (SP) reciprocal effect of STT is introduced. Guidance on the 
morphology of STT in non-uniform magnetic textures, as well as hints on the 
magnetoelectronic circuit theory widely used in current spintronics, will be given. 
This chapter ends with a section about spin-orbit torques (SOT), which result from 
the transfer of angular momentum between lattice and the electron orbitals (crystal 
field potential), electron orbital and spin (spin-orbit coupling), and spin and 
magnetization (sd-exchange interactions). 


Preface VII 


In chapter 5, we describe how the intrinsic crystal and spin symmetries, repre- 
sented by different symmetry groups, are particularly important notions, that have a 
significant impact on spintronics. With the examples of the different Hall effects, be 
they unquantized or quantized, we present the Berry formalism, and what it implies 
for intrinsic physical effects. The importance of breaking spatial or temporal 
inversion symmetries is detailed and illustrated. 

Chapter 6 provides a series of ten comprehensive exercises with solutions. 
'They are specifically designed to illustrate the ideas in the previous chapters and 
lead to other spintronic effects based on these ideas. The exercises deal with 
AMR (chapter 2), domain wall AMR (chapter 3), the drift-diffusion equation for 
spin accumulation (chapters 3-5), spin conductivity mismatch in CPP-GMR 
(chapter 3), the intrinsic intra-band scattering contribution to damping of 
magnetization dynamics (chapters 2 and 4), spin pumping and the inverse spin 
Hall effects (ISHE) (chapters 4 and 5), the extrinsic spin pumping contribution 
to damping (chapter 6), spin Hall magnetoresistance (SMR) (chapters 2-4), the 
harmonic analysis of the anomalous Hall voltage and related torques (chapters 4 
and 5), and the intrinsic anomalous Hall and Nernst effects (AHE, ANE) 
(chapter 5). 

Finally, chapter 7 concludes the book with a non-exhaustive presentation of some 
current topics, to which the readers may wish to turn, building on the previous 
chapters, which now allow them to go further, and understand other related or more 
elaborate spintronic phenomena. 


Vincent BALTZ 
Grenoble, November 2022 
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Chapter 1 
Introduction 


Spintronics consists of the investigation and exploitation of the electron spin 
degree of freedom in electronics. The multidisciplinary nature of spintronics 
requires familiarity with general aspects of condensed matter, notably in nano- 
magnetism |J. M. D. Coey, Magnetism and Magnetic Materials, Cambridge 
University Press (2010); E. du Tremolet de Lacheisserie et al. (eds) Magnetism I 
€& IL EDP Sciences (2002)] and spin transport, which are complementary 
disciplines. 

Spin 45 is one of the internal quantum properties of the electron, just like mass 
m, or electric charge —e. It has no classical equivalent, but, by analogy with 
orbital angular momentum, it is often equated with classical angular momentum 
(mer x v), based on the rotational movement of a uniform mass with velocity v on 
a circle of radius r (figure 1). The unit of spin is, therefore, the same as classical 
angular momentum, N.m.s ! or J.s. The spin of an electron is quantized and can 
take two values .s = +h/2 where f is Planck constant h divided by 27. These two 
E are related to the spin projection quantum number, m, = +1/2, with 
=+m,h, for a spin quantum number s = 1/2 and a corresponding amplitude, 


sz hy/ s(s -- 1). 

Note: throughout the text, bold is used for vectors, r = 7. All symbols and units 
are listed on page 150. 

Because an electron has a spin and a charge, it also has a magnetic moment m, 
just like a rotating electrically charged body in classical electrodynamics, which 
gives rise to a magnetic dipole moment, (x — $r x v) (figure 1). 

'The actual relation between the spin and the magnetic moment of an electron is 


derived from quantum mechanics 


€ 


m= — g8 (1) 


2m, 


Compared to classical mechanics, the result of quantum mechanics involves the 


Landé g-factor due to relativistic effects. For an isolated electron, g ^» 2. An electron 


thus takes values m ~ F up, with ug = i the Bohr magneton, in A.m? 


DOI: 10.1051/978-2-7598-2917-0.c001 
© Science Press, EDP Sciences, 2023 


2 The Basics of Electron Transport in Spintronics 
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m 


Fic. 1 — Illustration of the electron spin 8 and magnetic moment m, associated with the 


rotation of a uniform distribution of charge over mass +, represented by the arrow in the 


Me? 


horizontal plane. 


The expression above is often used as, 


m= ys = — 


|s (2) 


where y= —3--g«0 is the gyromagnetic ratio. For an isolated electron, 


D 28 GHz.T-!. 

In matter, the total spin angular momentum takes into account the spin and 
orbital angular momenta, as well as the coupling between the two, called spin-orbit 
coupling, and possibly, some quenching of the orbital angular momentum due to 
crystal symmetry. The resulting total spin angular momentum and the corre- 
sponding magnetic moment take on discrete values, the amplitudes and projec- 
tions of which are also defined by a principal j and a projection or secondary m; 
quantum number. The relation between the total spin angular momentum and the 
total magnetic moment takes the same form as equation (1), except that the value 
of g now takes into account spin and orbital contributions, and possibly, spin-orbit 
coupling and some orbital quenching, for example, gp,~ 2.09, gyi 2.19, and 
Ico ^» 2.15. 

It is important to note that, in the literature, the term spin is sometimes used 
interchangeably for spin magnetic moment m and spin angular momentum 8, which 
can lead to confusion as m and 4 are of opposite sign due to the negative charge — e of 
the electron. In this book, the term moment will refer exclusively to m and spin to 5. 
This distinction is crucial when dealing with the accumulation or transfer of angular 
momentum, particularly with regard to the sign of the coupling between several 
angular momenta, at the interface between two materials. 

Spin transport, especially when conduction electrons carry out spin, is the 
subject of this book. In fact, spin information can be conveyed in two direct 
ways (figure 2), via: 
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W electrons, i.e., via the spin of the conduction electrons, thus complementing the 
charge transport — although the net charge transport may be zero in some cases 
(chapter 2), 


WB magnons, i.e., via the exchange interactions between the spin of the localized 
electrons in the crystal lattice, giving rise to so-called spin waves, quantified by 
elementary excitations called magnons. Magnon-mediated transport may be 
considered to be a separate subject, beyond the scope of this book. The interested 
reader can find more information in S. Rezende, Fundamentals of Magnonics, 
Springer (2020). 


We note that starting from conduction electrons and magnons, there are 
several indirect ways to convert and transport spin. More generally, spin is a form 
of angular momentum tied to electrons but any excitation that possesses angular 
momentum (chapter 2), like phonons and photons, has the ability to be a carrier. 
Any angular momentum transfer process will be allowed via some interconversion, 
e.g., electron-electron (chapters 2-4), electron-magnon, electron-phonon, and 
electron-photon. 


3 


| T arm 


, / y 
' of 


Fic. 2 — Illustration of how spin information can be conveyed, by (Left) itinerant conduction 
electrons and (Right) collective excitations of the localized moments in the crystal lattice. J, 
represents the spin current density, and the arrow pointing up indicates the spin angular 
momentum 8, which gives the polarization of Js (see chapter 2). Note: in the illustrations, the 
net charge current is zero, which is a special case detailed in chapter 2 for itinerant conduction 
electrons. Reprinted with permission from Springer Nature: Y. Kajiwara et al., Nature 464, 
262 (2010). Copyright 2010. 


Js 


The field of spintronics builds on several success stories that have contributed to 
technology and to society. Two of them trace the scientific path from two funda- 
mental discoveries to industrial applications, involving magnetic sensors that are 
sold in the hundreds of billions: 
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E from the giant magnetoresistance (GMR) effect to sensors for automotive and IT 
applications 


The understanding of the basics of electron transport in relation to the GMR 
effect and the use of this fundamental effect in some applications will be dis- 
cussed in detail in chapters 2 and 3. The chronology of the success of GMR is: 


1988 Discovery of the GMR effect [M. N. Baibich et al., Phys. Rev. Lett. 61, 2472 
(1988); G. Binasch et al., Phys. Rev. B 39, 4828 (1989)] 


1991 Development and practical implementation in spin-valves |B. Dieny et al., 
J. Magn. Magn. Mater. 93, 101 (1991); S. S. P. Parkin, Phys. Rev. Lett. 64, 
2304 (1990)] 


1997-present Commercial product — sensors for automotive and information 
technology (IT) as billions of hard-disk drives — since supplanted by tunnel 
magnetoresistance (TMR) junctions in some applications [demonstration 
of TMR: J. S. Moodera et al., Phys. Rev. Lett. 74, 3273 (1995); the 
proposition of orbital-filtering for larger TMR values: W. H. Butler et al., 
Phys. Rev. B. 63, 054416 (2001)] 


2007 Nobel prize in physics awarded to A. Fert and P. Grünberg [https:/ / www. 
nobelprize.org/] 


Bl from the spin transfer torque (STT) effect to memories for IT applications 


The basics of electron transport in relation to the STT effect are discussed in 
detail in chapter 4. The chronology of the success of the STT effect is: 


1995 Discovery of the STT effect [J. C. Slonczewski, J. Magn. Magn. Mater. 159, 
L1 (1996), L. Berger, Phys. Rev. B 54, 9353 (1996)] 


2005 Practical implementation in TMR junctions 


2013-present Commercial product - MRAM (magnetic random-access 
memory) [B. Dieny et al. (eds), Introduction to Magnetic Random-Access 
Memory, Wiley-IEEE Press (2016)] 


'There are still more successes to come, not least because of the emergence of new 
fundamental phenomena, such as pure spin current physics (chapter 2), spin-orbit 
torques (chapter 4) and Berry curvature physics (chapter 5), and the commercial- 
isation of new MRAM devices of several types, including SOT-MRAMs by major 
microelectronics companies. The miniaturisation of electronics is hampered, in 
particular, by problems of power consumption. Microprocessor chips require about 
100 W.cm ? to operate, which is an order of magnitude higher than a kitchen stove. 
In addition, energy consumption in a chip is aggravated by leakage [C. Singh and 
R. Tangirala, IBS report (2013)] and by the volatile nature of the data storage. 
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Spintronics is ideally positioned to circumvent energy consumption issues, as it can 
combine non-volatility with dissipation-free transport of spin information when it 
uses pure spin currents. Energy harvesting spintronics is predicted to be a viable 
potential technology for many application areas [A. Hirohata et al., J. Magn. Magn. 
Mater. 509, 166711 (2020)], for example: 


E Information technology — IT, e.g., memory, processors, data security 
E Biomedicine, e.g., sensors 

E Telecommunication, e.g., transceiver 

E Artificial intelligence — AI, e.g., neuromorphic computing 


'This wealth of applications of interest to society promises a bright future for the 
discipline, whether approached from the point of view of fundamental physics (see 
chapter 7 for a description of some of the current fundamental topics) or from an 
applications perspective. To further emphasise how magnetism and spintronics are 
and will be present in our lives, I would quote M. Coey: *More transistors and 
magnets are produced in fabs than grains of rice are grown in paddy fields". 

It seems important to emphasize that fundamental research and the industrial 
developments that enable its use in applications are complementary and equally 
essential. While some advances are the result of responses to current societal issues, 
other advances are the result of Blue Sky research, the objective of which is to 
advance the state of knowledge, independently of the pressing issues of the day, and 
which will nevertheless be essential for responding to tomorrow's challenges. 

Motivated by the potential of spintronics, and because thinking outside the box 
requires a knowledge of what is already in the box, we are now ready to dive 
enthusiastically into the following chapters, which focus on the fundamental aspects 
of the discipline with the aim of providing the basic knowledge to understand most 
of the underlying physics. 


Chapter 2 


First Notions About Electron and Spin 
Transport — CIP-GMR, IEC, AMR 


When a system is brought out of equilibrium by an external stimulus, or force, 
dissipative transport processes cause it to relax back toward equilibrium. If the 
system is sustained out-of-equilibrium by a stationary force, it reaches a steady 
out-of-equilibrium state with steady transport processes. If not, the system relaxes 
back into equilibrium. Here, the word system combines electrons and their angular 
momenta. External forces can, for example, be caused by (electronic) potential 
(§2.2), chemical (§3.2), and thermal gradients. Electron and spin angular momen- 
tum relaxations are intertwined. These relaxations are due to the alteration of the 
crystal periodicity, resulting in scattering by crystal defects, lattice vibrations 
(electron-phonon interaction), electrons (electron-electron interaction, §2.4, §3.2 
and $4.1), and magnetic vibrations (electron-magnon). We recall that all the 
relaxation processes listed above can be spin-dependent because electrons, magnons 
and phonons can host angular momentum. In this book, we will mostly deal with 
(spin-dependent) electron-electron interactions. 

How electron and spin dynamics are described in spintronics and which type of 
relaxation mechanism prevails depends on the material's properties and device 
geometry, as will be presented throughout the text. We will first introduce the types 
of models that are used, their accuracy and their complexity. 


2.1 Types of Model 


'The models used to describe electron and spin transport are classified into four 
types, based on how close they are to a quantum mechanical description (figure 3): 


E Phenomenological models 


This type of model describes empirically the relationship between causes and 
effects. We will, for example, use such models to illustrate anisotropic 
magnetoresistance (AMR, §2.6) and giant magnetoresistance (GMR, §2.5), asa 
forenurmer to the more sophisticated models presented below. 
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W(r,t) AMR ‘CFI’ 82.6, STT ‘Slonczewski’ 54.2, Berry curvature 55.1 ... 
f(r,k,t) GMR drift-diffusion ‘Valet-Fert’ 3.3 ... 

n(r,t) Drift-diffusion ‘Drude’ 82.2 ... 

All: AMR, GMR, TMR ... 


Fic. 3 - Illustration of several types of models used in spintronics, based on how close they are 
to a quantum mechanical description. 


E Classical models 
'This type of model tracks deterministically the position of electrons at any time, 
like in the free electron model, where a position (r)- and time (t)- dependent 
function can be used, for example, to describe the density of electrons n(r, t). An 
example of a classical model for electron dynamics is the Drude model. It will be 
briefly evoked in $2.2 and $3.2 because the semiclassical descriptions of 
spin-independent and spin-dependent transport are connected to it. 


E Semiclassical models 

In semiclassical models, part of the system is described quantum mechanically. 
For example, the semiclassical Boltzmann model uses a space-, time-, and 
momentum(k)-dependent distribution function g(r,k, t), which corresponds to 
the Fermi-Dirac distribution function at local equilibrium /(r,k, t), therefore 
involving the dispersion of Bloch electrons e(k), where e is the electron energy. 
The various physical quantities, like the distribution of conduction electrons and 
the current, are obtained by performing integrals over the distribution function. 
'The semiclassical Boltzmann transport frame is used to model electron transport 
like, for example, in the Valet and Fert model for GMR. We will deal with such 
models in $2.2 and $3.3. Another example of the use of Boltzmann transport is 
found in $5.2 when deriving the transverse conductivity inherently related to 
crystal symmetry. 


E Quantum models 
A quantum model is self-consistent. The state of a system is described by a space- 
and time-dependent vector, via a wavefunction w(r, t) when electron transport is 
at play and via a spinor when spin angular momentum is also considered. Its time 
evolution obeys the Schrödinger equation, which yields the eigenstates y(r, t) 
and the energy eigenvalues e(k). Physical quantities, like charge and spin current, 
are represented by observables, and their values are associated with the eigen- 
values corresponding to the value of the observable in a given quantum state or 
eigenstate. A toy quantum model for spin transfer torque (STT) will be pre- 
sented in $4.2, in which, expressions for the above-mentioned observables are 
detailed. Quantum mechanics is also used when dealing with the notion of Berry 
curvature in $5.1, and AMR in $2.6. Other examples of quantum models used in 
spintronics are the Kubo formalism, which deals with the linear response of an 
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observable due to a time-dependent perturbation, and the Keldysh formalism for 
the description of the evolution of a system from a non-equilibrium state. 


Next, we will introduce why electron transport can be spin-dependent, using a 
semiclassical description. 


2.2 Two-Current Model, Ohm's Law, Spin-Dependent 
Mean Free Path 


Electrons are split into two categories according to the orientation of their 
moment /spin. Those two categories are called majority moment({)—spin({/) electrons 
and minority moment(|)-spin({}) electrons. It is common to consider the term 
majority for electrons carrying a magnetic moment 4n! ‘pointing up’, along the 
magnetization in a ferromagnet (as detailed in $2.3), and thus a spin angular 
momentum 8! ‘pointing down’ from s = —4n/|y|. Consequently, the minority 
moment(|)-spin(f) electrons carry a magnetic moment m!, and a spin angular 
momentum 8". 

We recall that, in the literature, the terms moment m and spin 8 are sometimes 
used interchangeably, and the 1 and | symbols sometimes refer to m or 8, which can 
lead to confusion as m and 4 are of opposite sign. In this book, the term moment refers 
exclusively to m and spin to 8; the 1 and | symbols refer exclusively to the orientation 
of m, and the fr and J) symbols refer exclusively to the orientation of 8. Note that, in 
the literature, the + and — symbols are sometimes used for the orientation of 5. 

For now, we also consider that majority moment(T)-spin()) electrons and 
minority moment(|)-spin(f)) electrons flow in separate channels, according to the 
two-current model introduced in 1936 by N. F. Mott (figure 4, left). Their (charge) 
current densities are j! and j!, respectively. 

We note that how the two spin channels mix is described in some of the next 
sections. Mixing can, for example, occur because of spin-orbit interactions (82.6), 
spin-flip scattering (83.3), and transfer of angular momentum ($4.1). 

Because electrons carry a charge and a spin (angular momentum), these two 
entities naturally flow in matter during electron transport. Therefore, both a charge 
(Je) and a spin (Js) current density must be considered 


e =j! +j! (3) 


J; - - (f -7') (4) 


In the literature, the spin current density Js = -( — j) is expressed either in 
A.m~’, units of charge current density, as here and throughout this textbook, or in 


J.m?, units of density of spin angular momentum am 
Js = -&(f = j^). resulting in a factor i difference is some of the equations 


depending on the articles and books consulted. In this textbook, we use Q = iJ. 


in which case 
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Fic. 4 — (Left) Representation of electron scattering in the two-current model. A, is the 
electron mean free path. Note: to make the representation easy to read, 4, is shown between 
two particular collisions. In practice, 4, is averaged over all collisions. (Right) Schematics of 
the types of charge and spin currents encountered in spintronics, reprinted with permission 
from Wiley: Y. F. Feng et al., WIREs Comput. Mol. Sci. 7, e1313 (2017). Copyright 2017. 


In the literature, the spin current is also defined either as J, —j! —j', or as 
J= -( j- i‘), resulting in a factor —1 difference in some of the equations 
depending on the articles and books consulted. Throughout this textbook, the term 
spin will refer to 8 and the spin 4 current is thus defined as J, = j? — j^ = — ( ji- i) ; 
because the f and | symbols refer to the orientation of the moment m, whereas the ff 
and |J symbols refer to the orientation of the spin 8, which points in the opposite 
direction to m (figure 1). 

Charge and spin currents can take very different values, particularly because the 
net motion of electrons in a given direction is the result of an averaging of instan- 
taneous motions, which may have components in all directions in space. We dis- 
tinguish the following types of currents (figure 4, right): 


E An unpolarized current, when the current is made of an equal number of 
majority moment(T)-spin(i) electrons and minority moment(|)-spin({}) 
electrons, all flowing in the same direction. In this case, Je Z 0 and J, = 0. 


W A spin polarized current, when the current is made of an unequal number of 
majority moment(T)-spin(4) electrons and minority moment(])-spin(1) 
electrons, all flowing in the same direction. In this case, Je Z 0 and J, Z 0. 


WB A fully spin polarized current, when the current is made of majority 
moment(T)-spin(|)) electrons only or minority moment(]|)-spin(f)) electrons 
only, all flowing in the same direction. In this case, Je #0 and J, is 100% 
polarized. 
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E A pure spin current, when the current is made of an equal number of majority 
moment(T)-spin(|.) electrons and minority moment(|)-spin(‘}) electrons, but 
flowing in opposite directions. In this case, Je = 0 and J, is 100% polarized. 


The main difference between a charge and a spin current relates to how they obey 
conservation laws and subsequently to how charge and spin accumulate, as detailed 
in §3.2. We note that most of the text is based on a one-dimensional description. 
This simplification was chosen to allow the reader to focus on the actual physical 
phenomenon at play. Three-dimensional descriptions and the associated complica- 
tions are discussed in §3.5. 

Next, we recall how transport is related to an external stimulus, such as an 
electric field E. It is possible to demonstrate that the semiclassical treatment of 
electron transport leads to a macroscopic description or Ohm’s law [equation (5)], in 
which the current densities are related to the electric field, and thus to the electric 
potential, via the macroscopic conductivity o, 


j-cE (5) 
with, from the Drude-Sommerfeld model, 
N Pta 
ga ee. N(er)e?D (6) 
Me 


where N(er) is the density of states (DOS) at the Fermi level, v, is the electronic 
relaxation time (related to the electron mean free path 4, via the Fermi velocity vp: 
Ae = vpt,), and D is the electronic diffusion coefficient. 

Using E = —V$, with $ the electric potential, and defining the electrostatic 
potential as 4, — —ed, the relation between current density and electrostatic 
potential can be obtained from equation (5) 
al, 


€ 


= (7) 
here, with u= ue. This type of macroscopic description remains valid when 
spin-dependent coefficients are involved [equation (8) and $2.3] and when electro- 
static potentials are extended to electrochemical potentials [equation (22)], when 
U = Ue us. The validity of this assertion in the most general case will be discussed in 
more detail in $3.2. The correct term is electrostatic potential energy because it is a 
potential energy, the unit of which is joule, but it is customary to use the abbreviated 
term electrostatic potential, just as the term electrochemical potential describes a 
potential energy. 

In the next part, we will explain why te (or 2e) is spin-dependent in ferromagnets, 
therefore making j! Æj! [from equations (6) and (7)], resulting in a spin-dependent 
generalized Ohm's law 


10) 
ot) YH 


€ 


(8) 


Notation: equation (8) is two concatenated equations: j! = c! ux andj! — c! n 
'This type of notation is used throughout. 
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2.3 Band Structures, Spin-Dependent Fermi Surface 
& Density of State 


In metals, electronic conduction occurs at or near the Fermi surface, i.e., for € = €f 
(figures 5 and 6). Due to band splitting, ferromagnetic metals have different Fermi 
surfaces for majority moment(T)-spin(lj)) electrons and minority moment(])-spin 
(fr) electrons (figure 6), making electron transport spin-dependent. We recall that 
band splitting causes a gain of Coulomb energy, because some electrons are 
unpaired, at the expense of a loss of kinetic energy, because electrons move up in 
energy. The overall energy balance is favorable for densely packed electrons when a 
small loss of kinetic energy leads to a large gain of Coulomb energy. In turn, the 
delocalized s-electrons have no tendency for band magnetism, in contrast to the 
of-electrons of the early actinides and the 3d-electrons of the late transition metals, 
which possess sufficiently narrow bands. Narrower bands favor localized (sometimes 
almost atomic) magnetism, like the late actinides and rare earth compounds. In this 
book, we focus on band magnetism [see also Stoner criterion in J. M. D. Coey, 
Magnetism and Magnetic Materials, Cambridge University Press (2010)]. 

The relation between exchange splitting and spin-dependent transport is pre- 
sented next, based on the partial densities of states DOS, N' (figure 7). A DOS 
represents the number of continuum states in an infinitesimally small energy interval 
e+ de. In 3d-transition metals, spins contributing to transport are split into two types: 


WB itinerant (delocalized) spins carried by light 4s-electrons, corresponding to the 
vertically elongated parts of the DOS in figure 7, which are not spin-split in 
energy, and 


W localized spins carried by heavy 3d-electrons, corresponding to the horizontally 
elongated parts of the DOS in figure 7, which are spin-split in energy for the Fe, 
Ni, and Co ferromagnets in the figure and unsplit for non-magnetic Cu, where 
the 3d-bands are full. 


Copper Aluminium 


3d!04s! 3s?4p! 


Fic. 5 — The Fermi surfaces of two non-magnetic metals: copper and aluminium. From 
http://www.phys.ufl.edu/fermisurface/. 
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Fic. 6 — The Fermi surfaces of (Left) bcc-iron (body-centered cubic — Fe) and (Right) 
fcc-nickel (face-centered cubic — Ni) for majority moment(1)-spin(4)) electrons and minority 
moment(|)-spin(f) electrons. From http://www.phys.ufl.edu/fermisurface/. The colors 
correspond to the contributions of several bands. 


As a result, for spin-split ferromagnets, the DOS at the Fermi level is 
spin-dependent, for example, N! (ep) « N' (er) for hcp-Co. 

The terms ‘light’ and ‘heavy’ come from the effective mass m7 theory used to 
simplify band structures by modelling the behaviour of a free electron with this 
mass: e(k) — FE E and mi = = WP / Zs. 

Note "E. ans spin polarization obtained from the DOS at the Fermi level 


ic depends on how the bands are split, explaining why hcp-Co and bcc-Fe 


have opposite DOS spin polarization (figure 7). Note also that, transport ( j) relates 
to processes occurring near the Fermi level whereas magnetism relates to properties 
up to the Fermi level, i.e., magnetization (M) relates " iu idi n the bottom 
of the band to the Fermi level M = up [^ f(e)N — fE f(e) N'(s)de, with 


fle) = 1 the Fermi-Dirac distribution A it. bci potential 


defined in $3.2. The majority moment(1)-spin(4J) electrons can be in the minority 
near the Fermi level compared to the minority moment(|)-spin(f)) electrons which 
are a majority, see e.g., hcp-Co in figure 7. 

How the vertically elongated unsplit partial DOS of the itinerant s-electrons and 
the horizontally elongated spin-split partial DOS of the localized d-electrons are 
represented in a simplified manner are shown in figure 8. 
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Fic. 7 — Energy(s — £p)-dependence of the partial density of state DOS, NU of 3d transition 
metals for majority moment(T)-spin(|)) electrons (left quadrant in all 4 graphs) and minority 
moment(|)-spin(1) electrons (right quadrant), obtained by ab-initio calculations. The Fermi 
level is represented by the horizontal dashed line at e — ep = 0. Note that, here, N is given in 
eV !.atom '.spin ! compared to J !.m ? in the text. Adapted with permission from Springer 
Nature: J. Mathon & A. Umerski, Theory of Tunneling Magnetoresistance in J. L. Moran-Lopez 
(eds) Physics of Low Dimensional Systems, Springer, Boston. Copyright 2001. 


0 
€ - £c (eV) 


-5 


e o 


e 

= 
mI 
DOS, N'Y (ev-1.atom.spin') 


Fic. 8 — Illustration of how the energy(e — £p)-dependence of the partial density of state 
DOS, N') is represented in a simplified manner. The colors illustrate that the s- and d-bands 
are filled up to the Fermi level indicated by the dashed line. The image on the left is adapted 
with permission from Springer Nature: J. Mathon & A. Umerski, Theory of Tunneling 
Magnetoresistance in J. L. Moran-Lopez (eds) Physics of Low Dimensional Systems, Springer, 
Boston. Copyright 2001. 
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2.4 Spin-Dependent Scattering, the sd Model 


We now introduce the sd model for which: (i) the current is carried by the 
selectrons, and (ii) elastic scattering between the itinerant selectrons and the 
localized d-electrons dominates. 

We consider electron-electron scattering only, and the two-current model stip- 
ulates that majority moment(1)-spin(4J) electrons and minority moment(|)—spin({}) 
electrons flow in separate channels (§2.2). It follows that a s!-electron can be 
scattered by a s!-electron or a d'-electron, and that a s!-electron can be scattered by 
a s!-electron or a d!-electron. 

According to Fermi's golden rule, the scattering probability P^ of an electron 
in an initial state |i) to a final state |f) is determined by the DOS N” (ep) allowed in 
the final state, after scattering P^ œ (f|#|i)?N/ (er), where H is the scattering 
Hamiltonian. Based on the largest DOS at the Fermi level N!(er) for the minority 
moment(|)-spin(f) electrons d!, it is possible to conclude that the most efficient 
channel for scattering corresponds to an itinerant s!-electron scattered by a localized 
d'-electron (figure 9), therefore making the electronic scattering rate 
spin-dependent: t! > t}. As an example, for Co, we have t! ~ 10r}, with A ^ 10 nm 
and 4. — 1 nm (Note: 4, = vpt,, with vp 2 x 109 m.s! ). 


— d! most efficient 
scattering channel 


Fic. 9 — Representation of electron scattering in the two-current model. Based on the DOS at 
the Fermi level: (i) scattering between the itinerant s-electrons and the localized d-electrons is 
more efficient than scattering between selectrons, no matter the conduction channel 
considered (f or |), and (ii) this type of sd scattering is most efficient for minority 


moment(|)-spin(fr) electrons, because there are more d!-states to receive them. 


Because 1, (or 4,) is spin-dependent, it naturally follows from equations (6) and 
(7), that j! Z j', resulting in the spin-dependent generalized Ohm's law anticipated 


in equation (8) and repeated here: j!) = ol!) 


vu ; 
—— (two equations). 
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Fic. 10 — (Top) Illustration of how asymmetric electron scattering, due to spin-dependent 
DOS, is represented in an equivalent circuit of two parallel resistors. (Bottom) Asymmetry 
parameter « = o! /o! experimentally determined for impurities of 3d transition metals in Ni. 
Adapted with permission from EDP Sciences: A. Fert, Reflets Phys. 15, 5 (2009). 


Spin-dependent electron scattering due to band splitting is also referred to as 
asymmetric electron scattering, and the resulting spin-dependent generalized Ohm's 
law [equation (8)] can be represented in an equivalent circuit of two parallel resistors or 
resistivities, wherein the incoming charge current Je = j! +j! splits in two (figure 10). 
A high resistivity value p! = 1/6! accounts for a large number of sd scattering events 
in the minority moment(|)-spin(f) electron channel hosting the j^ current, and a 
smaller resistivity value p! — 1/6! is used for less sdscattering in the majority moment 


(t)-spin({/) electrons channel hosting the j! current. The total resistance of this 
equivalent circuit is expressed as 


tol 
pp 
i gm (9) 
pi+p 


p 


It can be rearranged to account for the spin polarization f of the current at the 
Fermi level 


p=- B’) (10) 
with p! = 2p*(1— (+ )£), and 


pl-pl a-1 i 
~— pttpl atl un 
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where « is the asymmetric scattering parameter 


NA 
pi 


We will see other examples of the use of f later, in the context of CIP-GMR (82.5), 
and of the injection of a current across an interface ($3.2) and the subsequent transfer 
of angular momentum ($4.3). The spin polarization of the current relates to the 
parameter « also known as the spin asymmetry parameter as it relates to asymmetric 
electron scattering due to the spin-dependent DOS at the Fermi level. It is a 
material/element-dependent parameter, which can be determined experimentally 
(figure 10). 

Note that, the generic principle of the tunnel magnetoresistance effect (TMR) is 
based on spin-dependent DOS at the Fermi level, as it drives the tunneling probability 
across the tunnel junction of a ferromagnet /tunnel barrier/ferromagnet trilayer. A toy 
quantum model accounting for wavefunction propagation through a potential barrier 
depicts the phenomenon well. More subtle orbital-dependent tunneling is at play in the 
latest TMR devices. The principles of TMR are extensively described in the literature 
[S. Yuasa et al., J. Phys. D: Appl. Phys. 40, R337 (2007); W. H. Butler, Sci. Technol. 
Adv. Mater. 9, 014106 (2008)]. TMR devices are now extensively used in memory and 
sensor applications (82.5), and still require developments at the engineering level but 
not so much at the basic level, except for the newest materials (chapter 7). 


a (12) 


2.5 From Impurity Scattering to Heterostructures 
— CIP-GMR, IEC 


In the previous section, we discussed how asymmetric electron scattering influences 
transport in a mono-element ferromagnet. In this section, we will consider how 
asymmetric electron scattering influences transport, first, in ternary ferromagnetic 
alloys and then, in heterostructures made of a ferromagnet /non-magnet /ferromagnet 
trilayer [A. Fert, Reflets Phys. 15, 5 (2009)], and discuss a related effect known as the 
interlayer exchange coupling. 


E Element doping control over a stack’s resistance — the example of ternary alloys 


We consider a ferromagnet made of Ni, doped with Co and Rh elements having spin 
asymmetry parameters larger and smaller than 1, respectively: aco = Ph, / Peo -1 
and ag, = plui pls « 1 (figure 11, Type 1, Ni(Co-Rh) alloy). Let us compare it to the 
case of the same Ni ferromagnet but doped with Au and Co elements, both having 
spin asymmetry parameters larger than 1: aco = Ph, / ph, >1 and wy= 
pial p. > 1 (figure 11, Type 2, Ni(Au-Co) alloy). 


In the case of the Ni(Co,_,-Rh,) alloy, the electrons’ flow is altered in the two spin 
channels. In the majority moment(T)-spin(J) electrons channel, electrons are 
(additionally)-scattered by the Rh atoms (on top of the initial scattering by Ni 
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Type 1: the Ni(Co-Rh) alloy 
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Fic. 11 — (Left) Representation of spin-dependent scattering in ternary alloys, made of a Ni 
host and two impurities. Type 1: the two impurities (Co and Rh) have asymmetric scattering 
parameters larger and smaller than 1, respectively. Consequently, they alter electron transport 
in different spin channels, making both channels more resistive. Type 2: the two impurities 
(Co and Au) both have asymmetric scattering parameters larger than 1, consequently, they 
alter electron transport in the same spin channel, making electron flow efficient in the unal- 
tered channel — here, the channel for majority moment([)-spin(|j) electrons. (Right) 
Experimentally determined data. Adapted with permission from EDP Sciences: A. Fert, 
Reflets Phys. 15, 5 (2009). 


atoms), as (my = Pis / Pig « 1. In the minority moment(|)—spin({}) electrons channel, 
electrons are (additionally)-scattered by the Co atoms, as aco = Pb / Ph, >1 
(figure 11, Type 1). As a result, the overall resistivity of the ternary alloy is increased 
compared to expectations, if spin-dependent scattering were discarded. 


In contrast, in the case of the Ni(Au;-,-Co,) alloy, the electrons’ flow is altered in 
one spin channel only. In the majority moment(T)-spin(4)) electrons channel, elec- 
trons are not (additionally)-scattered, and in the minority moment(|)-spin(f) 
electrons channel, electrons are (additionally)-scattered by both the Co and the Au 
atoms, because aco = Poo / Do, > 1 and aay = Phyl Pha > 1 (figure 11, Type 2). Asa 
result, the overall resistivity of the ternary alloy is almost unchanged. 
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Note that, the reshaping of the DOS of the Ni ferromagnet due to the presence of 
dopants explains the use of asymmetric parameters op, and «x4, associated with 
non-magnetic dopants like Rh and Au. In fact, the presence of a non-magnetic 
impurity in a magnetic metal host will create a virtual bond state, with a DOS of 
finite width in energy. This state will hybridize with its magnetic environment. 
Depending on the position of the virtual bond state with respect to the Fermi level of 
the magnetic host, electrons are removed or added at the Fermi level, enhancing or 
diminishing the magnetic asymmetry of the ferromagnet, and conversely conferring 
the impurity a magnetic character, with momentum in the opposite or same 
direction to the host, respectively [J. Friedel, Nuovo Cim. Suppl. 2, 287-3 (1958); 
R. Zeller, J. Phys. F: Met. Phys. 17, 2123 (1987)]. 


Bi Relative direction of magnetization control over a stack’s resistance - CIP-GMR 


In the section above, we have seen that the spin-dependent resistance of a stack can 
be controlled by the nature of dopants. The initial idea underlying the giant mag- 
netoresistance effect (GMR) is to ‘replace’ the dopants with ferromagnetic layers 
(Fs) and the ternary alloy by a F/non-magnet(N)/F trilayer. The magnetizations of 
the two Fs can be manipulated independently in such a way that they can be set 
pointing in opposite directions, resulting in the antiparallel (AP) state, or in the 
same direction, resulting in the parallel (P) state (figure 12). The current is in-plane 
(CIP). Note that, the magnetization of a ferromagnet can be manipulated by an 
external magnetic field or an electric current through spin torques (chapter 4). In 
addition, the two layers can be manipulated independently, if for example, their 
coercive (reversal) field differ or if one of the two layers is pinned by exchange bias in 
spin-valve devices (82.5) [J. Nogues et al., J. Magn. Magn. Mater. 192, 203 (1999)]. 


Inthe F! /N/ F* AP state, the electron flow is altered in the two spin channels due 
to bulk scattering upon crossing the F! / N and N/F'! interfaces to enter the F! and 
F! layers (figure 12, top). In the majority moment({)-spin(J}) electrons channel, 
electrons are (additionally)-scattered near the N / F* interface since x pm pig d, 
and in the minority moment(|)—spin({}) electrons channel, electrons are (addition- 
ally)-scattered near the F'/N interface as ap; = p! /p! > 1. As a result, the overall 
resistivity of the AP state is ‘large’. This behavior is equivalent to the type-1 ternary 
alloy discussed in the previous section (figure 11). From the corresponding equivalent 
circuit (figure 12, top right), the resistivity of the AP state is 


Ap. pi pl 
FUTT 


(13) 


We considered that the two Fs are the same here. Different Fs lead to different 
asymmetric parameters and different spin-dependent resistivities, making the for- 
mula more complex, but leaving the physical explanation unchanged. 

In the F! / N/ F! P state, the electrons’ flow is altered in one spin channel only due 
to scattering near the F! / N and N/ F! interfaces (figure 12, bottom). In the majority 
moment(T)-spin(l)) electrons channel, electrons are not (additionally)-scattered, 
whereas in the minority moment(|)-spin(f) electrons channel, electrons are 
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Fic. 12 — (Left) Schematics of the current-in-plane (CIP) giant magnetoresistance 
(GMR) effect, which is a prototypical example of how asymmetric electron scattering influ- 
ences electron transport — due to spin-dependent DOS and subsequent spin-dependent scat- 
tering near interfaces in a ferromagnet /non-magnet /ferromagnet trilayer. Specular reflection 
is considered at the farthest interfaces. (Right) Corresponding equivalent circuit. Adapted 
with permission from Springer Nature: P. Grünberg et al., Metallic Multilayers, in Y. Xu et al. 
(eds) Handbook of Spintronics, Springer, Dordrecht (2015). Copyright 2015. 


(additionally)-scattered near both interfaces as «pı = o!/p! » 1. As a result, the 
overall resistivity of the P state is ‘smaller’. This behavior is equivalent to the type-2 
ternary alloy discussed in the previous section (figure 11). From the corresponding 
equivalent circuit (figure 12, bottom right), the resistivity of the P state is 
PO alll (14) 
p' p! 


'The magnetoresistance, here known as the current-in-plane giant magnetore- 
sistance, CIP-GMR reads 


AP _ oP tT _ nly 2 _1\2 
CIP-GMR =2 — ^. - (e : ) - E ;) (15) 
p p! +p! a+1 


CIP-GMR depends on spin-dependent asymmetric scattering x due to 
spin-dependent DOS. The important characteristic length for CIP-GMR is the 
spin-dependent electron mean free path, AO (82.3). Neither spin-flip scattering 
(chapter 3) nor spin mixing (82.6) is needed to create CIP-GMR. Those effects can 
alter CIP-GMR, but they are not necessary to trigger the effect. Conversely, we will 
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see in chapter 3 that current perpendicular-to-plane (CPP)-GMR is driven by 
spin-flip scattering and that the important characteristic length for CPP-GMR is 
thus the spin diffusion length l. The fact that I, > MY (roughly a factor of 10 
larger) marks an essential difference between CPP- and CIP-GMR. 

The terminology ‘giant’ refers to the fact that these types of magnetoresistance 
effects, demonstrated in the late 1980ies — early 1990ies, are much larger than the 
anisotropic magnetoresistance effect discovered much earlier, in the late 1850ies — 
modelled in the early 1970ies and due to spin-orbit spin mixing (82.6). A. Fert and P. 
Grünberg were awarded the Nobel prize in physics in 2007 “for the discovery of the 
giant magnetoresistance effect". 


E Note about the interlayer exchange coupling (IEC) 


We next present interlayer exchange coupling (IEC), which is a companion effect to 
CIP-GMR. It is also closely related to Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction between magnetic impurities. 

For CIP-GMR, the thickness dw of the non-magnetic (N) spacer layer is fixed 
and the relative alignment of the ferromagnetic layers is varied, for example, by use 
of an external magnetic field or an electric current vía spin torque (chapter 4), 
leading to a change in resistance (figure 12 and corresponding text). 

Conversely, for IEC, dy is varied, leading directly to a change in the relative 
orientation of the magnetization of the ferromagnetic layers and subsequently to a 
change in resistance (figure 13). Indeed, in the case of the P state (figure 13, top 
right), minority moment(|)-spin(f)) electrons can be trapped between the two fer- 
romagnets due to constructive scattering near the F!/N and N/F' interfaces. It 
results in a one-dimensional quantum well behavior. 

For a simplified explanation, consider quantum well states (QWS) of quantized 
energy levels 


Wer hrr 


2m, 2m,d^, 


En (16) 
where Kk, is the transverse (perpendicular-to-plane) wavevector and m is an integer. 
“For increasing dy, the levels move downward. When a level crosses the Fermi 
energy £p, the QWS is populated, and the total energy increases. When the QWS 
level moves further below ep, the energy again decreases until the next level 
approaches ep" — quoted from P. Grünberg et al., Metallic Multilayers, in Y. Xu 
et al. (eds) Handbook of Spintronics, Springer (2015). Consequently, for P alignment, 
the energy oscillates with dy. 

In the case of the AP state (figure 13, top left), no such quantization effect occurs 
because, for a given spin-electron channel, spin-dependent scattering at interfaces is 
not active in a collaborative manner at both interfaces. Therefore, electrons have a 
*classical behavior and the energy of the AP state is independent of dy. 

Because the energy of the P state oscillates around the energy of the AP state, 
and in order to minimize energy, the alignment switches between P and AP states 
when dy increases, meaning that the sign of the coupling oscillates between positive 
and negative. These oscillations result in subsequent oscillations of the CIP-GMR 
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Fic. 13 — (Top) Illustration of the interlayer exchange coupling (IEC). Adapted with per- 
mission from Springer Nature: P. Grünberg et aL, Metallic Multilayers, in Y. Xu et al. 
(eds) Handbook of Spintronics, Springer, Dordrecht (2015). Copyright 2015. (Bottom) 
Experimental data showing the oscillatory behavior of the GMR as a function of the 
non-magnetic layer thickness in a ferromagnet/non-magnet/ferromagnet trilayer, due to the 
oscillatory sign of IEC. Adapted with permission from the American Physical Society: S. S. P. 
Parkin, Phys. Rev. Lett. 66, 2152 (1991). Copyright 1991. 


between a highly resistive AP state and a less-resistive P state (figure 13, bottom). 
Note that, the present reasoning for IEC does not require a current to be applied in 
the trilayer. The application of a current will only allow measuring the resistance of 
the trilayer. It will not have any effect on the relative alignment of the magnetiza- 
tions of the layers at first order, unless transfer of angular momentum takes place, 
which will be the object of chapter 4. 

To describe real systems beyond the simplified picture presented above, refined 
models must take into account the reshuffling of the spacer layer DOS due to the 
change from the initial behaviour of the bulk to that of the quantum well [M. D. Stiles, 
Phys. Rev. B 48, 7238 (1993); P. Bruno, Phys. Rev. B 52, 411 (1995)]. This change is 
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dependent on dy, as the properties of the quantum well are determined by d y. It can be 
shown that electrons have to be taken out or added in at the Fermi level, alternatively 
depending on dy. Therefore, the total energy on filling the states up to the Fermi level 
oscillates with dy. The periods of the oscillations were calculated to be 7/kr es, with 
kF es some extremal spanning vectors of the Fermi surface of the spacer layer. We note 
that, most of the multiple oscillations due to kp cancel out on integration over the 
Fermi surface, leaving only the oscillations due to the extrema, kr es. 


E Numbers 


The purpose of this part is to provide some orders of magnitude for the effects 
discussed above, in the context of the applications in which they are at play. Because 
the present book is focused on basic effects, readers may complement their knowledge 
on spintronic applications by consulting for example E. du Tremolet de Lacheisserie 
et al. (eds) Magnetism: II. Materials and Applications, EDP Sciences (2002). 
Initial GMR values were around few percent (figure 14, left), but typical best 
values for GMR now reach few tens of percent at 300K. This is mostly due to the 
different mechanisms and related length scales involved, namely, the spin-dependent 


mean free path 410) for CIP-GMR ($2.4) and the spin-diffusion length Ik, for 
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Fic. 14 — (Left) Experimental data of the GMR effect in an Fe/Cr/Fe trilayer with in-plane 
anisotropy compared to data of the AMR effect for an Fe monolayer of equivalent total 
thickness. Reprinted with permission from the American Physical Society: G. Binasch et al., 
Phys. Rev. B 39, 4828 (1989). Copyright 1989. See also M. N. Baibich et al., Phys. Rev. Lett. 
61, 2472 (1988) for related results in [Fe/Cr] multilayers with out-of-plane anisotropy. (Right) 
Experimental data of the magnetic behavior and electric GMR effect in a spin-valve structure, 
wherein one of the ferromagnetic layers is pinned by exchange bias coupling to an antifer- 
romagnet. Reprinted with permission from Elsevier: B. Dieny, J. Magn. Magn. Mater. 136, 
335 (1994). Copyright 1994. 
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CPP-GMR (83.2). Larger values of few hundred of percent can be obtained in TMR 
devices, relying on more selective orbital-dependent tunneling. The other commonly 
encountered and historically much older anisotropic magnetoresistance effect known 
as AMR comes from spin-orbit scattering and gives magnitudes of few percent. The 
impact of spin-orbit interactions and the AMR effect will be discussed in the next 
section (§2.6). 

Note that, the application of those magnetoresistance effects was made possible, 
thanks to the invention of spin-valves (figure 14, right). One of the ferromagnetic 
layers is used as a reference, through pinning by exchange bias coupling to an 
antiferromagnet [J. Nogues et al., J. Magn. Magn. Mater. 192, 203 (1999)], and the 
other ferromagnetic layer is used for sensing, as it is free to rotate when subject to an 
external field, e.g., stray fields from encoded bits in a hard disk drive for data storage 
(figure 15, top). The discovery and deep understanding of those basic effects, along 
with thorough engineering, ended up, for example, in the fabrication and use of 
several billions of sensors. Coupled to the spin transfer effects (chapter 4) for writing, 
TMR is also in play for reading in some of the latest memory (MRAM) applications 
(figure 15, bottom). 
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Fic. 15 — (Top) Schematics of a magnetic disk with a read head on top for use in hard disk 
drives. (Bottom) Schematics of a magnetic memory made of STT-MRAM cells. Reprinted 
with permission from Springer Nature: S. S. P. Parkin et al., Nat. Nanotech. 10, 195 (2015). 
Copyright 2015. 


2.6 Spin Mixing and Spin-Orbit Interactions - AMR 


Up to now, we considered that the majority moment(T)-spin()) electrons and 
minority moment(|)-spin(f) electrons flow in channels that are strictly isolated 
from each other, without mixing. We will next consider mixing due to spin-orbit 
interactions. 
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Fic. 16 — Illustration of how spin mixing, due to the spin-orbit interaction, alters the equiv- 
alent circuit representation of asymmetric electron-scattering, due to spin-dependent DOS. 


At this stage, it is important to understand that the terminology ‘mixing’ means 
a mixing of eigenstates, where the electronic wavefunction becomes a linear com- 
bination of the |f} and ||) eigenstates — here due to spin-orbit interaction (see also 
$4.1 for a toy model describing mixing in the context of a transfer of angular 
momentum). We do insist on the fact that no spin-flip scattering needs to be con- 
sidered at this stage. The spin-flip scattering will be tackled in chapter 3. 

How spin mixing alters the equivalent circuit representation of asymmetric elec- 
tron scattering is shown in figure 16 (vs. figure 10), where the resistivity p!! accounts 
for spin mixing. The equivalent resistivity of a ferromagnet is therefore given by 


_ pla + pP (gl +p") 
p! + p! +4pt! 


p (17) 

Note that, in the absence of spin mixing, equation (17) naturally becomes 
equivalent to equation (9). 

The term spin mixing refers interchangeably to spin and moment mixing. 

The above spin mixing contribution p!! can be due to spin-orbit interaction, 
which couples the spin (8) and orbital (7) angular momenta of an electron, i.e., the 
electron’s spin and its motion. In 3d transition metals, spin-orbit interaction causes 
two effects: 


E lifts the energy states of the majority moment(T)-spin(|)) electrons and 
minority moment(|)—spin({}) electrons, and 


Bi mixes or reorients the d-orbitals. For 3d electrons, with principal quantum 
number n — 3, the orbital quantum number and its projections along the 
quantization axis are | = 2 and m; = —2, —1,0, +1, +2 (figure 17). 


In a quantum mechanical description, the corresponding Hamiltonian is 
H = Ego L.S (18) 


where £so is the spin-orbit coupling energy, with 
1 
LS = L,S,-- LySy+ L,5, = LaS; + zL S +L S*) (19) 


where L and S are the orbital and spin operators, and L* ^ and S*© are the 
corresponding ladder operators. These ladder operators involve the orbital and spin 
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Fic. 17 — Electron probability densities for the d-orbitals in solid. The color-code refers to the 
sign of the wavefunction: blue for negative and red for positive. From P. Flowers et al., 
Chemistry, OpenStax, Houston (2015) - https://openstax.org/books/chemistry / pages/ 
l-introduction — CC BY License. 


quantum numbers (l and s) and their respective projection quantum numbers 
(m; and ms) as follows 


L*O|L m) = /I(14-1) = mimi + (-)1)A|L mi + (—)1) 
S* ?|s, ms) = hls, m, + (—)1) 


(20) 


where |}, m;) and |s, ms} are the orbital and spin states. 


Consequently, a 3d! (mj) state will mix with a 3d!” (m; + 1) state by spin-orbit 
interaction. To set orders of magnitudes, sg can reach up ~10meV in some 
3d-metals, to be compared with band-spittings of a few eV, and Fermi energies of a 
few tenth of eV (1eV = 1.602 x 10 P J). 

A simplistic illustration of mixing of d-bands due to spin-orbit interactions 
considers that the majority moment(T)-spin(l)) electrons band acquires a minority 
moment(|)-spin(f) character and vice versa (figure 18). This is a practical way to 
describe the fact that wavefunctions become a mix of eigenstates. 

Based on the above description of spin mixing, we now describe the anisotropic 
magnetoresistance effect (AMR). For the sake of simplicity, we now consider a 
‘strong’ ferromagnet. This term is used when the Fermi level crosses only one of the 
3d-bands (see Ni and Co in figure 7). It is opposed to the term ‘weak’ ferromagnet 
when both bands are crossed at the Fermi level (see Fe in figure 7). Thus for a strong 
ferromagnet, the majority moment(1)-spin(4)) electrons band is full. It is below the 
Fermi energy and no d-state is a priori allowed for scattering of a T state: 
N t (s F) =0. 

Without spin-orbit interaction, when égo = 0 (figure 19, left), the s! to d! 
scattering channel is shut due to strong ferromagnetism, and the following three 
scattering channels are open: s! to sl, s! to sl, and s! to d!. 

With spin-orbit interaction, when €s9 £0 (figure 19, right), the mixing of 
d-bands modifies the spin-dependent scattering rates, mostly because it opens a new 
channel for scattering. Here, a s! to d! scattering channel opens up because the 
minority moment(|)—spin({}) electrons band acquires a majority moment(T)-spin(l)) 
character. Consequently, the resistivity and the total resistance of the ferromagnet 
are changed, here it becomes larger. 
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Fic. 18 — Illustration of how spin mixing, due to spin-orbit interaction (&5¢), alters the DOS 
of d-electrons. The illustration is simplified for the sake of clarity and does not mean that there 
is a loss of electrons. Adapted from S. Kokado et al., J. Phys. Soc. Jpn. 81, 024705 (2012) — 
CC BY license. 


Fic. 19 — (Top) Illustration of how spin mixing, due to spin-orbit interaction (£56), alters sd 
electron scattering by opening a scattering channel — with the example of a strong ferro- 
magnet, like Ni. (Bottom) Corresponding equivalent circuits. 


A practical example of the effect of spin-orbit interaction on a ferromagnet’s 
resistance is when the direction of an applied current I — J.A, with A the 
cross-section surface, is varied with respect to the orientation of magnetization 
M. The result of spin-orbit interactions, indeed, depends on the direction of the 
angular momentum of s-electrons (parallel to the wavevector hk) carrying the cur- 
rent relative to that of the delectrons (parallel to M) bearing magnetization. 
A simple picture is given in figure 20: 
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Fic. 20 — Representation of the anisotropic magnetoresistance (AMR) effect, which is a 
prototypical example of how spin mixing, due to spin-orbit interaction, influence electron 
transport. Here, spin mixing depends on the orientation between current (1 — linked to the 
angular momentum of s-electrons) and magnetization (M - linked to the angular momentum 
of d-electrons). It is, for example, adjusted through the application of a magnetic field. 


WB when I || M, the electronic orbits are L to I, offering a large cross-section for 
the selectrons to scatter, resulting in a high resistivity pj. This is equivalent 
to the case depicted in figure 19, right. 

E when 7 L M, the electronic orbits are || to J, and thus offer a smaller 
cross-section for scattering, resulting in a low resistivity p,. This is equivalent 
to the case depicted in figure 19, left. 


This anisotropic magnetoresistance (AMR) effect was demonstrated experi- 
mentally by W. Thomson (Lord Kelvin) in 1857 and explained theoretically in the 
1970ies by Campbell, Fert and Jaoul (CFJ model) [I. A. Campbell et al., J. Phys. C: 


Solid State Phys. 3, S95 (1970)]. AMR is expressed as follows 
^ a 
Ap _ py PL (21) 
p PL 


Exercise 6.1 is devoted to AMR and, in particular, to the calculation of its 
angular dependence. 

Some metals’ typical resistivities and AMR ratios are given in figure 21 (see also 
figure 14 for comparison with GMR magnitudes). The composition-dependence of 
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Fic. 21 — Experimentally determined values of AMR for several compositions of NiCo and 
NiFe alloys. Reprinted with permission from IEEE: T. McGuire and R. Potter, IEEE Trans. 
Magn. 11, 1018 (1975). Copyright 1975. 
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AMR partly relates to band filling, i.e., to how the Fermi level is moved with respect 
to the densities of states, which will alter scattering probability. 

Note that, the presence of a magnetic texture, like a domain wall (DW), may 
change the relative orientation between current and magnetization compared to the 
uniform state, and subsequently alter the resistivity of the system. This effect is 
known as the DWAMR and it will be evoked when discussing the general context of 
DWMR [equation (43) and corresponding text]. In addition, exercise 6.2 is devoted 
to DWAMR. 


Summary 


In this chapter, we have seen that spin-dependent transport in band ferromagnets 
are due to an asymmetric scattering probability due to spin-splitting of the 
majority moment({)-spin(J)) electrons band relative to the minority moment(|)— 
spin({t) electrons band, N'(er) 4 N'(er). This asymmetry leads to a 
spin-dependent mean free path AK and related Ohm's law j! = gl Vu’) fe, Tt 
is formulated in terms of the basic spin asymmetry parameter « = p! /p! , and is at 
the origin of the CIP-GMR effect. Adding spin-orbit interactions to the problem 
generates spin mixing o!! which can be thought of as a mixing of the majority 
moment(1)-spin(|J) and minority moment(|)-spin(f)) electrons bands, therefore 
reshuffling the spin-dependent scattering probabilities. The spin-orbit-related 
spin mixing mechanism is at the origin of the AMR effect. 


Chapter 3 


Spin Accumulation - CPP-GMR 


In chapter 2, we have seen that longitudinal transport is spin-dependent in ferro- 
magnets. It is formulated in terms of the basic asymmetry parameter «= p! /p!, 
which is due to distinct spin-dependent DOS at the Fermi level, leading to a 
spin-dependent mean free path, and thus causing the current density to be different 
for the majority moment(t)-spin(l) electrons and minority moment(\)—spin(f) 
electrons: j! Z j'. As a result, when electrons flow across an interface (figure 22), the 
distinct partial densities of current at the interface between materials of different 
types with different asymmetry parameters, for example, a ferromagnet and a 
non-magnet create a spin imbalance. This effect is called spin accumulation. 
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Fic. 22 — Schematics of a ferromagnet /non-magnet (F/N) metallic bilayer subjected to a 
perpendicular-to-plane charge current, illustrating the creation of a zone of spin accumulation 
at the interface, which accommodates for the difference in spin asymmetry (§2.4) between the 
two layers. 


How electron dynamics is described in such a case; which types of relaxation 
towards the steady state regime prevail; which conditions are necessary to maintain 
spin polarization across the interface; and how spin accumulation is used in spin- 
tronic devices, such as CPP-GMR sensors, will be detailed in the next sections. 
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3.1 Spin-Dependent Spin-Flip Scattering 


We first introduce the spin-flip scattering relaxation process (figure 23), by which 
the spin is eventually flipped upon scattering, after about a hundred non-spin-flip 
events. Spin-flip scattering makes the majority moment(t)-spin({) electrons 
channel communicate with the minority moment(\)—spin(f) electrons channel and 
vice-versa. 


Fic. 23 — Representation of electron scattering when a spin-flip scattering event is consid- 
ered. 4, is the (spin-dependent) electron mean free path. A, is the (spin-dependent) spin-flip 
length (total length of the red path). l is the (spin-dependent) spin-diffusion length; it is in a 
way A, reduced to a given direction, which is set by an applied current in a real 
three-dimensional experimental setting. Note: to make the representation easy to read, 4, and 
A, are shown between particular collisions. In practice, 4, and 4,; are averaged over all 
collisions. 


In addition to the electron mean free path Ae = vp, (82.2), a spin-flip length can 
be defined 4;; = vrtat. To avoid misleading the reader, we want to insist on the fact 
that, in figure 23, the spin-flip length corresponds to the length of the path colored 
in red. What matters for CPP transport, is the direction perpendicular to the 
interface. 4,; is thus barely used and another length called the spin diffusion length 
lf, is more relevant. Indeed, it describes the length between two spin-flip scattering 
events, along a given direction. Based on Einstein's relation, we have lẹ = 4/ Dis. 
Based on Fick's law, the diffusion constant D is written D = iv rA, and the 


spin-diffusion length can also be expressed as /, = 4/A,4,/3. This expression 
underlines that lf is the geometrical mean between 2e and A. 

Recall that lẹ > A,, which has consequences for the GMR in CIP vs. CPP 
geometries [$2.5, figure 25 and J. Bass and W. P. Pratt Jr, J. Phys.: Condens. 
Matter 19, 183201 (2007) for numbers]. In some cases, spin-flip events can be as 
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scarce as 0.196. It is also important to note that, because the mean free path is 
spin-dependent, the spin-diffusion length will naturally be spin-dependent. In 
addition, we admit here that the spin-diffusion length is also spin-dependent because 
spin-flip events can be spin-dependent. The spin-dependent spin-diffusion length 


therefore reads p" DO B" (two equations). 


We point out that, f and 4 represent the time for a majority moment()-spin 
(4.) electron and a minority moment(\)-spin(T) electron to be flipped, respectively. 
The notations t’! and t!! are sometimes used in the literature for spin-flip scat- 
tering. However, these latter notations can be misleading. In this book, the symbol 
Tl (page 151) is exclusively used to describe spin-mixing processes, i.e., mixing of 
the eigenstates, in the absence of spin-flipping (82.6). 


3.2 Drift-Diffusion Model, Spin-Coupled Resistance 
& Mismatch 


In the previous sections, we introduced the idea that the distinct partial densities of 
current induced at the interface between two materials of different natures create a 
spin accumulation (figure 22). Relaxation towards equilibrium conditions causes 
spins to diffuse near the interface. The average spin-diffusion lengths on either side of 
the interface D. rcv) are the characteristic lengths for this relaxation process. How 


uu 
F(N) 


[equation (27)]. We remind the reader that s notation Mr r(x) accounts for two 
n 


Up F(N) come into play and how they connect to ly will be detailed in this section 


lengths: [%, sr in the F layer and [ ssy In the N layer; and that » yy thus accounts for 


four lengths: Il, p, 15, ps Uy y, and I], y. 
The purpose of what follows is to determine, in the whole F/N stack: 


W the spin 5 accumulation or chemical potential imbalance: u, = — (u! — u^), 
WB the spin s current: J; = -(j —j'). 


Next, we show all the steps that allow us to find the expressions of u, and Js as a 
function of the layers’ parameters. At this stage, we skip the calculations, in order to 
go straight to a physical understanding of the phenomenon. Detailed calculations 
are given in exercise 6.3 and in P. C. von Son et al., Phys. Rev. Lett. 58, 2271 (1987), 
and T. Valet and A. Fert, Phys. Rev. B 48, 7099 (1993). 

In the literature, spin accumulation is defined either as u, = pm! — ul, or 
i= —(u! E u^), Or Ji, = x(u E ut) /2, resulting in a factor +1/2 difference in some 
of the equations, depending on the articles and books consulted. Throughout this 
book, the term spin refers to 8 and the spin 8 accumulation is thus defined as 
u, = pt — pl = - (ul = ut), because the T and | symbols refer to the orientation of 
the moment m, whereas the ¢} and |J symbols refer to the orientation of the spin 8, 
which points in the opposite direction to m (see symbols and units page 150). 
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Bi Step 1: generalized Ohm's law 


We first recall the generalized Ohm's law (4 equations) 


val 
AQ) _ LG) F(N) 
Jen) = PFW) ¢ (22) 


2 pl) 
with um ^ = uv N) 7 Ung + pA ni Drift, described by the electrostatic 
potential u, (82.2), is omitted here to facilitate reading. It will cause straightforward 
contributions to add up. In fact, in the most general case, u is defined as u = He + My 
and is called the electrochemical potential. 14, is called the chemical potential (linked 
to the diffusion process). It corresponds to the energy necessary to add one particle 
(here, electrons) to the system and thus equals the Fermi energy at equilibrium and 
zero temperature. It is defined as 4, = Hp 4 Na = fly, d £n on, where Hp, is the 


equilibrium chemical potential and ón is the electron density in excess compared to 
the steady state in the isolated layer. 

It is important to note that u„ is spin-dependent because o and ôn are 
spin-dependent parameters. It is also important to understand that, in the absence 
of spin-accumulation, the excess particle density as defined above is equal to zero, 
and Hp = Hm: So, in a single layer at steady-state equilibrium, “= ut, + Hn, and 
Vu = Vu, as used in chapter 2, see e.g. equation (7). 

Obtaining equation (22) and the equivalence to spin-dependent macroscopic 
transport equations is not as straightforward as it may seem. The accurate calcu- 
lation uses Boltzmann transport theory How the Boltzmann equation model 
reduces to macroscopic transport equations when Uy > 2e is extensively described 
by T. Valet and A. Fert in Phys. Rev. B 48, 7099 (1993). 

In the following, the F and N subscripts will be omitted to facilitate reading. The 
equations and steps discussed below, apply in the F and in the N layer separately. 


Bi Step 2: charge and spin currents 


Recall that the definitions of the charge and spin current densities [equations (3) and 
(4)] are 


J.—j tf (23) 


Js= -(j' - 5) (24) 


W Step 3: charge conservation 
The flux of charge current density is conservative 


V-Je=0 (25) 
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E Step 4: total angular momentum conservation 


The flux of angular momentum current density Q — EJ, is not conservative. 
However, total angular momentum is conserved. As stated by Lavoisier: “Nothing is 
lost, nothing is created, everything is transformed", meaning that any loss in spin 
current must be due to spin-flips (at this stage). This is expressed as follows 


SnI l 
V-J,= (S x (26) 


l 
Jum Ti 


At this stage, we only consider spin-flip scattering as the process by which spin 
current is lost. There are several other sources of spin angular momentum relax- 
ation, e.g., electron-phonon, electron-magnon, and electron-electron interactions, 
besides the spin-flip scattering process discussed here, but also the spin mixing 
process discussed in $2.4 and sd transfer of angular momentum through the spin 
torques described in chapter 4. 


Bi Step 5: obtain the diffusion equation for spin accumulation 


Using equations (22)-(26), it is possible to show (exercise 6.3) that CPP transport at 
a heterogeneous interface is governed by a simple diffusion equation describing spin 


accumulation (or chemical potential imbalance) u, = —(u! — ut) 
Vi, - I 0 (27) 
Uy 
where I5, can be viewed as an average spin-diffusion length, because = = 3 + na In 
of sf sf 


the literature, Ur is most often written as ly. The notation lj; may implicitly infer 
that the spin-diffusion length is spin-independent, which is not always the case. 
Hence, we prefer to use a superscript as a warning. 

When a time-dependent chemical potential a (§6.6) and external magnetic field 
H come into play, it is possible to show that equation (27) becomes 


On, 


i 2, — Hs 
ðt — 7| Us x (uH) T DV Hy Taf (28) 


where p, is a vector ($3.5 and $6.6) and —|y|u, x (ugH) accounts for the possibility of 
spin accumulation to precess around the magnetic field. The gyromagnetic ratio is 
y=-—»5,-9, and the Landé g-factor is g~2 for an isolated electron, making 


bl 28 GHz.T-!. 

Other effects may influence spin accumulation. An example of how the 
spin-dependent Seebeck effect (§5.5) alters the drift-diffusion equation for spin 
accumulation is shown in exercise 6.3. 


E Step 6: solve the diffusion equation for spin accumulation 


In a one-dimensional problem, in which the current flows along z, the solutions of 
equation (27) take the form (two equations) 
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Hs py (2) = Arn) el ron + Brn) ere) (29) 


where Ary) and Bron) are four constants, obtained from the boundary conditions. 
For infinite layers, we have 


im (usu) = 0 (30) 


z—-—(-)oo 


We note that exercise 6.6 proposes a problem in which the diffusion equation for 
spin accumulation is solved for finite-size layers. 
Continuity of u, (due to continuity of ju!) ) and J, across the interface gives 


Is, p (0) = ns (0) (31) 


Js. r(0) = sn (0) (32) 


Using equations (27)-(32), it is possible to calculate the spin accumulation in the 
F/N bilayer (figure 24, see T. Valet and A. Fert, Phys. Rev. B 48, 7099 (1993) for 
detailed calculations) 


-1 
i. Np tO 

Hoy (2) = 26] JL, vp (: dm e wm (33) 
sf, FPF 


with pp = pr/ (1 — p^). We recall that f — PERI is the polarization of the current in 


| 
the F layer, which is linked to the spin asymmetry parameter « = "F due to the 
F 


spin-dependent mean free path (§2.4). 
From figure 24, middle and equation (33), we see that the average spin-diffusion 
lengths on either side of the interface Uy, F(N) are the characteristic lengths for spin 


accumulation, as anticipated in the introduction of this chapter (see also figure 25 
for numbers). 

Note that u, py) is proportional to Je, meaning that increasing the current in a 
device will increase spin accumulation, though at the expanse of Joule heating and 
power consumption. In addition, because the spin accumulation is proportional to the 
charge current, it is odd in J, and consequently changes signs under reversal of the 
current direction, jt, p(y) (Je) = —Ms,p(w)(—Je)- This also implies that the spin accu- 
mulation changes sign under the reversal of the stacking order, 
Us p CP N) = — n, pay (N/F). These properties will come into play in the next 
section (83.3), when discussing F/N/F heterostructures in P and AP states. Typical 
values of charge currents used to estimate the effects related to spin accumulation are: 
Je~ 10? — 105 A.cm"?, and typical magnitudes of spin accumulation are 
LH, ^» 10 — 100 peV. 

Note also that, u, p(y) is proportional to f, reminding us that efficient spin 
accumulation is obtained by selecting appropriate materials, especially by selecting 
a F layer with a large spin asymmetry parameter, o (see figure 25 for numbers). Spin 
accumulation is thus odd in f and then changes sign under magnetization reversal, 
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Fic. 24 — (Top) Schematics of a ferromagnet /non-magnet bilayer, illustrating the zone of spin 
accumulation at the interface. (Middle) Corresponding spin 8 accumulation u, = — (ud = p), 
calculated for two orientations of magnetization M. The | and | symbols refer to the ori- 
entation of the moment 4n. Drift is omitted to facilitate understanding. (Bottom) Corre- 
sponding ratio between spin 8 current Js = -( — J) and charge current Je =j! +jt, 
calculated for the case of metallic and semiconducting non-magnets. Adapted with permission 
from EDP Sciences: A. Fert, Reflets Phys. 15, 5 (2009). 


as illustrated in figure 24, middle, u, p(y) (M) = ~us ron) (7 M), or written differently 
to ease further understanding, H, pw (F!) = -Hs rxy) (F^). This property will also 
come into play in §3.3. 
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Finally, note that u, p(y), which refers to the spin 8 accumulation, is negative for 
Je «0 (charge current going from N to F) and fj > 0 (magnetization ‘pointing up’), 
as the majority moment(T)-spin(l)) electrons are then in excess at the interface 
[figure 24 and equation (33)] and u, = ut — wh = —(u! — ut). 

Further using equations (22)-(32), it is possible to express the spin current in the 
F/N bilayer (figure 24, bottom) 


oe * = T 
Js r(x) = BJe 1 (: } att) er (34) 
sf, N 
E meee 
Js,.n(2) = BJ. ( | a ey. (35) 
sf, FF 


From figure 24, bottom and equations (34) and (35), it becomes clear that, as the 
result of spin accumulation, spin polarization diffuses from the F layer to the N layer. 
It is also said that spin is injected in the N layer. The same remarks as above can be 
drawn on the odd symmetry of J; roy) with Je and f. 


Bulk resistivity Bulk asymmetry Spin diffusion length 
reme B È [nm] 
Co 15 0.35 38 
Cu 2 0 100 
Ru 20 0 14 
CoggFe,o 20 0.7 25 
Nig,Fejs 20 0.6 4.5 
Cug4Pte 35 0 7 
Interfacial resistance Interfacial spin Interfacial spin memory loss 
r[mQ um?] asymmetry y ô 
Co/Cu 45 0.7 25% 
Co/Ru 48 -0.2 
CogpFe jo/Cu 45 0.75 25% 
Coog9Fejo/Ru 48 -0.2 31% 
Cu/Ru 48 0 35% 
NisgFeso/Cu 25 0.63 33% 


Fic. 25 — Examples of the bulk and interface (§3.4) parameters at stake in spintronic effects 
involving spin accumulation. Note: r, y and 6 are parameters defined in the next paragraph. 
They take into account spin-dependent process due to the interface itself. Reprinted with 
permission from the American Physical Society: A. Manchon et al., Phys. Rev. B 13, 184418 
(2006). Copyright 2006. See also J. Bass and W. P. Pratt Jr, J. Phys.: Condens. Matter 19, 
183201 (2007) for more numbers. 


W Spin resistance 


In the same manner, as p (c) is called the resistivity (conductivity) because it 
accounts for any process that opposes (favors) the flow of electrons, pli; (c/15.) is 
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naturally called the spin resistance/impedance (spin conductance) as it opposes 
(favors) the flow of spins. This notion is commonly encountered in the field of 
spintronics, especially when dealing with spin injection and transmission, as we are 
here. Note however that, it is customary to incorrectly use the term spin ‘resistance’ 
and ‘conductance’, although these quantities refer to a resistance area product in 
Q.m? and a conductance per area in S.m ?, respectively. In practice, spin conduc- 


. ofl, E 
tance is often used per quantum conductance Ga in m~?. It can be compared 


directly to the spin mixing ‘conductance’ g!! = ip (§4.2). 
Bi Spin impedance mismatch 


To further underline the importance of the choice of the layers’ parameters for 
efficient spin injection, we remark from equations (34) and (35) that, when 
Gs NPN > Uy ppp, for example, in the case of an F-metal/N-semiconductor bilayer 
(figure 24, bottom), spins accumulate and depolarize in the F layer. No spin 
polarization diffuses in the N layer, as J,w = 0. This issue is known as the spin 
impedance mismatch. Overcoming this matter requires adding a large 
spin-dependent resistance at the interface, e.g., a tunnel barrier. More about 
impedance mismatch is explained in exercise 6.4, based on A. Fert and H. Jaffrés, 
Phys. Rev. B 64, 184420 (2001). 


E Discontinuity of the potential and corresponding spin-coupled interface 
resistance 


As a final remark on this part, we discuss the notion of a discontinuity of the 
potential at the interface, often encountered in papers. First, to obtain the elec- 


trochemical potential I) the (drift) contribution of the electrostatic potential 
PEEM UH p(N) haereat enced 


He (N): Which we initially omitted to facilitate reading (see step 1, page 34), shall be 


added to the chemical potential ‘From equation (7) and the corresponding 
CHOC Mn F(N) 


text, it naturally follows that the electronic potential is linear with z: 
He pw) (X) = eJep rent. The plots of the electrochemical potentials for majority 
moment(T)-spin(|)) electrons and minority moment(|)-spin(f)) electrons are given 
in figure 26. 


We note that the out-of-equilibrium chemical contributions pio 


trochemical potentials Hi = He, F( EST N) vanish far from the interface. 


(N) to the elec- 


Indeed, in this case, the F and N layers behave as if they were isolated, when there is 
no excess density ôn of majority moment(T)-spin(4.) electrons and minority moment 
(1)-spin(f)) electrons. 

From figure 26, a simplification consists in considering only the linear depen- 
dence of the electrostatic potentials u, p(y) on either side of the interface, and how 
they mismatch at the interface (Au,) when extrapolated back there. In fact, this 
mismatch is due to the non-monotonic contributions of the chemical potentials 


BA N) subsequent to spin accumulation. These latter continuous contributions of 
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Fic. 26 — Illustrations showing how spin accumulation u, = — (ud — u^) can be reduced to a 


discontinuity —Aj, of the electronic potential, and henceforth, to a simple spin-coupled 
interface resistance area product r, = —Ay,/(eJ-) and corresponding voltage drop. 


T N) [equation (31)] which extend over finite sizes defined by the spin diffusion 


lengths Uy Fc) are thus relegated to a discontinuity of electrostatic potential Au, and 


subsequently to a simple voltage drop. It can be expressed in terms of a spin-coupled 
interface ‘resistance’ (in fact a resistance area product in Q.m?) 


e —An, o Ls i (0) 2 us n(0) 


s 36 
" eJ, eJ, eJ, (36) 
Using equation (33), one obtains 
j B 
Ts = 2Bly NPN | 1+ m (37) 
ET 


'This spin-coupled interface resistance can also be reworked as follows 
Ts = Bly NPN (38) 


where f* is the effective spin polarization at the interface. 


3.9 Heterostructures - CPP-GMR 


In the previous section, we studied how a CPP current creates spin accumulation at 
an F/N interface. We discussed which parameters are essential — considering 
spin-flip scattering only - and calculated the  thickness-dependence of 
spin-accumulation, as driven by a macroscopic (drift-)diffusion equation. Let us now 
consider the case of a CPP current in an F/N/F trilayer (figure 27), in which spin 
accumulation: (i) occurs at the F/N and N/F interfaces and (ii) depends on the 
relative orientation of the F layers. 

In the following, when commenting on figure 27, we will focus on the central part 
consisting of the F/N/F trilayer and discard the outside parts, for which periodic 
boundary conditions were set in order to ease calculations. 
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In the F! / N/ F! P state (figure 27, left), for Je > 0, the spin accumulation p, is 
positive at the F'/N interface and negative at the N/F' interface. It is a conse- 
quence of the odd symmetry of p, with stacking order: u, ~y)(F'/N) = 
—Hs F(N) (N/F') (page 36). 

Conversely, in the F! / N/ F* AP state (figure 27, right), the spin accumulation u, 
is positive at the F!/N interface and positive at the N/ F! interface. It is a conse- 
quence of the odd symmetry of u, with stacking order combined with the odd 
symmetry of Hs with magnetization reversal (page 36): 


Us ren (FT / N) = —(—)us re) (N/F!) = Us pny (N / FL). 


x 


Fic. 27 — (Top) Schematics of a ferromagnet /non-magnet/ferromagnet (F/N/F) trilayer set 
in the parallel (P) and antiparallel (AP) magnetic states. (Middle) Corresponding calculation 
of the spin 8 accumulation u,. Drift is omitted to facilitate understanding. (Bottom) Corre- 
sponding charge currents for majority moment(T)-spin(|)) electrons and minority moment 
(1)-spin({}) electrons. Periodic conditions are set outside the F/N/F trilayer, to ease calcu- 
lations. Adapted with permission from the American Physical Society: T. Valet and A. Fert, 
Phys. Rev. B 48, 7099 (1993). Copyright 1993. 
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Because spin accumulation in the P and AP states differs, the resulting 
spin-coupled resistances at the interfaces also differ and so does the overall resistance 
of the stack. For small enough thicknesses dp and dy of the F and N layer, respec- 
tively, such that drwy) < hr. F(N) it is possible to show that the expression for 


CPP-GMR reads 
AP | aP 20*.d 2 
CPP-GMR =" — =f Ger 7) 5 (39) 
p (os; dw + 2p%-dp) 


'The exact derivation of CPP-GMR can for example be found in Phys. Rev. B 48, 
7099 (1993). Note that CPP-GMR vanishes for dp(y > I; piv)? especially because 
the F/N and N/F interfaces become independent and disconnected when dy > L7, y. 

Examples of how typical parameters influence CPP-GMR are given in figure 28. 
Because l Au << Uc, it follows that CPP-GMR falls much faster in CoFe(1)/Au 
(dy) /CoFe(3) than in CoFe(1)/Au(dy) /CoFe(3) (nm) stacks [figure 28c, left]. When 
dp increases in NiFe(dz)/CoFe(1)/Cu(4)/CoFe(3) and CoFe(dz)/Cu(4)/CoFe(3) 
(nm) stacks, the absolute MR increases, as more polarization builds up, and it then 
levels out when dp ~ l}; p [figure 28b, right]. Above the threshold thickness dr ~ l}; p, 
adding more film makes no further contribution to polarization, but it increases the 
CPP resistance instead, explaining the non-monotonic dependence of CPP-GMR 
with dr [figure 28c, right]. 


3.4 Interface Scattering, Spin Memory Boost, and Spin 
Memory Loss 


In the previous sections on CPP-GMR, bulk spin-dependent scattering processes 
were considered; through the bulk spin-dependent mean free path a Ny and 
spin-flip scattering processes via Uy F(N): When the resulting spin accumulation can 


be modelled by a spin-coupled interface resistance, this type of resistance is not 
related to the spin-dependent properties that are inherent to the interface itself. 

Refined models must consider interfaces as well, because, for example, deco- 
herent roughness-induced stray fields or changes in the local DOS near interfaces can 
generate spin-dependent scattering as well. 

We will now discuss two possible ways to describe the spin-dependent properties 
inherent to interfaces. These types of modelling can be combined into complete 
models. Yet other ways to account for more spin-dependent interface properties are, 
for example, detailed in V. P. Amin and M. D. Stiles, Phys. Rev. B 94, 104419 (2016) 
and ibid 104420 (2016). 


E Spin memory boost 


It is possible to introduce the notion of interfacial spin-dependent electronic scat- 
tering, and model interface-related changes in the partial DOS. To do so, the 
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CoFe(1) / Cu(d,,) / CoFe(3) NiFe(d;) / CoFe(1) / Cu(4) / CoFe(3) 
vs CoFe(1) / Au(d,) / CoFe(3) vs CoFe(d,) / Cu(4) / CoFe(3) (nm) 


0.04 99-4 


R.A (Q.um?) 
R.A (Q.um2) 


AR.A (Q.im?) 
AR.A (Q.um2) 


AR/R (96) 
ARIR (%) 


40 50 
spacer layer thickness (nm) 


Name p P r y I op 
NigogCros 140 uQ cm 0 0.5 mO um? 0.5 20 nm 
NiggFe;o 25 0.7 0.3 0.7 4.5 
CosFeso 16 0.7 0.3 0.7 15 
Cu 5 0 0.3 0.7 50 
CosoFeso 16 0.7 0.5 0.1 15 
PtsoMnnso 180 0 0.8 0 2 
NigoCt9 140 0 0.5 0 20 


Fic. 28 — Calculations illustrating how the electron transport parameters in correlation to the 
layers’ thicknesses influence: the resistance area product in the P configuration, the absolute 
MR area product, and the relative CPP-GMR, in typical F/N/F trilayers. Note: y is a 
parameter defined in the next paragraph. It takes into account spin-dependent process at the 
interface. r is an interfacial resistance area product related to y. Adapted with permission 
from AIP Publishing: N. Strelkov et al., J. App. Phys. 94, 3278 (2003). Copyright 2003. 
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interface is modelled by adding an infinitesimally thin (ferromagnetic) layer dis- 
playing a spin asymmetry parameter y (see figure 29 for numbers). Interface 
spin-dependent resistance area products are thus introduced in the model as 
orte = 2M ntrtace(L — C7) (40) 
Tinterface — “interface Y 
with T terlasa = Tinterface/ (1 EA y 


These notations naturally resemble the ones used for bulk [equation (10) and 


corresponding text]: pio = 2p,(1— (--)p). This type of additional interfacial 


spin-dependent electronic scattering can give rise to an extra polarization of the 
current. It certainly modifies the spin-coupled interface resistance [equation (37)] 
and the corresponding voltage drop [T. Valet and A. Fert, Phys. Rev. B 48, 7099 
(1993); N. Strelkov et al., J. App. Phys. 94, 3278 (2003)]. 


Metal pair "TEN: 2AR* (fm?) YeN2AR* c (fO?) 
Co/Cu 0.87 1.0 0.9 [4] 

Co/Ag 0.85 1.1 0.9 [25] 

Fe/Cr —0.7; —0.59 1.6 1 [95, 96] 

Fe/Cu 0.55 L5 1 [89] 

Py/Cu 0.7 1.0 0.7 [71, 79] 

Ni/Cu 0.3 0.36 0.1 [90] 

Py/Al 0.025 8.5 0.2 [97] 
Coo9Fe1o/Al 0.1 10.6 1 [97] 

Fe/Al 0.05 8.4 0.4 [97] 

Metal pair Aa/a (%) Opn OF Opyp2 
Co/Cu 1.8 0.33: 963. [20] 
Co/Ni 0.05 0.35 + 0.05 [99] 
CosoyFe;o/Cu 1.8 0.19 + 0.04 [98] 
Co/Ru j 0.342 9% [103] 


Fic. 29 — Interface spin asymmetry y and spin-flip 6 parameters, which come into play in 
spintronics effects involving spin-accumulation. Reprinted with permission from Springer 
Nature: J. Bass, CPP-GMR: Materials and Properties, in Y. Xu et al. (eds) Handbook of 
Spintronics, Springer, Dordrecht (2016). Copyright 2016. Note that references from the figure 
refer to the latter citation and not to references from the present book. See also figure 25. 


B Spin memory loss 


An interface can also cause memory loss, or even worse, memory loss. 

More seriously, it is also possible to introduce the notion of interfacial 
spin-dependent spin-flip scattering to account for the fact that only a fraction of the 
spin current coming from the F layer effectively reaches the N layer. This effect is 
called spin memory loss (SML). The interface is considered as a thin extra layer with 


a finite thickness dinterface possessing a finite spin-flip length l; interface (figure 30). 
ke 


sf interface 


Spin-flipping is modelled by the interface spin-flip parameter 6 = dinterface/ 
(see figure 29 for numbers). 
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2 


x 
0 interface 


Fic. 30 - Illustrations showing how spin-dependent spin-flip scattering due to the interface 
(I) itself can be introduced in models, by adding a finite-size thin layer possessing a finite 
spin-diffusion length [7 


‘sf interface* 


'The following ratio is used as a prefactor to weight data obtained from bulk-type 
calculations 


Jin = dosertec) 


J.(x = 0) e) 


Hsu = 


More information on spin memory loss is available for example in 
J. C. Rojas-Sanchez et al., Phys. Rev. Lett. 112, 106602 (2014). 


3.5 Non-Collinearity and Non-Uniformity 
(Geometric and Magnetic) 


We first recall that spin accumulation results from the presence of excess spin 
density ôn’) compared to the steady state condition in isolated layers. It is for- 
mulated in the framework of electrochemical potentials by u, = — (u! — pt), with 


wO = ue + hn 4 xs) = p, + Up EDUN gU (page 36). In uniform sys- 


tems we have considered so far, ôn and subsequently u, are scalars. In the most 
general case, because spin is a vector, dn and consequently u, are vectors 
(exercise 6.8). Hence, spin accumulation is also often described in terms of 
magnetization imbalance m,, with m, = — Upg [Nl (ep)! — N(ep)gpt ]. This type of 
description implicitly recalls that spin accumulation is a vector property, which is 
crucial when dealing with non-collinearity and non-uniformity. To give an idea of the 
orders of magnitude, the ratio between spin accumulation and magnetization is 
m;/ M ^. 107? — 1074. 

An example of non-collinearity is the case where the magnetizations of the F 
layers in an F/N/F trilayer are neither exactly P nor exactly AP — either by design or 
due to thermal noise — giving rise to the angular momentum transfer between the flux 
of angular momentum related to m, and the layer magnetizations M (chapter 4). 

'The first example of non-uniformity is the case where the geometry of the system 
is such that it creates local excess of charge current, at hot spots, and subsequently 
of spin-current (figure 31). Because the spin accumulation is described by a vector in 
the most general case, i.e., 3 components of m,, the spin current is described by a 
3 X 3 2nd order tensor, with 3 components of J; flowing along 3 directions in space, 
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F (8nm) / Cu (2nm) / F (nm) - 


Fic. 31 — Simulations of the transport properties of a magnetoresistive F/N/F nanopillar 
sandwiched between two extended non-magnetic electrodes. (Left) The color map indicates 
the amplitude of the charge current Je (scalar). The arrows indicate the flow of charges Je 
(vector). (Right) The color map indicates the amplitude of the ycomponent of spin accu- 
mulation m, (scalar). The arrows indicate the corresponding flow of yspins J} (vector), with 
J = JEx-4- J”. We recall that, an overcaret is used for unit vectors, e.g. X = x/x. Reprinted 
with permission from IEEE: N. Strelkov et al, IEEE Magn. Lett. 1, 3000304 (2010). 
Copyright 2010. 


ip g^ qe 
J,= c= [J® JO Je Jm Jw J” (42) 
E P. 


where the i and j indices in the J‘ vectors and JË scalar coefficients account for the 
spin and space components, respectively [see also equation (50) for a quantum 
mechanical treatment]. 

'The second example of non-uniformity concerns the case of magnetic textures 
like domain walls (DWs), figure 32. The DW magnetoresistance is defined as 


Ap E Pwith DW — Pwithout DW = Ppw — Po (43) 
p Pwithout DW Po 


Spin accumulation and associated resistance can exist in two ways in magnetic 
textures: outside and in the texture. 


B Outside the magnetic texture 


A naive argument would consist of considering the wall as an abrupt transition when 
the spin-diffusion length /7, is much larger than the DW width wpw. Equation (33) 
obtained for an isolated interface (83.3) is then used to describe spin accumulation 
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Jis 


electron flow 
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Fic. 32 — Schematics of a Bloch magnetic domain wall illustrating the phenomenon of 
mistracking between magnetization M and the moment density m of itinerant s-electrons, 
which contributes to domain wall magnetoresistance. 


around the wall. The corresponding spin-coupled resistance [equation (37), with 
N — F| serves to estimate the wall resistance. However, this vision is most often 
wrong. Indeed, looking at it more closely, the wall itself needs to be considered, 
especially the ability of the moment (spin) density m of the itinerant s-electrons to 
follow magnetization M (figure 32). In the ultimate case when there is perfect 
tracking, the electron spin sees a homogeneous medium, resulting in no spin 
accumulation and no wall resistance. 

In the most general case, spin-tracking reduces the ability of the wall to generate 
spin accumulation around it. We note that the term spin-tracking is equivalent to 
moment-tracking. A full description considers the spin mistracking parameter €. In 
low spin-orbit materials, this parameter can be viewed as the square of the ratio 
between the exchange length lsa = y Dh/ J;4, and the spin diffusion length ly (see 
also $4.2 and $4.3). The change in resistance area product due to spin accumulation 
around the wall is [E. Simánek and A. Rebei, Phys. Rev. B 71, 172405 (2005)] 


8/3 


Atm, = -g Bee 2. (44) 


'This contribution is usually negligible compared to other contributions. 


W In the magnetic texture 


A more subtle ‘intrinsic’ contribution to DW magnetoresistance is directly related to 
spin mistracking as this gives rise to a local and gradual mixing of eigenstates 
(figure 32), and therefore opens channels for scattering. This contribution can be of 
the order of a few percent. The corresponding change in resistance area product is 
[P. M. Levy and S. Zhang, Phys. Rev. Lett. 79, 5110 (1997)] 


12 ya 5 
Ar; == f pp "ws (1+ s - e) (45) 
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The case of large spin-orbit materials like semiconductors and the way spin-orbit 
interactions influence spin mistracking and spin transfer torque in DW (reciprocal 
effect to DW magnetoresistance) will be mentioned in §4.5. 


E Note about the AMR contribution in magnetic textures 


We remark that the spin-accumulation-induced DW magnetoresistances coexist 
with yet another contribution: the one due to AMR (DWAMR). Its origin is the 
same as the AMR effect described in §2.6 for a uniform system. It is due to spin-orbit 
interactions and depends on the angle between magnetization and current. In the 
case of the Bloch DW shown in figure 32, the DWAMR is zero, Apaypr = 0, as the 
current (along x) is perpendicular to magnetization (in the yz plane), with and 
without a DW: ppw = pọ. Exercise 6.2 is devoted to DWAMR. 


Summary 


In this chapter, we have seen that spin accumulation occurs when electrons flow 
across an interface, due to the distinct partial densities of current between 
materials of different types, for example, a ferromagnet (f) 4 0) and a non-magnet 
(B — 0). Spin accumulation is described as the difference between electrochemical 
potentials u, — — ( ul — ph). It can be calculated in the frame of the Valet and Fert 
drift-diffusion model. Relaxation to a steady state and the spatial extension of 
spin accumulation is regulated by spin-flip scattering, modeled by an average 
spin-diffusion length b. Spin accumulation gives rise to a spin-coupled interface 
resistance rs. A good match between the spin impedance plf; of the materials on 
either side of the interface is needed for efficient spin injection. The odd symmetry 
of spin accumulation with magnetization orientation and stacking order in 
heterostructures is involved in the GMR effect. Intrinsic contributions due to the 
interface (spin asymmetry and spin-flip parameters, y and 6), non-collinearity, 
and non-uniformity (geometric and magnetic) alter spin accumulation and must 
be considered for a complete description. 


Chapter 4 


Transfer of Angular Momentum — STT, 
Spin Pumping, SOT 


In chapter 3, we described how electron transport across an interface creates spin 
accumulation. In this chapter, we will detail the process by which spin angular 
momentum (or equivalently magnetic moment) can be transferred from current to 
magnetization, and the type of spin transfer torque (STT) that this effect generates. 
'The key parameters involved, such as the spin mixing conductance, will be discussed 
in detail, using a toy quantum mechanical model. How STT terms alter the trans- 
port equations, and subsequently the spin accumulation, by providing an additional 
relaxation channel for the spins (moments) (83.2) will be tackled. The ability of STT 
to trigger oscillations and magnetization reversal will be discussed, using the 
dynamic equations of magnetization. In fact, spin transfer torque STT and giant 
magnetoresistance GMR. ($3.3) are reciprocal effects. In the GMR effect, magneti- 
zation alters electron transport due to spin-dependent relaxation. In the STT effect, 
the electron transport alters magnetization for the same reasons, except that the 
current necessary to observe the effect is larger, e.g., an order of magnitude larger for 
STT-induced magnetization reversal than for GMR. Next, guidance on the mor- 
phology of STT in non-uniform magnetic textures will be given, echoing $3.5. In a 
penultimate section, we briefly present how STT is modeled in the frame of the 
magnetoelectronic circuit theory, which is widely used in current spintronics, and 
how this theory compares to the toy model presented earlier. A final section 
introduces spin-orbit torques (SOT), which result from the transfer of angular 
momentum from the lattice to magnetization, mediated by spin-orbit coupling 
between the current and the lattice. 


4.1 (sd-)Coupling, Spin Transfer Torques - STT 


Consider a Fyinnea/N/Firee trilayer, in which a current flows perpendicular to the 
interfaces. The different energy contributions in Fpinnea are such that the ‘pinned’ layer 
is not influenced by the current (the layer may have a large coercive field or should be 
exchange biased by a third-party antiferromagnetic layer). Conversely, the Free 
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magnetization is supposed to be ‘free’ to oscillate/rotate due to the action of the 
current. 

In line with the description in chapter 3, the s-electrons flowing from the pinned 
layer to the free layer acquire a spin polarization in the Fyinnea layer of magnetization 
Mpinnea and accumulate at the N/Fg,, interface as shown in figure 33. We describe 
the incoming spins by a magnetic moment density Mpinned = Mpinned / V, related to 
the spin angular momentum density Spinned = 8pinnea/ V = —Mpinnea/|}|. Naturally, 
Mpyinnca relates to the current polarization in Fyinnea and to spin accumulation at the 
N/Fg«, interface. 


Ferromagnet Non-magnet Ferromagnet 


ggg 
BPR 


M, pinned | Mypinned Msc 


electron flow 


Fic. 33 - Illustration of a pinned-ferromagnet /non-magnet/free-ferromagnet trilayer, when 
the magnetizations Mpinnea and Mg. of the ferromagnets are misaligned. The moments 
incoming on the free-ferromagnet are polarized in the pinned-ferromagnet. They are described 
by their magnetic moment density mpinneq related to the spin angular momentum density 


Spinned AS Mpinned = —|y|Spinnea- 


The incoming mhi (related to selectrons) and the magnetization M free 
(related to d-electrons) are coupled via the sd exchange interactions, which is 


expressed by the following Hamiltonian: Hg = —J sif pinnca M free, with the 

exchange constant J,;. Consequently, in the reference frame of the layer magneti- 

zation, Mfee experiences an effective field H Mico X J sdMpinncea and a torque 
l dM iro 


T» so Mirco X rmx called the spin transfer torque, STT. In a more detailed 


description, we will see that, in fact, what matters for spin transfer torque is the flux 
Of Mpinnca, Or, equivalently, up to a sign, the flux of spin angular momentum Spinned- 
Note that in the following we give the expression of —T because M goes like —|y|T 
(see also equation 65). 

Only projections of mpinnca perpendicular to M« give rise to a torque. Therefore, 
the STT is often written in terms of the following two orthogonal components: the 
damping-like (DL) torque -Tp = uM x (m x M) and the field-like (FL) torque 


Mpinned 


—Tr = ym x M, in J.m?. The terms tpi and tp, are the DL and FL magnetic 
flux density (or magnetic B-field) amplitudes, in tesla, for direct comparison with 
the action of a flux density oH oriented along (M x m) and m, respectively 
[figure 34, also equation (46). We took M = Moree and m = Mpinnea for simplifica- 
tions. Actually, the expressions in terms of the angular momentum density s — 


ym of the incoming current are also encountered in the literature: -Tp = 


ly| vM x (M x s) and —Tpy = |)| j£ M x s. Most often, these expressions are used 


Transfer of Angular Momentum — STT, Spin Pumping, SOT 51 


| 
- 
I 


JULMx(mxM) + CmxM (46) 
Mz Ms 


= ly| 2E M x (M x s) + ly| ZM x s 
M2 Ms 


aj bj 

In-plane Perpendicular 
Slonczewski Transverse 
Damping-like (DL) Field-like (FL) 


electron flow 


Fic. 34 — Representation of the damping-like (DL) and field-like (FL) contributions to 
dM /dt = —|y|T due to torques acting on a magnetization M, originating from transfer of 
angular momentum with incoming electrons of moment density m, and spin density 
s = —m/|y|. The DL torque orients M towards m, or equivalently away from s, and the FL 
torque induces a precession anticlockwise viewed from above. 


with normalized quantities. In this book, we used the non-normalized values to make 
it clear in terms of units and conversion from one notion to another. 

The above type of expressions for the STT are the most general ones, as they 
consider the universal case of polarized electrons incoming on a magnetization M, be 
the polarization of the electrons expressed as s or m, and whatever the effect at its 
origin. 

In fact, the direction of s or m depends on the effect that gave rise to spin 
polarization. Therefore, the directions of FL and DL torques can be reversed 
depending on the initial direction of m. For STT, m || Mpinnea, but in some other 
structures for spintronics, effects arising from spin-orbit coupling, like the Rashba 
effect, can give rise to a different orientation of the spin polarization and, subse- 
quently, to a different orientation of the torques. These torques are called spin-orbit 
torques (SOT) and will be introduced in $4.7. 

In articles on STT, in GMR structures, because the polarization of the current 
Myinnea relates to Mpinnea, it is common to find the torques expressed in terms of 


TE 
M pinned: —Tpy, = Tor M tree x (M pinned x Mirec) and -TFL = M M pinea x M tree- 
E f 


When more than one terminology can be found, the names damping-like 
(DL) and field-like (FL) are now accepted. The origin of such a name is direct from 
the magnetization dynamics framework [equation (65)], which will be detailed in 
$4.4. The sign convention will be detailed in $4.2. 

We note that, the total magnitude of the torque terms depends on materials 
properties (e.g., magnetic volume M sdp, and polarization fi), the geometry, and the 
current amplitude and direction (Je) [equation (65)]. About the orientation of the 
torques with regards to the current direction, the rule is as follows: when an electron 
of magnetic moment density m (and spin density s pointing in the opposite direc- 
tion) flows toward a magnetization M, the DL torque orients M towards m, or 
equivalently away from s, and the FL torque induces a precession anticlockwise 
viewed from above (figure 34). 
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Finally, we remark that, in the literature, m is used either to define a property, in 
A.m ^, units of magnetization, as throughout this textbook, or for the dimensionless 
magnetization unit vector, in which case m — o resulting in some possible 
misunderstanding and a factor M s difference is some of the equations depending on 
the articles and books consulted. In this textbook, we use an overcaret for unit 
vectors, i.e., M = in (see symbols and units page 150). 

The purpose of what follows is to determine, by use of a toy quantum mechanical 
model: 

Bl what governs the transfer of angular momentum by STT, and 


E how magnetization is influenced by the torque T. 


4.2 Toy Quantum Mechanical Model for STT, Angular 
Momentum Flow 


B Reminder 


Before presenting the quantum mechanical model, it is important to recall the 
quantum mechanical treatment of charge and spin, which we will use. 

Consider a spinor y for the full description of the spin-electron state, i.e., a 
spin-half electron wave function. 

The observable density of charge p is described by the product between the 
electron charge — e and the density probability, and the observable density of charge 
current Je additionally involves the velocity operator v 


p = —eQlv) (47) 
Je = —-eReQlv|)) = P am(QIV IV) (48) 


€ 


Regarding the spin, the density of spin angular momentum is described by 


s = (ISI) = 5 lel) (49) 


and the density of spin current by 


2 


h h 
O = J. = Re(WiS vip) = -= 


SIm((wla & V|y)) (50) 


€ 
Both expressions involve the use of the spin operator S. It is related to the Pauli 


E m PE M EUN ; [l1 ON. 
operator a by S= $a, with e^ = (1 ose 0 , and o^ = os sp 


z being the quantization axis. The expression for the spin current density involves a 
tensor product & between the spin operator S from the spin Hilbert space, and the 
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velocity operator v from the coordinate Hilbert space (see also $3.5 for the classical 
definition of the spin current tensor). 


E Toy quantum mechanical model 


We now consider an incoming spin density at the N/Feee interface, resulting from an 
initial spin polarization in the Fpinnea layer of a Fyinnea/N/F ‘ree trilayer (figure 33). 
The magnetization of the Fyinnea layer is angled by 0, whereas that of the Fg, layer 
is set along the z axis (making it the quantization axis). For the full description, 
please refer to M. D. Stiles and J. Miltat, Top. Appl. Phys. 101, 225 (2006). 
Dynamic spin transport is represented by an incoming electron-spin with a state 
|U;) — carrying a density of angular momentum s (figure 35). [V;) is a linear com- 


bination of the |f) and |J)) eigenstates: |V;) = Jr (cos? fr) + sin £ |J). We recall that 
the f and | symbols refer exclusively to the orientation of m, and the f and 4} 


symbols refer exclusively to the orientation of s. 


10 
l 


0 6 Pe iS, .. - 9 
Ii) « cosz |f) + sinz |U) g s y U el $ lpt) x cos; |f) 


0 à Ty —1 
Ir) x sin |U) Y * 7 
y X 


Fic. 35 — Schematics corresponding to a toy model for STT (see text), illustrating how an 
incoming electron with wavefunction |V;) — carrying angular momentum 8 — transfers angular 
momentum with magnetization — having angular momentum density Sy. |w,) and |w,.) are the 
transmitted and reflected wavefunctions. t} and ry are the transmission and reflection 
coefficients for the |f) and |J) eigenstates, respectively. Normalization to volume 1/4 /V and 
propagation terms e!" and e^" were omitted to simplify the scheme. 

Upon impinging on the Fpee layer, represented by the macrospin angular 
momentum VSy = -i VM, with V the volume, the incoming wavefunction is 
filtered for the |f) and |J) eigenstates (§2.4). Filtering is expressed in terms of 
transmission and reflection coefficients for each state, ty) and rya). These 
parameters take into account the sd exchange interactions between s and Sy, which 
are at the origin of the transfer of angular momentum. In this toy model, we consider 
that the |f) eigenstate is transmitted, ty) = 1(0) and the |J) eigenstate is reflected, 


rya = 0(1). The resulting transmitted and reflected states are thus expressed as 
follows: |w,) = Jy cossi) and |y,) = Fr sin|U). 
Note that, the use of ? and not of 0 comes from the fact that the eigenstates of 


opposite (x) spins are represented on an orthogonal (Z) basis. 1/V/V is used to 
normalize the volume. 
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Combining equation (49) with the incoming |y/;), transmitted |V/j) , and reflected 


(2) ) 


|W.) states, it is possible to express! the incoming s”, transmitted gO, and 


reflected gt?) densities of spins. They are decomposed along the z quantization axis 
and its orthogonal direction z by use of EUN = (Vio Wann) 


One obtains s? = Je sind and s; = É cos0 for the incoming density of spin. The 


final density of spin after interaction through transfer of angular momentum between 
the incoming electron of density of angular momentum s and the ferromagnet mac- 


rospin Sy, is the sum of the transmitted and reflected densities of spins. One obtains 


T — of m cs — Z— o z2. 165 20 Jom: 
sf = 8; +s, — 0-- 0 = 0 and s; = 8; + S; = zy COS 5 — 5p Sin 3 = gy cos. 


The loss of angular momentum for the incoming electron can be expressed as 
As? = (sf +s?) — s? = —dosinü and As? = (sj + s*) — s? = 0. Using conservation of 
total angular momentum, we obtain As-F AS, = 0. The correlated gain of angular 
momentum for the magnetic layer is thus expressed as AS m = AS), = d sind. 

This angular momentum transfer from the spin of incoming electrons to mag- 
netization can be rewritten as 


URS 
ASy = AS), x = jy nox =s; (51) 


It is said that the transverse component of spin angular momentum is absorbed. 
This eventually alters magnetization (figure 36), as the result of a torque T y, acting 


Fic. 36 — Follow-up schematics corresponding to figure 35 after transfer of angular 
momentum. The transverse component of the spin angular momentum density s, has been 
absorbed by the macrospin Sm, which has thus gained AS4,x = s1. In this example, the torque 
is damping-like. 


In this footnote, we detail the calculation of the incoming density of spin along z, as an example. 
Calculation of the other incoming 57, transmitted git?) and reflected sio) densities of spins is then 
straightforward. The Bra and Ket state vectors are represented by column vectors and their 
0 

1 ? 
as it is the use in quantum mechanics. ‘Therefore, — 57 = (y;|S"|y;) = 


conjugate transpose, respectively: (}|=(1 0), WJ)=(0 1), It) = a and |J) = ( 


211 0 0 2 


sins 
2 
and the propagation terms were omitted to facilitate calculations. Those terms can be included 
afterwards. 


0 
9 " 0 1 COS 5 . j i n 
(cos$ sin?) fi ( 2 = £ 2cos$sin$ = É sind. We note that normalization to volume 
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on the moment VSy. In turn, the spin transfer torque T y can be directly calculated? 
from the spin current flow (here along x) entering (incoming Q;) and leaving 
(reflected Q, and transmitted Q,) the volume V, and acting on an area A. It is 


: ARK op ARES (s $ ) (53) 
gg PU MX |S XƏM 


Comparing equations (53) and (46) we see that Ty is damping-like only. 
E Refined quantum mechanical model 


In refined models, other ingredients that contribute to the process of transfer of 
angular momentum are introduced (figure 37). We will now describe these ingre- 
dients and how they contribute to both a damping-like and a field-like STT. 


— Pn ot Tas 
T.i | € Tr, 
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Ferromagnetic 
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Fic. 37 — Illustration of some of the ingredients involved and introduced in models in order to 
describe STT more accurately. Those ingredients are: spin filtering due to spin-dependent 
scattering ($2.4); spin rotation and dephasing upon reflection; and spin precession upon 
transmission in the ferromagnet. When the ingredients are introduced, both damping-like 
(Tp) and field-like (Tfr) terms contribute to the torque. Reprinted with permission from 
Springer Nature: M. D. Stiles and J. Miltat, Top. Appl. Phys. 101, 225 (2006). Copyright 
2006. 


Some of the ingredients introduced to explain the physical origin of STT more 
accurately are (figure 37): 


"The density of spin current is given by Q@=—43m/((p|S-V\y)) with 
Q” = —£ Sm((y|S" - Vily) = — 2 (y|S"|y) = — *5 s", where the element Q” represents the 


flow along the direction 6 of spins polarized along a. Based on the geometry of the problem, we 
u s" s! sg 

have Q = Q = — ike | 0 0 0 | and ¥-Q= ne (s* + s! + s?)x. The torque is therefore given 
0. 0 0 


by Ty = AF - (Q;+0,- Q0) = Ík (As* + As! + As?)x = ATE: sin 3. 


V2m. 
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Wi spin filtering (82.4): spin-dependent transmission (f4(j)) and reflection (rpy)) 
coefficients. Initially, these are real numbers that are not necessarily set to 0 or 1 
like in the toy model discussed above. They satisfy tig) + "iu = 1 


E spin rotation ($) and dephasing (0$) upon reflection: e?+°%), leading to a 
mixing of the eigenstates and represented nowadays by a complex number 
parameter known as the spin mixing conductance (G^). The terminology mixing 
is naturally the same as that used in $2.6, but the origin of the mixing differs. The 
term spin mixing refers interchangeably to spin and moment mixing (G"! = G^). 


Bi spin precession upon transmission in the ferromagnet: ei h)s. influencing 
torques and spin accumulation. 


'The ingredients mentioned above are taken into account in models, in a com- 
bined manner, by considering complex number transmission tyy) and reflection ray) 
coefficients. This now enables the spin to exit the (x, z) plane. As an example, the 
calculated density of reflected spin current in this refined model [M. D. Stiles and 
J. Miltat, Top. Appl. Phys. 101, 225 (2006)] reads 


hk ; 
Xr o To: * id 
Q7 == sinRe(rirye ) (54) 
hk 
yo To. ~ * ib 
QF =- a sinüSm(ri rye ) (55) 


where Q” represents the flow along the direction b of spins polarized along a, and 
e)». 


with, by definition, rry = |rsry 


Bi Spin mixing conductance 

This latter term r5rj describes spin dephasing upon reflection, due to reflection 
along axes transverse to Sy || Z. Actually, rı = Re( rari) and rx = Sm(rsri) 
describe reflection along the two axes transverse to Sy. They are called the trans- 


verse reflection coefficients. We will see in a few lines that this complex number is 
involved in the definition of the parameter called the spin mixing conductance 


Guia (1 — rn) [equation (80)| which is encountered whenever there is an 


sd-type transfer of angular momentum. The real Gi = Re( G") and imaginary 
GH = Sm( G") parts of the spin mixing conductance describe the spin current that 


1 
is polarized along the two axes transverse to Sy and that is absorbed by the fer- 
romagnet, therefore creating torques. We recall that no spin-flip process is involved, 
and that the terminology ‘mixing’ accounts for a mixing of the eigenstates. We also 
note that the spin mixing conductance is most often used per unit area A per 


quantum conductance per spin channel Gy/2: g! = aur Typical values are a few 
tenth of nm ?. 
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Complete calculation [M. D. Stiles and J. Miltat, Top. Appl. Phys. 101, 225 
(2006)] shows that, when dephasing and precession are introduced by considering 
complex transmission and reflection coefficients, both DL and FL terms contribute 
to the torque 

AWk 


Ty= Tm. sino [(1 - Re(tt; + nn))s — Sm(tt; + rary) 3] (56) 


_ Ak 
| V2m, 


Ty DET x (5 x Sa) + "Sir x3] (57) 
with gi! 21— (tr ty + rary) the spin mixing coefficient. 

From equation (56), we note that Ty = 0 if t4 = ty and ry = ry. Indeed, in this 
case, there is no spin filtering. We also note that, Ty = 0 if 0 = 0 or 0 = n, when there 
is collinearity between the polarizing layer F yinnea and the free layer Fg, (figure 33) of 
the Fyinnea/N/F free trilayer. In practice, thermal noise always alters collinearity. 

The torque is maximum for 0 = 2/2. However, equation (56) considers the tor- 
que created by one type of electron wave only. The actual torque created by all 
electrons can be obtained by summation over the Fermi surface of the N layer, 
corresponding to all possible incident wave vectors. With this complexity, reflected 
(dephased) spins interact with incoming spins and this interaction affects the initial 
spin polarization. This results in a non-trivial 0-dependence of STT, and especially 
shifts the maximum torque to an angle larger than 0 = 1/2 (figure 38). In addition, 


F/N/F trilayer 
0.5 


T (a.u) 


— 


0.0 


0 90 0 180 


Fic. 38 — Angular(0)-dependence of the torque T. Spin rotation and dephasing upon 
reflection affect spin polarization, contributing to non-trivial 0-dependence of T (full line). 
This trend deviates from the sin@-dependence (dashed line) anticipated from a simplified 
model [equation (56). Adapted with permission from Springer Nature: M. D. Stiles and 
J. Miltat, Top. Appl. Phys. 101, 225 (2006). Copyright 2006. 
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with this complexity, it is possible to explain experimental data showing that, in 
metallic structures, the FL term is small, less than 596 of the DL term. In this case, 


some transverse components average out: Re(tt;) = Sm(tt;) & 0 and 


Sm(ry J « Re ty ty + rp rj) . This contrasts with tunnel junctions, where the 


k-selection upon tunneling reduces the effect of averaging. In this latter case, the FL 
term is about 30% of the DL term. 


Bi Sign convention for the torques 


From equation (57), we get 


v—=T, (58) 


dSy A hk 


"dt =y Im. [Su x (s x Sa) +g! Sm x Q (59) 


The torque Ty is in N.m = J, because it is the torque of a spin angular 
momentum density s acting on a macrospin VSm. The actual torque T on Sy «x M 
requires dividing by V, in which case the torque is indeed expressed in J.m~*. We 
also note that the torque is effective over a finite length governed by sd-exchange and 
spin-orbit coupling (83.5, $4.3, and $4.5). Therefore, the volume V is in fact an 
effective volume and not necessarily the total volume of the ferromagnetic free layer. 
Sometimes, the notation V, is preferred. 

The = of T with respect to magnetization is obtained from equation (59) with 


Sm = WM, 1.€. Su =-M. One gets 


dM A Rk 


Ad lih DE (M x3) + 9M xS] (60) 


Equation (60) thus takes the form 


dM 
——— = — |) T — |} T 61 
di |y] DL |y] FL ( ) 
with 
TDL TDL 
-Tor = |y M x (M xs) = EM x (m x M) (62) 
-Tri = [yl E (M x s) = C (m x M) (63) 
Ms Ms 


where tp, and tp, are the DL and FL B-field amplitudes, in tesla. 
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In the literature, the torque T = — +“ is expressed either in J.m~*, as here and 


[y| dt 
E dM 


throughout this textbook, or in A.m7!.s~!, in which case T = — ^j resulting in a 
factor |y| difference is some of the equations depending on the articles and books 
consulted (see symbols and units page 150). 

Note that, the coefficients tyy) and ry) account for transmission and 
reflection of a spin angular momentum #8 by a ferromagnet of macrospin density 


Sm = -aM . Equivalently, the coefficients tiq) and ry) account for transmission 


and reflection of a spin magnetic moment m = —|y|s by a ferromagnet of 
magnetization M = —|y|$,,. Therefore, tj) = ty) and riy =r: and 


git =1- («t T nnt) eje (trti + rary). 


4.8 STT in CPP Transport Equations 


In the previous section, we have described the transfer of angular momentum from 
the incoming spins to magnetization. We have also detailed how several physical 
effects, which are combined in a single parameter, the spin mixing conductance, 
govern the transfer efficiency. The transfer of angular momentum represents a loss of 
momentum for the spin current, which must be introduced in the diffusion equations 
of transport. More specifically, it is added in the equation describing the flux of spin 
current, in the same manner as the loss of angular momentum for the spin current 
through the spin-flip process was introduced in equation (26). 
Taking into account STT, the flux of spin current reads 


o*(1— f^) o 


VJ, = dE 
Aet) "5^ de sh, 


(M x p,) (64) 


with Uy the average spin-diffusion length defined in §3.1, /,; the ‘exchange’ length 
accounting for the absorption of the transverse angular momentum, and 
o = c(1 — p’). 

Basically, equation (64) accounts for the fact that any loss in spin current must 
be due either to spin-flips, represented by the 1st term, also equation (26) 
(demonstration in exercise 6.3), or to STT, involving damped precession of the spin 
around the local magnetization due to sd exchange interaction, represented by the 
2nd term. In metals, /,; is typically equal to a few nanometers. In the most general 
case, a generic length (see also $3.5 and $4.5) is defined rather than /,4. It also takes 
into account precession around the spin-orbit field. In materials with large spin-orbit 
interactions, like semiconductors, this latter effect is no longer negligible. 

We note that we consider here the three-dimensional case because precessions are 
taken into account. Therefore, in equation (64), the spin current J, is a tensor, and 
its divergence and spin accumulation u, are vectors. This case contrasts with the 
one-dimensional case of equation (26), in which J, was a vector, and its divergence 
and the spin accumulation u, were scalars. 
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Typical critical current densities for switching magnetization by STT are 
J «ait ~ 10" A.cm"?, hence the need for nanometer-size lateral dimensions (figure 39, 
top-left) — in addition to nm-thick layers, because the layers’ thicknesses need to be 
smaller than D. 


From figure 39, we observe that switching is polarity-dependent and that the 


positive critical current density is larger than the negative one: J RES pa. ui M 


fact, for a positive electron flux J spins travel only once across the N 


ecrit? 
(Cu) layer before acting on the Free layer, whereas for a negative electron flux 


Jia Spins travel twice across the N (Cu) layer (figure 39, bottom). As a result, 
+ 


e,cri 


e,crit 


the exponential increase of J/ is twice faster than J (figure 39, top-right). 
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Fic. 39 — (Top-Left) Experimental result showing the current induced STT-switching of the 
magnetization of the free-ferromagnet, and the corresponding change of resistance due to the 


GMR effect; inset: schematics of a representative pillar used for STT experiments, making use 
of a pinned-ferromagnet /non-magnet/free-ferromagnet trilayer, reprinted with permission 
from the American Physical Society: J. A. Katine et al., Phys. Rev. Lett. 84, 3149 (2000). 
Copyright 2000. (Top-Right) Dependence of the critical switching current on the thickness of 
the non-magnet, for two cases: Jý when electrons flow from the free- to the 
pinned-ferromagnet, and J, in the opposite direction. Reprinted with permission from the 
American Physical Society: F. J. Albert et al., Phys. Rev. Lett. 89, 226802 (2002). Copyright 
2002. (Bottom) Illustration of the STT process in the two cases presented, with 
dM / dt = —|)|T. 
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4.4 STT in Magnetization Dynamics Equations, 
Reciprocal Spin Pumping 


E The STT terms in magnetization dynamics equations 


Naturally, the STT terms can be introduced in the Landau-Lifshitz-Gilbert 
(LLG) equation, which describes magnetization dynamics [basics about the LLG 
equation can, for example, be gained in J. M. D. Coey, Magnetism and Magnetic 
Materials, Cambridge University Press (2010)]. The result of adding the STT terms 
leads to 


dM 

E = —|y|T Larmor = IE Damping EE \y|T pt a \y|T er 

dt 

dM a dM TDL 

— = -—|)|M A. M M M 

de M x (HoHete) + We ae ll im e (m x M) (65) 


+ lil ype (m x M) 


The first term in equation (65) represents the Larmor precession around the 
effective field Hepp. It includes the contributions of the external magnetic field and 
internal fields due to anisotropy, exchange and dipolar interactions. We recall that 
the gyromagnetic ratio y = — 5£- g, with g the Landé g-factor, is negative because of 


2m, 
the negative charge —e of the electron, making the precession anticlockwise viewed 
from above. In the literature, the gyromagnetic ratio y = — 57 g is used either as is 


in the equations, 2M ^ic = YM x (ugHar)..., or in absolute value |y|, as throughout this 
textbook, M- = —|y|M x (uHeft)..., or a by Yọ = —yuo; to make it positive 
and include fly to match the cgs Dd system, 4M — = —yM x Her.. ., resulting in a 
factor —1 or cji difference is some of the equations depending on the articles and 
books consulted (see symbols and units page 150). 

The second term in equation (65) represents damping «, due to internal relax- 
ation ascribed to intraband and interband scattering (exercice 6.5), and to external 
relaxation, for example, due to spin pumping (exercise 6.7). 

Finally, the last two terms represent the DL and FL components of STT, 
respectively. How the two contributions to STT were named is clear from the LLG 
equations and corresponding representation (figure 40), because the DL torque is 
along the damping or anti-damping torque, and the FL torque is along the Larmor 
or anti-Larmor torque. 

Because the FL term is most often negligible in metals ($4.2), we focus on the DL 
term to estimate the critical current needed to trigger magnetization dynamics 
(precession or switching, figure 40). The amplitude of the DL torque is expressed as a 


function of the material's parameters as follows: tpr, — AE. T 1 , with dp the thickness 


of the magnetic layer, f the polarization, and Je the amplitude of the charge current 
density. 

The sign of the torque depends on the current direction. We recall that, in 
regards to the orientation of the torques, the rule is as follows: when an electron of 
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Heff Damping 


Stable precession Switching 


dM Torque = damping Torque » damping 


dt 


Anti-damping 
Torque 


Fic. 40 — (Left) Schematics of magnetization dynamics when the magnetization M is sub- 
jected to several torques: the Larmor torque created by the (effective) magnetic field (Her), 


the Damping torque due to energy transfer, for example, to electrons, magnons and phonons, 
and the Anti-Damping component of the spin transfer torque (STT). (Right) Calculations 
illustrating two cases: stable precession and switching, depending on the relative magnitudes 
of the torques involved. Reprinted with permission from Elsevier: D. C. Ralph and M. 
D. Stiles, J. Mag. Magn. Mat. 320, 1190 (2008). Copyright 2008. 


magnetic moment density m (or equivalently of spin density s pointing in the 
opposite direction) flows toward a magnetization M, the DL torque orients M 
toward m (away from s). 

The critical current density is obtained when the (anti)DL torque compensates 
the damping torque in each period of oscillation of M 


a 
|y|TDL crit = 2n = 2naf (66) 
"S ejto Ms drz 
J ecrit = ^f (67) 


'The resonance frequency of a ferromagnet is in the GHz range. For an isotropic 
thin film with a planar magnetization subject to an applied field in the plane, 
Kittel’s law describing the frequency-field relation is: 2zf = |y|ug/ H(H + Mg) 
[J. M. D. Coey, Magnetism and Magnetic Materials, Cambridge University Press 
(2010)]. For the calculation, we took the frequency corresponding to the zero-field 
value of the high-field asymptotic dependence: 2zf = hine Ms. 

The reasoning is better when considering energies. The critical current density is 
obtained when the energy lost due to damping is compensated by the energy gained 
due to antidamping by ST'T in each period of oscillation of M. 

The parameters for NiggFes) permalloy are « = 0.008, Ms = 0.8 MA.m"!, and 
p — 0.3. If we consider dr — 3nm, and the values of the physical constants given 
page 150, we obtain J ecrit ~ 107 A.cm ?, indeed one order of magnitude larger 
compared to the current needed to observe the GMR effect (83.3). We also get 
tp [Lu T] ~ 40J [107 A.cm~?]. 

ST'T is now used in commercial MRAM applications, for writing, when TMR is 
used for reading. How the fundamental effects described in this book are used in 
MRAM applications, is detailed in B. Dieny, R. B. Goldfarb, K.-J. Lee (eds) Intro- 
duction to Magnetic Random-Access Memory, Wiley (2016). 


Transfer of Angular Momentum — STT, Spin Pumping, SOT 63 


E STT - spin pumping reciprocity 


Spin pumping and spin transfer torque are reciprocal effects [Y. Tserkovnyak et al., 
Rev. Mod. Phys. 77, 1375 (2005)]. An intuitive picture consists in comparing spin 
transfer torque to a water flow (spin current) moving the blades of a watermill 
(magnetization) and spin pumping to moving blades (magnetization) creating a 
water flow (spin current) (figure 41). The spin pumping effect will be described in 
detail in exercises 6.6 and 6.7. 


Spin transfer torque Spin pumping 
Js puts M into motion The motion of M induces J, 


Fic. 41 — Schematics showing (Left) spin transfer torque and (Right) the spin pumping 
reciprocal effect. The drawings are courtesy of D. Gusakova. 


In short for now, the spin pumping effect refers to the ability of a magnetic 
material, when brought out-of-equilibrium, to generate a spin current (J?) at an 
interface, due to transfer of angular momentum. The technique usually involves 
inducing resonance in a ferromagnetic (F) spin injector — e.g., a NiFe layer — which is 
adjacent to a non-magnetic material (N) known as the spin sink — e.g., a Pt layer 
(figure 42). Spin pumping is therefore treated as a loss of momentum for the spin 
current, due to transfer of angular momentum at the F/N interface, thus impacting 
spin accumulation through the addition of a relaxation channel [equation (26)] and 
magnetization dynamics through the addition of a STT term [equation (65)]. This is 
similar to what we have described in the previous sections, but here, relaxation and 
STT are induced by a pure spin current rather than a spin-polarized charge current. 

Transport: understanding which material's parameters drive spin pumping (and 
a subsequent spin-charge conversion) is the object of exercise 6.6. In this exercise, 
based on the diffusion equation for spin accumulation [equation (27) we will 
understand why the spin current propagating in the N layer, along y, reads 


AY — 4s sinh (dy/l; y) 


and why this gives rise to a transverse voltage, due to the inverse spin Hall effect 
(ISHE) 


(68) 
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NiFe 


Fic. 42 — (a) Illustration of the spin pumping effect due to sustained out-of-equilibrium 
magnetization dynamics. (b) Illustration of spin current diffusion across the N layer (here, 
Pt). (c) Equivalent circuit considering spin current to charge current conversion due to the 
inverse spin Hall effect in the N layer. Adapted with permission from AIP Publishing: K. Ando 
et al., J. Appl. Phys. 109, 103913 (2011). Copyright 2011. 


wOssl tnb (dy /(20%)) ais ol! 
dyon + dFOFp 8zo2c 


Visug = sin(0y)T (69) 
with w the width of the layers, Üsug the spin-Hall angle, h,; the rf excitation field 
used to reach F resonance, and o the resonance angular frequency. T is a parameter 
accounting for the trajectory of the magnetization precession in the zy plane. It can 
be viewed as the elliptic- to circular-trajectory ratio. g!! is the real part of the spin 
mixing conductance per unit area per quantum conductance per spin channel, 
accounting for the ability of the F/N interface and N layer to absorb the spin 
component along M x (m x M). 


Magnetization dynamics: how spin pumping gives rise to extra damping for mag- 
netization dynamics is the object of exercise 6.7. We will show that the additional 
torque term due to spin pumping reads [R. Igushi et al., J. Phys. Soc. Jap. 86, 
011003 (2017)] 


h a dM 
-TSE = , : = p 
v (oh) = pg a me) 
where 
= I| h TL (71) 


is an extrinsic damping term showing that spin pumping creates a loss of angular 
momentum for M, with a the real part of the effective spin mixing conductance 
per unit area per quantum conductance per spin channel, that takes into account 
back reflections of the spin current at the outer interfaces of the F/N bilayer, 
compared to g!! which does not. 
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E Note about the general expression of the torque due to a spin or spin-polarized 
current 


Another way of finding the expression of the torque created by a spin current or 
spin-polarized current on magnetization is shown in figure 43. It consists in looking 
at the conservation of total magnetic moment (without spin-flip relaxation) by 
analogy with the conservation of charge [equations (72) and (73), where p is the 
charge density defined in equation (47)]. From this, we get equations (74) and (75). 
Consequently, the torque created by a flow of spin 8 angular momentum density 
Q= ÍJ on a magnetization M is - T = V- (£Js), with J, = -( -j'), where the T 
and | symbols refer to the orientation of the moment m. 

We note that, naturally, the torque created by a flow of spin 8 angular 
momentum density Q — iJ on a macrospin angular momentum density Sm is of 
opposite sign, as Sọ and M are of opposite sign, in agreement with equations (52) 
and (58) from the toy quantum mechanical model. 


nr - ff «s (72) 


-VJe (73) 


Il com f t ones s 


dM h 
Ea bv ($a) = -ir (75) 
HBJ 4 
` Je aa - 
Ni x Oe- A, N 
ap dM HBJ 
-_— — HBI . 


Fic. 43 — Illustration of (Left) charge conservation, and (Right) magnetic moment 
conservation. 
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4.5 STT in Magnetic Textures 


Up until now, we considered STT in the case of uniform magnetization. In this 
section, we will deal with the case of non-uniform magnetization in magnetic tex- 
tures like domain walls (DWs). This part echoes $3.5 in which we discussed DW 
magnetoresistance. 

Taking aM = —|y|T, the STT in the case of a DW is expressed as follows 


np ABC 


T= +— Mx (Je VM (76 
2eM$ (1+ &) 2eM$(1-- &) SIM Se 


M x [M x (J.- VM] 


We recall that here, for consistency, f is the polarization of the current, not to be 
confused with the ‘beta term’ found in the literature, which takes into account the 
adiabaticity of the process. 

In fact, two contributions are usually considered, depending on the ability of the 
spins 8 (or equivalently moments m = —|y|5) of the itinerant s-electrons to follow 
magnetization M (figure 32): 


Bi the adiabatic (DL) contribution, when there is perfect tracking, represented by 
the 1st term in equation (76), and 

Bi the non-adiabatic (FL) contribution, when there is mistracking; 2nd term in 
equation (76). 


'The adiabatic contribution is very similar to the mechanism described in $4.2, in 
which the spin current or flux of spin density transfers angular momentum to 
magnetization due to sd-interactions, therefore generating a predominant DL tor- 
que. In the non-adiabatic contribution, spin accumulation takes place due to 
spin-mistracking. It then precesses around the magnetization and generates a FL 
torque upon relaxation. 

We recall (83.5) that the process of spin-mistracking is driven by the 
spin-mistracking parameter č. In low spin-orbit materials, € = 23 (4). with 
lea = V Di,, and lI, = \/Dt.. The precession of itinerant electrons around the 
spin-orbit field enhances €. In weakly spin-orbit-coupled dense-moment F metals like 
Co, č « 1, and the DL term dominates the FL one: DL > FL. In contrast, in 
strongly spin-orbit-coupled dilute-moment F semiconductors, like (Ga,Mn)As, 
č X» 1, and the FL term dominates: FL > DL. 

In the same way that STT can trigger magnetization dynamics in a uniform 
system, it can also move magnetic textures. The torques contributions distort the 
initial Bloch configuration (figure 44), therefore redistributing the initial conditions 
for the calculation of T. This point underlines that there is constant feedback between 
VM and T. 

A possible application of STT in magnetic textures is the racetrack memory 
(figure 45) [S. S. P. Parkin et al., Science 320, 190 (2008)], in which it is possible to 
store information in a 3-dimensional way, unlike current storage technologies. The 
interested reader of 3-dimensional spintronics applications may consult focused 
reviews like A. Fernández-Pacheco et al., Nat. Commun. 8, 15756 (2017). 
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Fic. 44 — Schematics illustrating how current-induced (Left) damping-like (DL) and (Right) 
field-like (FL) contributions to dM/dt distort and move a Bloch magnetic domain wall. 
Adapted from T. Janda et al., Nat. Commun. 8, 15226 (2016) — CC BY license. 


DW injectors 


Sensors 


3D racetrack memory 


Fic. 45 — Schematics illustrating the principle of racetrack memory. Reprinted with per- 
mission from Springer Nature: S. S. P. Parkin et al., Nat. Nanotech. 10, 195 (2015). Copyright 
2015. 


4.6 Magnetoelectronic Circuit Theory 


In this section, we want to briefly mention the magnetoelectronic circuit theory of 
transport [A. Braatas, G. E. W. Bauer, and P. J. Kelly, Phys. Rep. 427, 157 (2006)]. It 
is a generalization of circuit theory to non-collinear magnetization. The findings of 
this theory are now commonly used in the field of spintronics, be it related to 
CPP-GMR, STT, spin pumping, or other effects like the spin Hall magnetoresistance 
(SMR), which is treated in exercise 6.8. Because CPP-GMR, STT, and spin pumping 
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are companion effects, it is possible to measure the same type of parameters, like the 
spin mixing conductance G!! and the spin diffusion length D in different ways. 
While this section is intended for advanced readers, the comparison of the findings of 
the model with the toy model for STT discussed in $4.2 is within reach of all readers. 

In this theory, the spatial dependence of spin accumulation is disregarded. 
Spin-accumulation is reduced to a vector-potential V; y and the N layer is reduced to 
a node with a set number of net spins that have accumulated Ne|V, y| (figure 46). 
The initial model is therefore appropriate when interfacial resistances prevail. We 
note that, in this book, for the sake of consistency, V; y refers to the vector potential 
representing electrons with spin s in the N layer, meaning that in the F layer, 
V.p—8V.p—-—mV,. 


Spin accumulation 
Charge vector i 
-potential 
potential ~> P 
Ven, Von 


reservoir reservoir 


V, 


conductor node conductor 


Fic. 46 — Schematics of the equivalent circuit used to model STT in the frame of the 
magnetoelectronic circuit theory of transport. Adapted with permission from Elsevier: 
A. Braatas, G. E. W. Bauer, and P. J. Kelly, Phys. Rep. 427, 157 (2006). Copyright 2006. 


The transverse spin 8 current and current density, at the interface, arriving from 
the layer N to the layer F; (figure 46) have the form 


Gu el 
J,A=2 em x (Mi x Vs) + 15 X Ven (7T) 
E 
2G! By 2G) a, 
aay VT a UE *) | zu 75 E (78) 


where Gi = Re( G) and G! = Sm(G") are the real and imaginary parts of the 
spin mixing conductance, accounting for the spin current that is polarized along the 
two axes transverse to M, and that is therefore absorbed by the magnetic material. 
We recall that the absorption of the transverse angular momentum creates a torque. 
Given the transverse direction to which they refer, an results in the damping-like 


torque and GI! results in the field-like torque. In the literature, the spin mixing 
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conductance is defined either as G'!, in siemens as here and throughout this 
textbook, or as G!! / A, in which case it is in units of S.m7?, resulting in a factor A 
difference in some of the equations depending on the articles and books consulted. 
We also recall that the spin mixing conductance is most often used per unit area A 


per quantum conductance per spin channel Go/2: g!! = deme It is customary to 


incorrectly use the term spin mixing conductance for g!!, although its unit is m^? 


(see symbols and units page 150). 
The torque acting on F,, due to absorption of the transverse spin current is then 


h h | at ql 
-T = s. = r 1 
V (5.4) OV ma x (M x py) + We X Ms (79) 


where 


neN 


Gt = >J | -Nnrí(n) ] (80) 


The integers n, m label what are called the propagating modes or scattering 
channels. They actually relate to the wave numbers kpm at the Fermi level. The 
elements THY) are the reflection coefficients for a transverse mode m with spin ft (4) 


at the N side reflected to a transverse mode n with spin ft (4)). 
We recall that the coefficients ty) and rpq) account for transmission and 


reflection of a spin 8 by a ferromagnet of macrospin density Sy = — uM. Equiva- 

lently, the coefficients ty) and ryq) account for transmission and reflection of a 

y 
2 : * 

hy = tay and ry =r, and G'-TYE4b-Yuan(n)|- 

2 * 

ED nen [1 m went nm | : 


Combining equations (79) and (80), it is possible to express the spin transfer 
torque in the magnetoelectronic circuit theory 


T= -— (s e (Y vn) x (m x m) 


meN 


— w(X (n) ju x n 


meN 


moment m = — 


4 by a ferromagnet of magnetization M = —|)|S,,. Therefore, 


(81) 


where N represents the number of scattering channels in the N layer. 

This expression is very similar to that obtained with the toy quantum mechanical 
model [equation (57)], with its one type of electron wave — which can be refined with a 
summation over the Fermi surface. For the comparison, we recall that the notations 
M and Sy refer to the magnetization and macrospin density of the ferromagnet 
receiving the flux of angular momentum, and gt, s, and m to the spin s accumulation, 
spin 8 current density, and magnetic moment m density of the incoming electrons 
bringing a flux of angular momentum (see symbols and units page 150). 
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4.7 Spin-Orbit Torques - SOT 


In this section, we now turn to another type of torque that exists in spintronics: 
spin-orbit torque (SOT). We describe several types of SOT and show how they result 
from the transfer of angular momentum from the electron current to magnetization, 
through interactions between lattice and electron orbital (crystal field potential), 
electron orbital and spin (spin-orbit coupling), and spin and magnetization 
(sd-exchange interactions). 


$46 — X 


coupling 
orbital spin 
crystal field exchange 
potential interaction 
% d. Q 
9e. 9-0 
Q W Q 
9 [o 
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lattice local magnetic moment 


Fic. 47 — Schemes of four reservoirs of angular momentum in crystals — spin and orbital of the 
electron, the crystal lattice, and the local magnetic moment — and the sources of interactions 
and transfer/interconversion between them — (sd)-exchange interactions (see also chapter 4), 
spin-orbit coupling (82.6), and crystal field potential due to Coulomb interactions (chapter 5). 
From D. Go et al., Phys. Rev. Res. 2, 033401 (2020) - CC BY license. 


The concept of SOT refers to a torque created by an electric current on a 
magnetization, the origin of which involves a spin-orbit coupling mechanism. More 
generally, this type of torque involves interaction processes between different sources 
of angular momentum in crystals (figure 47): 


W sd-exchange interactions between magnetic moments due to localized electrons 
and the spin of the itinerant conduction electrons (see also chapter 4), 


WB spin-orbit coupling, between the spin and orbital angular momenta of the itin- 
erant electrons (see also $2.6), 


B Coulomb interactions (crystal field potential), coupling the lattice and the orbital 
angular momentum of the itinerant electrons. This interaction has given rise to a 
recent research subfield called orbital spintronics [D. Go et al., Europhys. Lett. 135, 
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37001 (2021)], in which no spin-orbit coupling is initially required to trigger the 
spintronic effect (see also chapter 5). In the case of torque, spin-orbit coupling is 
needed in a second step, to couple the orbital angular momentum and the mag- 
netization to be set in motion. This type of torque is also classified as an orbital 
torque. 


SOT is most often found in non-magnet(N)/ferromagnet(F) bilayers for a cur- 
rent in-plane (CIP) configuration. In contrast to the STT structures described in 
§4.1-4.6, SOT structures do not require an additional F-layer to create a 
spin-polarized current (i.e., a flux of angular momentum). In SOT structures, not 
only spin currents (flux of spin angular momentum, J,), but also non-equilibrium 
moment densities m (equivalently spin densities s = —m/|y|) are generated in the N 
or F layer, due to spin-orbit coupling. 

Spin-orbit effects can be readily understood by considering the motion of an 
electron in a potential gradient V $ and how fields are transformed between inertial 
frames. The net electric field, created by the potential gradient becomes a magnetic 
field in the electron's rest frame H so(k). This momentum(K)-dependent magnetic 
field couples to the magnetic moment of the electron 4n, through a Zeeman term, 
-Mm : uH solk). In magnetic materials, this interaction interconnects the direction 
of electron flow and the magnetic order parameter. It is said that the spin angular 
momentum is locked to its wavevector. Atomic spin-orbit interactions, where œ is 
the potential with lattice periodicity, can be distinguished from effective spin-orbit 
interactions, which result from combining the atomic spin-orbit interaction with 
particular crystal or multilayer systems. In crystals or structures lacking inversion 
symmetry, the k dependent magnetic field becomes odd in momentum, resulting, for 
example, in Dresselhaus or Rashba fields (85.4). [Text adapted from V. Baltz et al., 
Rev. Mod. Phys. 90, 015005 (2018)]. 

Several mechanisms have been identified as being at the origin of spin-orbit 
torque, including the bulk spin Hall effect (SHE) and bulk or interfacial inverse spin 
galvanic effects (iSGE), like the Rashba-Edelstein effect (REE) (figure 48). In 
prototypical models, beyond the origin of the effects, the main difference between 
SHE and iSGE is that the former generates a flux of angular momentum between the 
N and F layers, whereas the latter generates a non-equilibrium spin density at the 
N/F interface. As a result, SHE produces a damping-like torque, through absorption 
of the flux of angular momentum by the F layer, similar to STT (84.1—4.6), and 
iSGE produces a field-like torque due to a direct sd-exchange mediated by con- 
duction electrons in the F layer, which impinges on and is reflected at the interface. 
In refined models, and in practice, each torque has both a damping-like and a 
field-like component, whether induced by SHE or iSGE. They are of the form 


-T = TM x (m x M) + Tem x M (82) 


where m || J. x u is the non-equilibrium moment density, with 4 a unit vector 
determined by the symmetry of the system. In the case of a magnetic multilayer 
perpendicular to Z, and experiencing SHE or REE as depicted, i = Z. Two types of 
these SOT are detailed below. 


72 The Basics of Electron Transport in Spintronics 
i —lylT pL 
M mLa M 
- ues t TTS m 
=~ ise 


—lylTrL 


Bulk spin Hall effect Interfacial Rashba-Edelstein effect 
(SHE) — induced torque (REE) — induced torque 


Fic. 48 — Schematics of various origins of the CIP-induced spin-orbit torque in a non-magnet 
(N)/ferromagnet(F) bilayer. (Left) Spin Hall effect (SHE) takes place in the bulk of the N 
(heavy metal) layer and creates a flux of spin s angular momentum J; across the interface. 
(Right) Inverse spin galvanic Rashba-Edelstein effect (REE) arises at the interface with the N 
(heavy metal) and creates a non-equilibrium moment density m (equivalently a 
non-equilibrium spin density s = —m/||) at the interface. 


E Bulk SHE-induced torque (figure 48, left) 


Due to bulk spin-orbit coupling in a metallic N layer, electrons of opposite spin 
scatter in opposite directions [N. Nagaosa et al., Rev. Mod. Phys. 82, 1539 (2010)]. 
This effect is called the spin Hall effect (SHE), which we admit for now. Thus, a 
longitudinal charge current is transformed into a transverse spin current in the N 
and diffuses in the F. This flux of angular momentum creates a torque on the 
magnetization of the F. The torque is thus similar to that described in $4.1—4.6, 
although the origin of the flux of angular momentum differs. T'herefore, in the case of 
a pure SHE-induced spin torque, the torque itself is of the STT type. Exercises 6.6 
and 6.8 on ISHE and SMR, respectively, give examples of similar physics and how 
similar torques can be calculated (see also figure 56 and corresponding text). 


Damping-like component of the bulk SHE-induced torque. 
For G!! = 0 [equations (60)-(63)], the SHE-induced torque takes the form: 


SHE 
ay) dM TDL 
[SHE = MOD M = M M 83 
STT \y|M 5 x dt M? x (m x ) ( ) 
From this equation, we see that 1 is damping-like (same as in §4.1-4.6) and 
SHE 


thus creates an additional damping term ocyrop. This type of torque enables mod- 
ulation of the F’s damping [K. Ando et al., Phys. Rev. Lett. 101, 036601 (2008)]. 

Linearizing the LLG equation for magnetization dynamics [equation (65)], it is 
possible to show that 
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1 h 
SHE si 0 
a sin(0 J 84 
oe (Hoa + “ tet) toM g dp ( oe j (84) 


The prefactor before sin(0) i J? accounts for the fact that the trajectory of the 
magnetization is hindered by effective fields, including anisotropy and applied field. 
The sin(0)-dependence, where 0 is the angle between Je and M, results from the 
vector product at the origin of the torque. 

Since the spin accumulation and the corresponding transverse spin current Y 
arise from the longitudinal charge current Je, through the SHE process, we have 
J? = Üsux1J.. The spin-orbit mechanism is included in the spin Hall angle sug (see 
figure 56 and corresponding text, also exercise 6.6). 

The parameter 7 is the result of spin accumulation. It takes into account: 
(i) spin flip process in the non-magnetic layer, i.e., that only a fraction of the 
converted spin current matters as the conversion far from the interface will be 
depolarized by the time it reaches the interface, and (ii) the efficiency of the 
transfer of angular momentum, through the spin mixing conductance (per unit 
area per spin channel) g!! = G!!/(AG)/2). Based on the drift-diffusion equation 
for spin accumulation [equation (16) for z = 0 in Y. T. Chen et al., Phys. Rev. B 
87, 144411 (2013), see also exercise 6.8], and considering that the spin accumu- 
lation in the F layer is negligible, and that gi — 0 (84.2) it can be shown that 
n= g| tab (5 J + gj! coth (79 )]. 


Uy N 


The additional damping due to SHE-induced STT is thus 


dy 
SHE sin(0) h a J {tanh (3) 
io = SHE” e 
MOD - (Ha m Ma) e 


(85) 


Experimentally, it is most often possible to measure gl he rather than git, for 
example, from the additional damping due to spin pumping [equation (71)]. How those 
two parameters are related is the object of exercise 6.7. We recall that His takes into 
account back reflections of the spin current at the outer interfaces of the F /N bilayer. 


Field-like component of the bulk SHE-induced torque. 


A field-like contribution to the SHE-induced torque appears when ai # 0, i.e., when 
considering the spin current that is perpendicular to both the magnetization and the 
spin accumulation, similar to what was discussed in the section on STT $4.1—4.6. 


Bi Interface REE-induced torque (figure 48, right) 


Field-like component of the interface REE-induced torque. 

We remind that in systems lacking bulk or interfacial inversion symmetry, the 
spin-orbit coupling becomes odd in momentum. Examples are linear and cubic 
Dresselhaus spin-orbit coupling in strained Zinc-Blende semiconductors, Rashba 
spin-orbit coupling in bulk Wurtzite and at interfaces between dissimilar materials 
as well as cubic Rashba spin-orbit coupling in oxide heterostructures. The main 
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result of such a type of spin-orbit coupling is that the spin angular momentum is 
locked perpendicular to the momentum k, as shown in figure 49, left and middle for 
the case of the interfacial Rashba coupling. 


A naive yet pedagogical derivation of the Rashba Hamiltonian Hp considers that, 


due to the Lorentz transformation uH so = — (XE) electrons experience a magnetic 


field H so in their reference frame when moving in the electric field E = E,Z created by 
P-breaking at the surface (figure 48, left). The corresponding Hamiltonian is 


Hr = —m H so = #3 (v x E) - o = HP (v x e,) -o = arlk x 0) Z, where xg = 


dà 
gig Eoft 
2m,c? 


is the Rashba parameter. In a refined model, the Dirac gap at play in the naive 
model m,c?-1MeV is replaced by an energy term of ~1eV, which results from 
splittings in the energy bands, and ag ^ 0.1 — 1 eV.À. 

At equilibrium, when J, = 0, the moment density is zero (figure 49, left and 
middle). In contrast, a non-equilibrium (Edelstein) moment density m x Je x U x 
ary builds up on the application of a current J, = J,x (figure 49, right). This 
moment density is proportional to the Rashba spin-orbit coupling parameter ag. It 


leads to a field-like torque of the form —TS57 = fm x M, due to direct sd-exchange 


mediated by conduction electrons in the F layer [A. Manchon and S. Zhang, Phys. 
Rev. B 78, 212405 (2008)]. 


J sax, or 


Fic. 49 — Illustration of the Rashba-Edelstein effect (REE). (Left and Middle) 3D-view and 
cross sections at the Fermi level of the band structures for spin(|J) and spin({}) electrons at 
equilibrium, for Je = 0, corresponding to spin-orbit induced parity symmetry breaking of the 
Rashba type, resulting in angular momentum locking. The Hamiltonian reads H = 
hk’ /(2m,) + ap(k x o) -2 (85.4). (Right) Edelstein effect in which a non-equilibrium excess 
spin density s (equivalently moment density m = —|y|s) along —y (y) builds up for 
J. — J,X z 0. 


Damping-like component of the interface REE-induced torque. 

'The ideal picture described above considers a two-dimensional spin density, localized 
at the interface. In practice, spins may accumulate in three dimensions instead. As a 
result, the spin accumulation can diffuse from the N to the F and thus create a flow 
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of angular momentum, whose absorption by the F layer generates a DL torque of the 
spin transfer type (§4.1-4.6). Another source of DL torque may come intrinsically 
from the lattice and can be formulated in the language of the Berry curvature 
(chapter 5). It is described intuitively by considering the spin-dependent (due to 
angular momentum locking) electron spin dynamics governed by H so(k) during the 
acceleration phase, between the scattering events [H. Kurebayashi et al., Nature 
Nanotech. 9, 211 (2014)]. 

It is now recognized that SOTs can be classified according to whether the torque 
originates from the CIP electric current flowing in the N or F, or whether it is due to 
spin-orbit coupling in the N or F layer [V. P. Amin and M. D. Stiles, Phys. Rev. B94, 
104420 (2016)]. We recall that the SHE-induced torque originates from the electric 
current flowing in the N, whereas the REE-induced torque originates from the current 
flowing in the F. In both cases, the torque is due to spin-orbit coupling in the N, in the 
bulk for the SHE, and at the interface for the REE. The contribution of spin-orbit 
coupling in the F-layer can also create torques. Because of their common origin, by 
analogy with the anomalous Hall effect, they are called AHE-induced torques. We will 
also see in chapter 5 that intrinsic effects may give rise to the so-called spontaneous 
anomalous Hall effect (SAHE) which can also give rise to SAHE-induced torques. In 
this last case, the effects arise from the spin-lattice asymmetry, sometimes indepen- 
dent of spin-orbit interactions, and the spin-orbit interaction is only involved in 
coupling the spin and the crystal lattice, as a perturbation. Other torques, including 
orbital Hall-effect-induced torques, may also contribute to some structures. Interested 
readers may consult for example the focused articles by D. Go et al., Phys. Rev. Res. 2, 
033401 (2020), and A. Manchon et al., Rev. Mod. Phys. 91, 035004 (2019). 

We note that the torques discussed above can contribute concomitantly. 
Examples for the case of Fe/W and Ni/W bilayers are in D. Go et al., Phys. Rev. Res. 
2, 033401 (2020). Using insulating or conducting stacks, or varying the layers' 
thicknesses may allow a particular type of torque to be selected. We also note that 
some of these torques can also be encountered when the F is a metal, an insulator, a 
topological insulator, a 2D material, and if it is replaced by some ferrimagnets or 
antiferromagnets. 

A standard technique to measure spin-orbit torques is based on the harmonic 
analysis of the periodic Hall voltage following the application of a periodic charge 
current. This technique is introduced in exercise 6.9. 


E SOT in magnetic textures: the example of the skyrmion Hall effect 


Similarly to STT (84.1—4.6), and thus already discussed in detail, SOT terms can be 
added to the transport and magnetization dynamics equations and influence the 
motion of magnetization, be it uniform or textured [D. Go et al., Phys. Rev. Res. 2, 
033401 (2020); A. Manchon et al., Rev. Mod. Phys. 91, 035004 (2019)]. To mention 
just one, a striking example is the SHE-induced torque motion of a skyrmion in an F 
layer adjacent to a N layer, in which the effect known as the skyrmion Hall effect 
deviates the skyrmion laterally from the electron flow. 

Magnetic skyrmions are nanoscale chiral spin textures that are categorized 


based on their topological index Q=7/f M. (om x 3M ) dady 
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Skyrmion trajectory 


Fic. 50 - Illustration of the SHE-induced torque motion of a skyrmion, in which the skyrmion 
Hall effect deviates the skyrmion laterally from the electron flow. Reprinted with permission 
from Springer Nature: G. Chen, Nat. Phys. 18, 112 (2017). Copyright 2017. 


[A. Soumyanarayanan et al., Nature 539, 509 (2016)]. As an example, a Néel-type 
skyrmion is shown in figure 50. By considering a steady and rigid behavior, the 
in-plane (x, y) motion of the center of mass of a skyrmion with velocity v is best 
described by a modified Thiele equation (equation 86). It is derived from the 
known LLG equation (§4.4) in which transport-specific effects are accounted for, 
such as the damping contribution of the electron spin due to coupling to the 
skyrmion by the Berry phase and the forces on the electrons arising from the Berry 
phase they pick up when their spin adiabatically follows the topology of the 
moments in the skyrmion (see chapter 5 for a description of the Berry phase and 
the topological Hall effect, companion to the skyrmion Hall effect). Detailed cal- 
culations showing how the transport-specific effects for a skyrmion are included in 
the LLG equation can, for example, be found in K. Everschor et al., Phys. Rev. B 
86, 054432 (2012). The Thiele equation thus writes: 


Gxv—aD-v+F=0 (86) 


The first term in the equation describes the topological Magnus force. For a 
Néel-type two-dimensional skyrmion, one gets G = — Mete Am Qz, known as the 
gyromagnetic coupling vector, indeed related to the skyrmion topological index Q. 
'The second term represents the dissipative force. It is governed by the dyadic tensor 


D, with Dy — Mate [ (ðM - 0M) didj. Here, D= D, = Dy, and Dry = Dy, = 0. 
'The last term is the force at the origin of the motion. For the case of a SHE-induced 
torque, F relates to the force on the electrons arising from Berry phase inherent to 


the magnetic topology of the skyrmion. It writes F = 41B - Je, where B is a tensor 
related to the emergent topological magnetic flux density, and Je is the current 
applied in the adjacent N layer. We recall that the corresponding transverse spin 
current at the interface creating the torque is J : = Ospen J e- 
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For a current applied along X, it can be shown from equation (86) that the 
deviation of the trajectory from the electron flow is: 2 == 2 [W. Jiang et al., Nat. 
Phys. 13, 162 (2017)]. 


Summary 


In this chapter, we have seen that sd-coupling between an incoming moment 
density m (equivalently spin density s = —m/|y|) and a magnetization M results 
in a transfer of angular momentum. This angular momentum transfer alters spin 
accumulation and gives rise to a transverse spin transfer torque STT, which is 
decomposed in a damping-like component along M x (m x M) and a field-like 
component along m x M. The process of angular momentum transfer is governed 
by a complex number called spin mixing conductance G'!, which relates to spin 
rotation and dephasing upon reflection and transmission, leading to a mixing of 
the eigenstates. It is most often used per unit area A per quantum conductance 
per spin channel Gy/2: g! = e and its real and imaginary parts relate to the 
damping-like and field-like components, respectively. The transfer of angular 
momentum alters spin accumulation, through the addition of a relaxation 
channel, and magnetization dynamics, through the addition of a STT term. With 
a sufficiently large current, determined by the material's properties and system 
geometry, STT can trigger magnetization dynamics, precession or switching, of a 
uniform magnetization, and propagation of magnetic textures. Importantly, 
many effects in spintronics, like CPP-GMR, STT, spin pumping, and SMR, 
involve a transfer of angular momentum of the same type. They can thus be 
combined to accurately determine key spintronics parameters, like the spin 
mixing conductance g!! and the spin diffusion length lj. Other types of torques, 
called spin-orbit torques SOT, can influence transport and magnetization. These 
torques originate from spin-orbit coupling mechanisms. They are most often 
found in N/F bilayers for a current in-plane CIP configuration. In SOT mecha- 
nisms, not only a flux of angular momentum like in the STT mechanism but also 
non-equilibrium spin densities can be generated. 


Chapter 5 


Berry Curvature, Parity and Time 
Symmetries — Intrinsic AHE, SHE, QHE 


In this chapter, we discuss an alternative approach to the previous chapters, which 
relied on extrinsic scattering. Here, we take advantage of the intrinsic symmetries of 
crystal and spin structures, represented by the different symmetry groups. 
Symmetries are particularly important because they have an impact on physical 
properties — and thus on spintronics — whether they are electrical (Hall effect family), 
thermal (Nernst effect family), or optical (Kerr effect family) [K. Brading and 
E. Castellani, Symmetries in Physics: Philosophical Reflections, Cambridge 
University Press, Cambridge (2003)]. This causal link can be formulated in the 
language of Berry formalism and is related to the breaking of spatial or temporal 
inversion symmetry. One of the objectives of this chapter is to demystify the Berry 
formalism ($5.1) which is now commonly used in condensed matter physics. With 
the example of the various Hall effects, whether unquantized (85.2) or quantized 
(85.3), we will present what the Berry formalism implies for intrinsic physical effects. 
The importance of breaking spatial or temporal inversion symmetries will be 
explained and illustrated by the example of several dispersion spectra commonly 
encountered in spintronics ($5.4). In a final section (85.5), additional information 
will be provided on thermal analogs of the Hall effects. 


5.1 Berry Curvature, Berry Phase 


The Berry formalism stems from the adiabatic theorem in quantum mechanics, 
which states that a physical system remains in its instantaneous eigenstate, up to a 
phase, throughout the process of cyclic evolution, if a given perturbation is acting on 
it slowly enough. 

'The electronic transport in a material subjected to an external potential gradient 
is in fact determined from the Schrödinger equation H|,,(k)) = e,(Kk)|v, (k)) , where 
the subscript indicates the nth band. This equation returns not one but two essential 
pieces of information to describe the dynamics of electrons: 
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E the energy band dispersion, directly from the energy eigenvalues e, (k(t)) 


E the Berry curvature dispersion ©,,,(k(t)) = Vy x (w,,(k)|tVelW,,(k)), from the 

Bloch wavefunction eigenstates |y,,(k(t))) to within one phase e”»), From this 
definition, we can already infer that the Berry curvature is large wherever the rate 
of change of spin in k-space is large, for example, close to the Dirac point of a 


Dirac cone [D. Xiao et al., Rev. Mod. Phys. 82, 1959 (2010)]. 


In the above expression, k(t) is the gauge invariant crystal momentum. It 
includes the time-varying external potential in the electron frame. It stems from the 
following transformation explained in detail in the review article by D. Xiao et al., 
Rev. Mod. Phys. 82, 1959 (2010): “A uniform E means that the electric potential 
(r) varies linearly in space and breaks the translational symmetry of the crystal so 
that Bloch’s theorem cannot be applied. To avoid this difficulty, the electric field can 
enter through a uniform vector potential A(t) that changes in time.” 

We first give a set of definitions in k-space, followed by an example of a classical 
equivalent for a better understanding of this type of physics. 

It is important to understand that the Berry curvature OQ, (k) is a local quantity 
defined at all points in k-space, 


OQ, (k) = Vk x QI, oii], (K)) (87) 


It is a manifestation of the geometric properties of w,, in k-space. It is related to 
the Berry phase, which is a global quantity defined as the flux of OQ, (k) through a 
surface in k-space 


i) = d A (de f| o tok (88) 


where An(k) = (W,,(k)|iVklw,,(k)) is the Berry connection or Berry potential. The 
Berry curvature arises from a cyclic adiabatic evolution, such as electron dynamics 
in a periodic crystal structure. 

In specific cases, the Berry phase involves the integration of the Berry curvature 
over a closed surface. In these specific cases, mathematics tells us that the Berry 
phase is a multiple of 27, leading to quantum effects like those we will describe in 
§5.3. The Chern number is then defined as 


Ch = = $ Q, (Kk) d? k = integer (89) 


An equivalent of the Berry curvature, potential, and phase in real(r)-space are the 
magnetic flux density B(r), vector potential A(r) and phase g(r) = $ A(r)dr = 
Br) dr in the Aharonov-Bohm effect when electrons flow around a ring enclosing 
the magnetic flux density B(r) [C. Cohen-Tannoudji et al., Quantum Mechanics, 
Volume 1, Wiley (1997)]. In this case, the equivalent of the Chern number is the Dirac 
monopole q, = f$B(r)d?r = 2óyi, where ©) = x is the magnetic flux quan- 
tum and i is an integer. 

The way equations (87)-(89) are derived is detailed in D. Xiao et al., Rev. Mod. 


Phys. 82, 1959 (2010). We admit them here and start the reasoning from there. 
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We note that, in electron dynamics, the Berry phase relates to variations of k and is 
only dependent on the trajectory of this parameter. The Berry phase actually belongs 
to the broader class of geometric phases, grouping classical and quantum physical 
effects. It addresses the following question: what trajectory did the system follow? 

An example of a classical geometric phase is presented below, to better 
understand what a geometric phase is. 


Fic. 51 - (Left) Representation of classical parallel transport on a closed contour on a sphere, 
showing how a frame spanned by two vectors a,b can gain a geometric phase, here 8° = 5, 
which depends only on the trajectory followed. (Right) General closed contour. 


A frame (a,b) experiencing parallel transport on a closed contour C (figure 51, 
left) can gain a geometric phase °°, which depends only on the trajectory followed. 
8° corresponds to the solid angle subtented by C. For the trajectory in figure 51, 
left, we observe that q**? = §. The solid angle subtended by the contour is indeed 
g4n =}. 

For a general formalism for parallel transport on a closed contour C (figure 51, 
right), it is possible to define an initial state (a;,b;) and a final one (ay, by) for the 


frame (a,b), linked by 


ay = ajcosQ*? — bysing®° 


by = ajsing®” + bjcosq*?? (90) 


with g®° the solid angle subtended by the closed contour C. If we consider the 
following initial state for the system 

y; = ai + ibi (91) 
then the final state is 


V, — ap + iby = e^ y, (92) 


Further information on geometric phases can be found, for example, in 
J. Dalibard's lecture, www.college-de-france.fr. 
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5.2 Unquantized Hall Effects — Intrinsic AHE, SHE 


“It is now well recognized that information on the Berry curvature is essential in a 
proper description of the dynamics of Bloch electrons, which has various effects on 
transport and thermodynamic properties of crystals.” [quoted from D. Xiao et al., 
Rev. Mod. Phys. 82, 1959 (2010)]. In this section, we will detail why Berry formalism 
has played a key role in understanding intrinsic effects in spintronics, such as the 
intrinsic Hall effects. Extrinsic contributions will be briefly presented at the end. 

The Hall effects are manifested by the creation of a transverse current by 
applying a longitudinal potential gradient. Taking this statement into account, 
Ohm's law [equation (5)] reads 


Orr Ory Orz 
Je=0-E= | -omy Oy Cy | -E (93) 
Orz Oyz Ozz 


where oj with i Æ j are the transverse conductivity coefficients accounting for the 
Hall effects. In the most general case, the Hall current density Jy takes the form 


Ju —-hxE (94) 


with h = (Czy, Ors Bu) the Hall vector. We will consider that o,, = oy, = 0 in the 
next parts, for simplicity. 

Obtaining a non-zero Hall effect is subject to the breaking of a set of spatial and 
temporal reversal symmetries (PT, with P the parity or space inversion symmetry 
and 7 the time reversal symmetry) [N. Nagaosa et al., Rev. Mod. Phys. 82, 1539 
(2010)]. These asymmetries can be achieved by several means, such as the appli- 
cation of external stimuli, coupling with an arrangement of magnetic moments, and 
being intrinsically linked to the magnetic point- and spin-symmetry groups. The 
Hall effect is then described by the following transverse resistivity, made of several 
contributions that will be defined in the following paragraphs 


pa, = RN mH + Ro M + p$,- pay (95) 


Bl Ordinary or normal Hall effect (OHE), [d = R” aH. 


The OHE (figure 52, left) occurs in the presence of a perpendicular magnetic flux 
density B = uyH, when the trajectory of the electrons is deviated by the Lorentz 
force —er x B. More generally, the change of momentum reads 


hk = —cE — er x B(r) (96) 


with k the wavevector. In this case, the application of a magnetic field or flux density 
ensures 7 -breaking, and the Lorentz force connects electrons to 7 -breaking. 

Note that the acronym OHE for ordinary Hall effect is also used to refer to the 
orbital Hall effect which is of a completely different nature [D. Go et al., Phys. Rev. 
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(b) AHE (1881) *, (c) SHE( ya E, 
Dy whe 
Z í Z 


Fic. 52 — Representation of three types of Hall effects converting a longitudinal flow into a 
transverse one: (Left) the ordinary Hall effect (OHE), (Middle) the anomalous (also called 
extraordinary) Hall effect (AHE), and (Right) the spin Hall effect (SHE). In the last case, a 
net charge current is converted into a transverse spin current, but there is no net transverse 


(a) Hall (1879) 


charge current (see also $2.2). Reprinted with permission from IoP Science: C.-Z. Chang et al., 
J. Phys.: Condens. Matter 28, 123002 (2016). Copyright 2016. 


Res. 2, 033401 (2020)]. We also recall that an overdot indicates a time-derivative: 
t = Or/Ot. 

In addition to the OHE, when an electron travels through specific spin structures 
(described by the magnetic point- and spin-symmetry groups) and spin textures 
(like skyrmions, $4.7) these can exert on it an effective or fictitious magnetic field or 
flux density and an associated effective Lorentz force. It thus produces the 
anomalous (or extraordinary) Hall effect gs = R^u,M, the spontaneous anomalous 
Hall effect pa and the topological Hall effect fan depending on the nature of the 
spin structure and texture, and the interactions considered. 


E Extraordinary or anomalous Hall effect (AHE), m = R^ mM. 


The AHE (figure 52, middle) describes the behavior of ferromagnetic materials with 
a non-zero net magnetization M. A companion effect to the AHE is known as the 
spin Hall effect (SHE), in which a transverse spin current is generated under the 
application of a longitudinal electric current (figure 52, right), whether the material 
is magnetic or not. In certain cases, the AHE can be considered as being the SHE for 
a spin-polarized material, and its magnitude equals the SHE magnitude times the 
polarization. An accurate description of the intrinsic contribution to AHE and SHE 
calls for a sophisticated model. 

In the most general case, the effect can be linked to the spin-dependent electronic 
energy band structure, and formulated in the language of Berry curvature, as a geo- 
metric effect in reciprocal space. More generally, the change of position of an electron is 


Oe; (Kk) 
hok 


with ¢,(k) and Q,(k) the energy- and Berry curvature-dispersion of the nth band. 
How equation (97) is obtained is extensively described in the review article by 
D. Xiao et al., Rev. Mod. Phys. 82, 1959 (2010). The first term of equation (97) 
relating velocity to the band dispersion is the one usually encountered. The second 


Un(k) = — k x Q,(k) (97) 
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term results from internal P- or 7-symmetry breaking of the crystal,’ creating a 
non-zero Berry curvature Q,„(k), and leading to an effective Lorentz force and a 
related transverse velocity —k x O,(k) after coupling the spin and lattice via a 
spin-orbit perturbation. 


E Spontaneous anomalous Hall effect (SAHE), n 


Intrinsically, we do not need spin-orbit coupling to obtain a non-zero value for 
OQ, (Kk). The SAHE is specific to certain crystals belonging to spin-symmetry groups 
that allow internal 7-breaking despite zero net magnetization M = 0 (hence the 
name spontaneous). It is, for example, associated with a newly demonstrated 
unconventional d-wave magnetism, in analogy to d-wave superconductivity. This 
type of magnetism manifests itself as localized peaks of 2,,(k) and the associated 
electron dynamics is also described by equation (97). Those d-wave magnets open 
broad opportunities for designing and exploring unconventional quantum mag- 
netism without geometrically frustrated lattices (see THE below) and with common 
elements (large spin-orbit interactions are not needed) [H. Reichlova et al., 
arXiv:2012.15651 (2020) and L. Smejkal et al., Nat. Rev. Mater. 7, 482 (2022), see 
also figure 61]. 


Topological Hall effect (THE), po: 


The THE is specific to certain spin/moment (M) chirality, which produces an 
effective magnetic field of the type M. (om x 3M) . Those spin chiralities are, for 


example, encountered in specific antiferromagnets and in skyrmion textures. The 
phenomenon can also be formulated in the language of Berry curvature Q, [equation 
(97)] as the result of a geometric effect in real space, as opposed to the geometric 
phase in k-space discussed in the previous paragraphs. Note that, for antiferro- 
magnetic materials, the term topological used for the effect is related to the topology 
of the spin structure but does not imply a correspondence with a topological 
invariant, unlike the case of skyrmions (see $4.7 for a description of the companion 
skyrmion Hall effect). The THE is out of the scope of this book. The interested 


'Non-zero £2, requires breaking of P- or T-symmetry for the reason detailed below. 

Under  P-symmetry, we have k kk k; vu, (—k) = —v, (k); &n(—k) = en(k); 
Ie, (K) ey (Ke) 
hək ? — Wok. 
OQ, (—k) = OQ, (Kk). 

Under  7-symmetry we have k- —k kk;  w(-k)- —vw(Kk);e;(—k) = e (K); 
Oe, (k) Oc, (K) 


AB — — 45k" Combining the above transformations with equation (97) imposes that 


OQ, (—k) = —.(K). 

Hence, for a PT symmetric crystal, i.e., for a crystal not breaking either P- or 7 -symmetry, 
OQ, (k) = —Q,(k) = 0. 

We note that more subtle symmetry considerations need to be taken into account, as some 
operations, like translation can restore symmetry, resulting in O,(k) = 0 despite initial P- or 
7 -breaking (figure 61). 


Combining the above transformations with equation (97) imposes that 
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reader can find more information in Topology in Magnetism, J. Zang et al. (eds), 
Springer Nature (2018). 

We next show how the contribution of Berry curvature to v,(k) [equation (97)] 
connects with the Hall transverse conductivity o,, [equations (93) and (94)]. From 
Boltzmann's semiclassical transport theory [C. Kittel, Introduction to Solid State 
Physics, John Wiley & Sons (2004)], the charge current is 

dk 
L=- [ states) C (98) 


where d is the dimension of the system, BZ is the Brillouin zone (the integral is over 
the d dimensions) and g(én(k)) = f(én(k)) + óg(e,(k)) the distribution function. We 
remind that this function accounts for the number of electrons with momentum k in 
the infinitesimal element d^k whose position r is located in the infinitesimal element 
d'r. f(en(k)) is the Fermi-Dirac distribution at equilibrium, and dg(¢,(k)) accounts 
for the deviation of the distribution from its equilibrium. 

Combining equations (97) and (98), with k=- ; E, it is possible to show that 
the total charge current takes the form 


Oe, (k) d 
Je= 23 TE yO óg(e, (k)) d? k 


+ - x DE |, «(0 f (e (Kk)) d? k 


(99) 


The first term in equation (99) is the common one resulting directly from the 
energy dispersion. The second term represents the intrinsic contribution to the 
transverse conductivity. From J, = F - E this contribution is expressed as 


e 1 
Ory = he Gael f «€ f (e, (K)) d* k (100) 


The above equation describes the unquantized version of the Hall effect [see also 
equations (101) and (102)]. 

We note that, per definition, for a three-dimensional system, when d — 3 in 
equation (100), c; is in S.m^!, whereas, for a two-dimensional system, as in $5.3, 
when d = 2, Oy is in S, unit of conductance. 

It is clear at this point that the electronic fingerprint of a material must contain 
two parts: the dispersion of the energy bands e,(k) and that of the Berry curvature 
or total Berry curvature Q(k) = $7, f (e, (k))Q, (k), as shown for bcc Fe in figure 53. 
'The peaks and valleys in the distribution of the total Berry curvature are due to 
pairs of unequally populated spin-orbit coupled bands, in a small k-interval, i.e., 
wherever the rate of change of spin in k-space is large. 

Fcc Pt is another prototypical material that illustrates the importance of the 
energy band topology and thus the need to provide the Berry curvature fingerprint. 
From figure 54, we observe that the intrinsic spin Hall conductivity c7, peaks just at 


the Fermi level, at e — £p = 0. This makes Pt a material of choice for charge current 
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T(000) H(100) 
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Fermi loop (first term) 7 0 -4 
Berry curvature 753 -2203 477 
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Fic. 53 — (Top) Table of the anomalous Hall conductivity czy, computed in several ways, for 
bcc Fe, fcc Ni, and hcp Co. The third row shows the results obtained from calculations using 
the theoretical framework of Berry curvature. The other approximations are beyond the scope 
of the present book. The last row shows experimental data. The disagreement between theory 
and experiment for fcc Ni is due to the moderate onsite electron-electron correlations, which 
must be taken into account with more precision. Reprinted with permission from the 
American Physical Society: X. Wang et al., Phys. Rev. B 16, 195109 (2007). Copyright 2007. 
(Bottom Left) Electronic band structure and z-component of the total Berry curvature 
Q^(k) = 35, f(e,(K))Q; (k) along symmetry lines of the reciprocal space, calculated for bulk 
Fe. (Bottom Right) Corresponding Fermi surface in the (010) plane, represented by solid lines, 
and Berry curvature in atomic units, represented by the color map. Atomic units mean 
x (n/ ay? with a the lattice parameter, where the proportionality factor depends on the 
number of k-points computed. For the calculations, the spin-orbit strength was 5.1 mRy 
(1 Ry ~ 13.6 eV), and spin-orbit interactions were thus treated in a perturbative manner. 
Reprinted with permission from the American Physical Society: Y. Yao et al., Phys. Rev. Lett. 
92, 037204 (2004). Copyright 2004. 


to spin current interconversion (figure 52, right, see also $4.4 and exercise 6.6). 
The large value of c7, is attributed to the double degeneracies at the L and X points 
of high-symmetry, near the Fermi level, which are present even without spin-orbit 
coupling. From figure 54, it can be seen that when ep is lowered from its initial 
position, 67, reduces, then changes sign and peaks again at about —4.5 eV. 

Let us now crudely consider that the Berry curvature fingerprint of the Pt 
material holds for all materials in the 5d row (Ta, W, ... , Au). The change from Au, 
Pt to W, Ta reduces the band filling and is equivalent to a lowering of £p. This rough 
approximation helps understand the dependence of o7, on band filling, and in 
particular, the change in sign between o7, (Au, Pt) and c;,(Ta, W) (figure 55). 
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Fic. 54 — (Left) Electronic band structure along symmetry lines of the reciprocal space, 
calculated for bulk fec Pt, when spin-orbit coupling (SOC) is turned off and on. (Right) Spin 
Hall conductivity (Tiy) calculated for the same material, using the theoretical framework of 
the Berry formalism. Reprinted with permission from the American Physical Society: 
G. Y. Guo et al., Phys. Rev. Lett. 100, 096401 (2008). Copyright 2008. 


E Note about extrinsic contributions 


The first term of equation (99) accounts for the longitudinal current. It may also 
contain transverse components in the presence of skew (or Mott) and side jump 
scattering on defects (figure 56), which impacts the non-equilibrium part óg(e,(k)) 
of the distribution function g(en(k)) [N. Nagaosa et al., Rev. Mod. Phys. 82, 1539 
(2010)]. Simply put, because the skew and side jump scattering mechanisms involve 
a one-step and a two-steps scattering process, respectively, their contributions to the 
total Hall resistivity um are linearly and quadratically dependent on the longitu- 


skew 


dinal resistivity, respectively: pS" x p,, and päde-iua? ox o?.. In addition, because 


the intrinsic Berry curvature contribution is Dub qd: of scattering, we have 


gintrinsic — — constant, and therefore oo x p from ginirinsic — — pas We 


zy ry 

note that the relation between pne and yes is dup obtained from inversion 
of the resistivity matrix? Ggusidbring that p,, > JR i.e., when the determinant 
reduces to p?.. From a simple scaling relation, the total transverse resistivity reads 
piel = = ap,, + bp?,, where the coefficient a contains the skew scattering contribution 


and the coefficient b the side-jump scattering and intrinsic contributions. 


-1 
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The intrinsic spin Hall conductivity 


SOI 
of Ad and Sd metals Metals Structure Electron number (Ry) 
Nb bee 5 (Ad'55!) 0.006 

Mo bcc 6 (4d55!) 0.007 

Te hep 7 (4d°5s') 0.009 

Ru hep 8 (4d'5s!) 0.01 

ay Rh fee 9 (Ad*5s!) 0.011 
g Pd fcc 10 (449559) 0.013 
E Ag fcc 11 (4455!) 0.019 
a Ta bcc 5 (5d'652) 0.023 
p w bee 6 (Sd*6s?) 0.027 
c Re hep 7 (Sd°6s*) 0.025 
Oo Os hep 8 (Sd°6s*) 0.025 
T Ir fee 9 (Sd'65) 0.025 
Pt fec 10 (Sd?6s!) 0.03 


Au fec 11 (5d6s') 0.03 


"9 6 7 n 8 9 10 1l 
4d Nb Mo Tc Ru Rh Pd Ag 
5d Ta W Re Os Ir Pt Au 


Energy (eV) 


Fic. 55 — (Left) Spin Hall conductivity (SHC) at the Fermi level versus electron number (n), 
calculated for the 4d and 5d transition metals, using the theoretical framework of Berry 
curvature; (Right Top) table listing several parameters, for the 4d and 5d transition metals. 
SOI stands for spin-orbit interaction (1 Ry ~ 13.6 eV). Reprinted with permission from the 
American Physical Society: T. Tanaka et al., Phys. Rev. B 17, 165117 (2008). Copyright 2008. 
(Right Bottom) Spin Hall conductivity (oj) calculated for bulk fcc Pt, using the theoretical 
framework of Berry curvature. The vertical lines indicate the spin Hall conductivity at the 


Fermi level by considering that the calculation of o2, vs. energy in Pt stands for Ta, W, and 


r4 
zy 
Au, and that band filling moves the Fermi level. Adapted with permission from the American 


Physical Society: G. Y. Guo et al., Phys. Rev. Lett. 100, 096401 (2008). Copyright 2008. 


Experimentally, a and b are obtained by varying the contribution from impurity 
scattering via changing the impurity content or temperature. The exact form of the 
scaling relation is actually more complex. It naturally involves residual 
(impurity-related) resistivity p,, r-o The fact that the multiple scattering mecha- 
nisms compete with each other must be considered [D. Hou et al., Phys. Rev. Lett. 
114, 217203 (2015)]. 

The mathematical description of extrinsic mechanisms is tedious. The models 
that describe these mechanisms go beyond the simplified scope of this book and are 
described in detail in, for example, the review article of N. Nagaosa et al., Rev. Mod. 
Phys. 82, 1539 (2010). 

We note that the spin Hall angle is defined as OsHE = p,,/p,,. Considering the 


intrinsic 


intrinsic mechanism only, we thus obtain  Üggg = em Po, from 


gi rimsie — pinimmsic / p^. In this case, the Hall angle is directly proportional to the 


longitudinal resistivity. This explains why the Hall angle for Au is much smaller than 
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Side jump 


The electron velocity is deflected in opposite directions by the opposite 
electric fields experienced upon approaching and leaving an impurity. 
The time-integrated velocity deflection is the side jump. 


Skew scattering 


Asymmetric scattering due to 
the effective spin-orbit coupling 
of the electron or the impurity. 


Fic. 56 — Schematics of two extrinsic mechanisms for the spin Hall effect: side jump and skew 
scattering. Adapted with permission from the American Physical Society: N. Nagaosa et al., 
Rev. Mod. Phys. 82, 1539 (2010). Copyright 2010. 


for Pt, because Au is one order of magnitude less resistive than Pt, although the 
intrinsic transverse spin conductivity of Au and Pt are similar (figure 55). The factor 
of merit for spin-charge interconversion is therefore Osyp lsp, because I, ox 1/p,,,. 


5.3 Quantum Hall Effect - QHE 


In the previous section, we have seen the most general case of the unquantized version 
of the Hall effect. Quantization comes into play for certain types of energy band 
structures, like in some 2D electron gas and band insulators under large magnetic 
fields, especially when the energy levels or the bands become discrete Landau levels 
due to quantization of the cyclotron orbits of the particles. How the quantum 
mechanics of particles moving in a magnetic field connects to Landau levels with a 
finite density of states is, for example, demonstrated in D. Tong, arXiv:1606.06687. 
In such conditions, it is also possible to fold the Brillouin zone (BZ) to a torus 
[J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) (2007)], and projection or 
contraction onto a closed sphere under symmetry operations shows that points in the 
northern hemisphere are conjugate under 7 to those in the southern hemisphere. 
Therefore, the Berry phase or the integration of the Berry curvature on the closed BZ 
[equation (89)] becomes a multiple of 2z. Equation (100) simplifies as 


e 1 2 
Ory = aoa: fe k (101) 


Combining this last equation and equation (89) leads to 
2 
Ory = 22 C, (102) 


where C, is the Chern number. 


90 The Basics of Electron Transport in Spintronics 


We recall that, per definition [equation (100)], for a two-dimensional system, as 
here, G, is in S, unit of conductance. 

'The Hall effect then becomes quantized and is called the quantum Hall effect 
(QHE) [H. Weng et al., A-APPS Bulletin 23, 3 (2013)]. The quantum Hall con- 
ductivity can be considered as a counter of the number of fully filled Landau levels. 

The QHE was discovered by K. von Klitzing et al. [Phys. Rev. Lett. 45, 494 (1980)] 
in silicon-based 2D electron gas MOSFET (metal-oxide-semiconductor field-effect 
transistor) stacks. Under a strong magnetic field, the 2D electron gas bands are 
quantized into discrete (yet with finite width) Landau levels, which move up in energy 
when increasing the field. As a result, the number offilled Laudau levels is reduced with 
increasing field, every time a Landau level crosses the Fermi level. Consequently oj, 
(Pry) reduces (increases) with field in a series of discrete jumps (figure 57). 


10 


V. (mV) 


Vui(mV) 


Magnetic field (T) 
Von Klitzing constant : R, = h/e? 


Fic. 57 — Quantum Hall effect demonstrated through quantized steps in the magnetic 
B-field-dependence of the transverse Hall voltage Vy. Each step corresponds to an incremental 
change in transverse conductivity: Cry = mc where i is an integer [equation (102)]. Reprinted 
with permission from Springer Nature: K. von Klitzing et al., Nat. Rev. Phys. 2, 397 (2020). 


Copyright 2020. 


The resistance Rg = E was subsequently named the von Klitzing constant. It 
depends only on physical constants and is more stable and reproducible than any 
other resistance: Ry = 25812.807 Q. For that reason, the QHE was essential for the 
redefinition of the kilogram in the international systems of units in 2019 (figure 58). 
Nowadays, SiC-substrate/graphene stacks are used as the QHE standard. The von 
Klitzing constant Rx relates to the conductance quantum Gy as follows: 
Go — gc x As an aside, we note that the Josephson effect [B. D. Josephson, 
Phys. Lett. 1, 251 (1962)] in alternating current allows measurement of a similar 


quantity: Ky = 2e, named the Josephson constant. 
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Fic. 58 — Diagrams of the base units of the old (outside) and new (inside) international 
systems. The new one, since 2019, is based on physical constants. From BIPM (bureau 
international des poids et mesures), A concise summary of the International System of Units, 
SI, https://www.bipm.org/en/home — CC BY License. 


It is important to note that the QHE is one of the main manifestations of 
topology in condensed matter physics. For its discovery, K. von Klitzing was 
awarded the Nobel Prize in Physics in 1985. He opened the way to many other 
important related discoveries, including other types of quantum Hall effect [H. Weng 
et al., A-APPS Bulletin 23, 3 (2013)] for which other Nobel prizes have been 
awarded, for: “the discovery of a new form of quantum fluid with fractionally 
charged excitations” in 1998 (Nobels H. Stórmer, D. Tsui, and R. Laughlin), 
“groundbreaking experiments regarding the two-dimensional material graphene” in 
2010 (Nobels A. Geim and K. Novoselov), and “theoretical discoveries of topological 
phase transitions and topological phases of matter” in 2016 (Nobels D. J. Thouless, 
F. D. M. Haldane, and J. M. Kosterlitz). 


5.4 Parity and Time Reversal Symmetries — P, T 


We have seen in §5.2 that obtaining a non-zero Berry curvature and subsequently 
non-zero intrinsic Hall effects require breaking of P- or 7 -symmetry. In figure 59, we 
give some examples of how P- or 7-symmetries alter the spin-dependent energy 
band dispersions. We next develop how some of those symmetries connect with the 
Berry curvature. 


B Kramers — degenerate 


For free electrons (figure 59, top left), the band structures for majority moment 
(T)-spin(]j) electrons and minority moment(|)-spin(f) electrons are degenerate, 
and no-Hall effect is expected. The corresponding Hamiltonian contains the kinetic 
energy for free electrons and reads 


Ne 
H = 


= 1 
jn. (103) 
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Kramers - degenerate Rashba - P-breaking 
(nm 
> > 
E 2 
o o 
c = 
LU LU 
k k 
Zeeman - J-breaking Rashba + Zeeman - P+J-breaking 


Energy 


Fic. 59 — Band structures for majority moment(1)-spin(|)) electrons and minority moment 
(l)-spin(f) electrons, corresponding to several parity and time symmetries, which are 


reflected in the Hamiltonian H involved. 


E Zeeman — 7-breaking 


For a ferromagnet (figure 59, bottom left), the band structures for majority moment 
(T)-spin(]) and minority moment(]|)-spin(f) electrons are split in energy (82.3), 
leading to a non-zero magnetization and a Hall effect upon spin-orbit coupling. The 
corresponding Hamiltonian includes the Zeeman term 


— - AG; (104) 


E Rashba - 7-breaking 


An example of Rashba parity breaking is when symmetry is extrinsically broken at 
an interface (figure 59, top right). The band structures for majority moment(T)-spin 
(4) and minority moment(|)-spin(f) electrons are thus shifted in momentum, 
leading to a Hall effect upon spin-orbit coupling. The spin angular momentum is 
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then locked in a direction perpendicular to the linear momentum k. The corre- 
sponding Hamiltonian includes the Rashba term (see also figure 49, left and middle 
and corresponding text) 

Wk 


H = — +ap(k x a): Z (105) 
2m, 


Contributions triggered by spin-orbit coupling are called ‘relativistic’ because 
OXE), Conversely, the contribu- 
tions triggered intrinsically with no need of spin-orbit coupling are called 
non-relativistic — although spin-orbit coupling is needed to couple spin and lattice in 
a final step. We note that the magnitude of relativistic contributions is limited by 
the ratio between electron velocity and speed of light, — 10 ^ in typical metals. 
'These contributions thus require a large spin-orbit interaction to be significant. By 
contrast, non-relativistic contributions can be of the order of a few eV, at least one 
order of magnitude larger [Topology in Magnetism, J. Zang et al. (eds), Springer 
Nature (2018)]. 


they involve a Lorentz transformation, HoH 5g = — 


Bi Spin-polarized 2D electron gas with Rashba spin-orbit coupling 
— P- and T-breaking 


A prototypical example illustrating Berry curvature physics is the case of a 
spin-polarized 2D electron gas with Rashba spin-orbit coupling (figure 59, bottom 
right) [S. Onoda et al., Phys. Rev. Lett. 97, 126602 (2006) and D. Xiao et al., Rev. 
Mod. Phys. 82, 1959 (2010). The combination of Rashba P- and Zeeman 
T -breaking results in a shift in the momentum of the band structures for majority 
moment(T)-spin(i.) electrons and minority moment(|)-spin(f) electrons, accom- 
panied by a gap opening at k — 0. The spin angular momentum is locked on the 
momentum k. The Hamiltonian is 


VE 
H = zp t+ar(k x o) -Z — Ao; (106) 


The associated eigenvalues (figure 60, top left) can be expressed as (exercise 6.10) 


sa(k) = 1 + ag? +A? = solk) + (A)| (107) 


| 2m. 


The corresponding eigenvectors are (exercise 6.10) 


A 
Ib(&)| 


|J) 4- ie? [14 ipao (108) 


where $ is defined from k,+ iky = ke”, and |J4) = xe"(s) and |y) = 


Z eer ( 1) are the eigenvectors for free electrons. 
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From equations (87) and (108), it is then possible to show that the Berry cur- 
vature of the two subbands is (exercise 6.10) 


arg ^ 


2 (ue? m 3/2 (109) 


The Berry curvature is non-zero around the gap, near k = 0 (figure 60, bottom 
left) if both Rashba (ap # 0) P- and Zeeman (A Æ 0) T- breaking take place. 


Energy 


2 


e 2 d?k 
Hall conductivity x] 2,0500 Cr)? 


Spin polarized 2D electron gas 
with Rashba SO coupling 


h2 2 
fe te aa aes 


Berry curvature 


2(ar?k? + A2) 


Fic. 60 — (Left) Band structure and Berry curvature for a spin-polarized 2D electron gas with 
Rashba spin-orbit (SO) coupling, described by the Hamiltonian (H). Note: beware the color 
code, when red/blue is commonly used to distinguish majority moment(T)-spin(4J) electrons 
and minority moment(|)—spin({}) electrons bands (see figure 59, bottom right panel), it is also 
commonly used to distinguish the upper/lower band, no matter spin (like here). (Right) 
Corresponding Hall conductivity obtained from integration of the Berry curvature over the 
wavevector k for given energies. 


The Hall conductivity o,, is further calculated numerically based on equation 
(100). From figure 60, right, we observe that for e< — A, o,, first increases when 
increasing the energy, as the contribution to the Berry curvature of the lower band 
Q (k) is integrated. In the gap, for CA < e < A, Gry saturates as the lower band has 
contributed fully. Finally, above the gap, for € > A, oz, reduces as both contributions 
Q? (k) +94 (k) are now integrated, and the contribution to the Berry curvature of 
the upper band opposes that of the lower band. 
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Symmetry of the Antiferromagnetism Hall effect 
spin structure n-m,-m; h = (Czy: Oxz: 0.) 


tT (n) = -(-n) =n tT (h) = +(-h) = -h 


= = Even zz Odd 
E Compatible Incompatible 
i= 
os 
£2 PT (n) = —(-n) =n PT (h) = +(—h) = -h 
LES = Even = Odd 
a E Compatible Incompatible 
oc 
E E ta(n) =n = ty2€2y(h) = 
E — t6 (n)--(-n)-n Even - Compatible along y RuO,, CrSb, 
> E Mni 
NIU Æ ty€r,T(h) = -(-h) = -h Dg 
oR Odd — Incompatible along x 
* 
= typM,R;(n) =n tM. h)-h 
5 E x ead ) a5 a m Epi-Mn.Sis 
2 J Ala h l 
S z Compatible Compatible (hexagonal) 
BS 
zd 
3t 
a 
2 nN 
3 
Š$ typM,RsT(n) =n typM,RsT (h) = -h : 
E = Eon zz Odd Bulk-MnsSis 
E E - (orthorhombic) 
Rü Compatible Incompatible 
[:] 
N 


Fic. 61 — Illustration of a symmetry analysis to establish which of the antiferromagnetic 
structures shown on the left column may host the spontaneous anomalous Hall effect. Here, 
m; and m» are the sublattice magnetizations, normalized, and n is the Néel vector, which 
defines the magnetic order of an antiferromagnet. The tı /2 half-unit cell translation, C2,,y) 
2-fold rotation-axis symmetries, M, mirror symmetry and Rs spin rotation allow describing 
the antiferromagnetic structures. Applying the corresponding symmetry combinations to the 
Néel vector n is naturally compatible with the proposed structures. However, applying these 
symmetry combinations to the Hall vector h is compatible in some structures only, meaning 
that only these structures may host both antiferromagnetism and the spontaneous anomalous 
Hall effect. Note that M, and Rg are unlocked, and description via spin groups is therefore 
relevant for the two corresponding cases (see also figure 62). The illustration is made after the 
information contained in L. Smejkal et al., Nat. Rev. Mater. 7, 482 (2022) and H. Reichlova 
et al., arXiv:2012.15651. 


'To conclude this section, we note that the Curie-Neumann principle of symmetry 
states: ^when certain causes produce certain effects, the elements of symmetry of the 
causes must be found in the effects produced". Based on a symmetry analysis, it is 
therefore possible — and now systematically done in practice — to establish which type 
of symmetry is compatible with one physical effect, such as the Hall effect, and which 
type of related crystal can a priori host the effect (figure 61). The reciprocal principle 
is not true as the effects produced are often more symmetrical than the causes. 
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We also remark that a thorough symmetry analysis calls for the use of group theory, 
whether the relevant group is the crystallographic point group (no magnetism), the 
magnetic point group (in which magnetic and crystal rotation operations are locked; 
it corresponds to the case of large spin-orbit coupling) or the spin group (in which 
magnetic and crystal rotation operations are partially independent of each other 
(figure 62) [H. Reichlova et al., arXiv:2012.15651 (2020), L. Smejkal et al., Nat. Rev. 
Mater. 7, 482 (2022), P. Liu et al., Phys. Rev. X 12, 021016 (2022)]. 


Magnetic order 


Nonmagnetic 
w/o SOC Oo 


~ fF 

H= 2m * V 

Space group 
{Cn (0) |7} 


Nonmagnetic D Ò Magnetic 


SOC SOC 


2 = E 1 
Re E V) — Ove) x p)-@ A= E + VP) + 2 (VG) x B) - 8 (9-8 
Double space group Magnetic double space group 
{Un(8)||Cn(8)|z} {Un(8)|1Cn(8) |} 


Magnetic order 


Fic. 62 — Illustration describing the symmetry of solids. S(F) and V(F) in the Hamiltonian 7, 
stand for exchange field due to the distribution of local magnetic moments, and crystal 
potential energy, respectively. 6 is the Pauli operator associated with the electron spin and p 
the momentum operator associated with the linear momentum. Spatial and spin rotations are 
described by C, (0) and U,(,(0, $), where n and m denote the rotation axes, and 0 and $ are 
the rotation angles. The notations w and w/o SOC correspond to with and without spin-orbit 
coupling. w SOC (Bottom panels), the spatial and spin rotations are locked. w/o SOC (Top 
panels), the spin and spatial rotations are completely or partially decoupled. In nonmagnetic 
cases (Left panels) spatial rotation only is considered or an unconstrained spin rotation is 
added. In magnetic cases (Right panels), spin rotation is constrained by the magnetic order of 
the system. From P. Liu et al., Phys. Rev. X 12, 021016 (2022) — CC BY license. 


5.5 Thermal Nernst Counterparts 


Electrical and thermal transport are closely related. Therefore, there are thermal 
counterparts to electrical effects in spintronics. Analogously to Ohm's law 
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J.=—@-V¢ with F the electric conductivity tensor, and $ the electric potential 
[equation (5)], which describes the magnitude of the induced electric current in 
response to a potential gradient, the magnitude of the induced thermoelectric cur- 
rent in response to a temperature gradient is as follows 


Je=-a:-VT (110) 


with x the thermoelectric conductivity tensor. This tensor takes into account, among 
others, the Seebeck effect in the diagonal terms, with æ = ø- S, where S are the 
Seebeck coefficients, and the Nernst effect in the off-diagonal terms. 

Analogous to the Hall effect, the Nernst effect accounts for the creation of a 
transverse voltage in response to a longitudinal thermal gradient. The thermoelec- 
tric transverse resistivity o9, (with on = p%,/ p2,?) can be described similarly 
[equation (95)] by a sum of normal, anomalous, spontaneous anomalous and topo- 
logical contributions as 


pi,— RN mH + Re^ uM + pu du (111) 


In the most general case, the thermoelectric and electric coefficients are related 
by the Mott relation [equation (112)]. It has been shown that the validity of these 
relations extends to intrinsic contributions based on the Berry curvature formalism 
[D. Xiao et aL, Phys. Rev. Lett. 97, 026603 (2006)]. Therefore, the intrinsic 
anomalous Nernst conductivity is also related to the intrinsic anomalous Hall con- 
ductivity, calculated in $5.2, via the following Mott relation 


9, P 99a er) 


Gs, de (112) 


ay = 


2 ] à RT a 
where Go = 27 = E is the electrical conductance quantum and gj = 2 37 is the 


thermal conductance quantum, to which quantized thermoelectric effects relate 
[J. Noky et al., J. Phys. Commun. 5, 045007 (2021)], in a manner analogous to the 
link between Gp and quantized electric effects like the quantum Hall effect (85.3). 
From equation (112), we see that «, is related to the derivative of os, with 
respect to energy, and thus to the derivative of the DOS with respect to energy, in a 


similar way to other thermoelectric properties. Therefore, «sy will be maximized as 
065, 


the edges of the bands are approached (figure 60, from >"). At the microscopic 
scale, the application of a temperature gradient populates the higher energy states 
and depopulates the lower energy states, on the hot side. The system returns to 
equilibrium by diffusion. A net current is created if and only if the diffusion of the 
high-energy and low-energy states is asymmetric, hence the derivative in the Mott 
relation. 

More generally, the links between thermoelectric and electrical coefficients are 
valid whether the contributions are intrinsic or extrinsic, and whether the effect 
involves diagonal or off-diagonal terms. For example, how the spin-dependent Seebeck 
effect alters the drift-diffusion equation for spin accumulation is shown in exercise 6.3. 
All this contributes to the richness of the cross-effects that exist (figure 63), 
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x 


Fic. 63 - Illustration of the wide range of charge, spin, and heat-related effects and how these 
effects are intertwined. In this illustration, £” is a tensor that connects a charge (i = c), spin 
(i= s), or heat (i = q) current to a force of an electric Vf (j = c), spin Vu, (j = s), or heat 
VT (j = q) nature. The number on the bottom right of each panel indicates the order of the 
corresponding tensor. In this book, we discussed the AHE, ANE, ISHE, and SHE. Courtesy of 
H. Ebert. See also M. Seeman et al., Phys. Rev. B 92, 155138 (2015) for a full description. 


the corollary being that these effects can be difficult to distinguish from each other in 
experiments [S. R. Boona et al., Energy Environ. Sci. 7,885 (2014); K.-i. Uchida et al., 
J. Phys. Soc. Jpn 90, 122001 (2021); M. Seeman et al., Phys. Rev. B92,155138 (2015)]. 
'To complete this comment, in the same way that it is possible to pump spins and then 
act on the magnetization via the application of an electric current through a potential 
gradient, it is also possible to induce similar effects using a temperature gradient, and 
cross-check the results to confirm the findings. 
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Summary 


In this chapter, we have seen that Berry formalism, which derives from the 
adiabatic theorem of quantum mechanics, is appropriate to describe the cyclic 
adiabatic evolution of Bloch electrons in a periodic crystal structure. In addition 
to the conventional energy band dispersion ¢,,(k), the Berry curvature dispersion 
Q,,(k) is an essential ingredient for a correct description of electron dynamics. As a 
thumb rule, the Berry curvature is large, wherever the rate of change of spin in 
k-space is large. From its impact on the electron velocity in k-space 


v,(k) = d — k x Q, (Kk), the Berry curvature gives rise to intrinsic spintronics 
effects, like the AHE, the SHE, the THE, and Nernst thermal counterparts. It also 
leads to quantized physical effects, like the QHE, under specific dimensionality 
and cyclic conditions. P- or T-symmetry breaking is essential for the Berry 
curvature to be non-zero, therefore underlying the importance of considerations 
on symmetry groups. Finally, thermal and electric effects are related via the Mott 


x OC ry EF : k , g 
relation % = — 2e 9 Malte ) which can be used in spintronics. 


Chapter 6 
Exercises and Solutions 


'This chapter provides a series of ten comprehensive exercises with solutions. They 
are specifically designed to illustrate the ideas in the previous chapters and lead to 
other spintronic effects based on these ideas. 

'The exercises deal with: 


Bl Anisotropic magnetoresistance (AMR) (chapter 2) 

ll Domain wall anisotropic magnetoresistance (DWAMR) (chapter 3) 

Bi Drift-diffusion equation for spin accumulation (u,) (chapters 3-5) 

li Spin conductivity mismatch (CPP-GMR) (chapter 3) 

E Intra-band scattering — intrinsic damping (25) (chapters 2 and 4) 

Bi Spin pumping (SP) and inverse spin Hall effect (ISHE) (chapters 4 and 5) 

ll Spin pumping (SP) — additional damping («,) (chapter 4) 

Bi Spin Hall magnetoresistance (SMR) (chapters 2-4) 

B Harmonic analysis of the anomalous Hall voltage (V$, vy (chapters 4 and 5) 
E Intrinsic anomalous Hall and Nernst effects (AHE, ANE) (chapter 5) 


We recall that a List of symbols and units, as well as the corresponding formulas, 
is provided on page 150. When there are several possible definitions in the literature, 
the definition used in this book is made explicit, and the differences in prefactor 
depending on the articles and books consulted are also made explicit. 


6.1 Anisotropic Magnetoresistance (AMR) - chapter 2 


Bi Introduction/Reminder 
The anisotropic magnetoresistance (AMR) effect (82.6) refers to the dependence of 
the electrical resistivity p on the relative angle 0 between the applied electrical 
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current I and the magnetization M of a magnetic material [W. Thomson, Proc. Roy. 
Soc. 8, 546 (1857); T. McGuire and R. Potter, IEEE Trans.Magn. 11, 1018 (1975)]. 
It is a bulk property caused by anisotropic mixing of majority moment(T1)-spin(l)) 
electrons and minority moment(|)—spin({}) electrons conduction bands, induced by 
the spin-orbit interaction. 

In the reference frame set by M (figure 64), the electric field E and the current 


density Je are linked by the resistivity tensor p 


E — pJ. 
een 0" 
E, 0 PL Jı 
with pj and p, , the resistivities for Je || M and Je L M, respectively (82.6). 


E Questions 


(a) Find the expressions of E, and E, vs. J», Jy, pj, p, and 0, in the experimental 
reference frame defined by (X, y), (figure 64). 


y 


525 


Fic. 64 — Illustration showing the two reference frames considered. 


(b) In practice, how would you proceed experimentally to measure py and p,? 
What is the order of magnitude of the values of field, current, voltage, etc. that 
you think you would use or measure? 

(c) Experimental measurements of Ni films at room temperature returned p, — 8.2 
uQ.cm and p, = 8 uQ.cm. Calculate the AMR ratio ^p = um for Ni. Can you 
think of any practical use of the AMR effect? 

(d) We now consider a ‘Union Jack’ device shown in figure 65, a magnetization 
pointing along X, a film of thickness d, and a time-dependent square wave 
current densityJ;;(1). Plot the time(t)-variation of: Vis, and Va; for Jis; Vga, 
and Vos for Jg4; Vz4, and Vig for J73; Vg», and Vg, for Jeo. 

(e) An antiferromagnet is used instead of a ferromagnet. This antiferromagnet has 
two collinear sublattices with magnetizations pointing towards opposite direc- 
tions (M; = —M;»). The total magnetization is M = M, + M»; = 0. Do you think 
that using AMR is appropriate to characterize this type of magnetic material? 
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Fic. 65 — Illustration of the ‘Union Jack’ device considered. 


E Solutions 


(a) A change of basis can be done by using the appropriate rotation matrix 
E cosÜ sin E. J cos? sin J 
|) — x | \ — z 
(2) B ee m tu and (3) B Ge an @ pd 
Combining equations (113) and (114) gives 
( cosÜ A @ " E 0 ) ( cosÜ a) & 
—sind cosO/\E,/ XO p,/\—sin@ cos0)\ Jy /’ 
(2) 7 ( cosÜ T (^ 0 ) ( cosÜ =) (5) 
E, © V-—sinÜ cosÜ 0 p, —sinÜ  cosÜ Jy)? 
es 1 ee ed E 0 ) ( cosÜ sind ) ( 5l 
E, ~ cos?0--sin?0 \ sinÜ cos 0 p, —sinÜ  cosÜ Jd 
(5) E cosO py —sindp , ( cosÜ oe ( 7) 
E, i sinp — cos0p, —sinÜ cos da 


a cos?0p, 4- sin?0p costsinü('p, —p. ) ( 7) 
cosÜsinO (o, -p1) sin’ Op) + cos?0p , 


(&) i e - (p= p. )si?0 ^ cosósint(p,—p.) t3 "TD 


cosósinü('p, —p. ) pL,—lp,— p; sin?o Jy 
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For J, = 0, we obtain 


E, = |^ = (o, = p.) sin?0]| Ja, 
EQ-— (6-01) cosOsinO Jz. 


Note that, in this case, a transversal voltage drop x E, is obtained when 0 4 0 
and 0 4 2/2. By analogy, this effect is called the planar Hall effect. It is note- 
worthy that the terminology planar Hall effect is misleading because the phe- 
nomenon at play is unrelated to the Hall effects. The terminology transverse 
AMR is sometimes preferred. 

Possible option 1: use 4-point resistance measurements and a known geometry 
[I. Miccoli et al., J. Phys. Cond. Mat. 27, 223201 (2015)]. 

Set J, # 0 and I, = 0. 

Monitor Vz, to get R;. 

Use a magnet to apply an external magnetic field at 0 = 0 and get pj; and at 
0 = 90° and get p,. 

Possible option 2: use 4-point resistance measurements and a known geometry. 
Use a magnet to apply an external magnetic field and set 0 = 0. Set I, Z 0 and 
I, =0 and monitor V;, to get R, and then pj. Set I; — 0 and I, #0 and 
monitor V,, to get Ry and then p,. 

Orders of magnitude: mm device, mA, tenth of Ohms, mV. 

For Ni, Ae ~ 2.5% at room temperature. 

The AMR effect was used in the first generation of sensors, like read heads in 
hard disk drives. Nowadays, most sensors are based on tunnel magnetoresis- 
tance with orbital filtering. Typical values of magnetoresistance are now greater 
than 100%. The temperature dependence of AMR sensors can still be a plus for 
some applications. 

From equation (115) we can deduce the following: 

For Jj5, we have 0 = 0 


Eis = py J15 


Visa 


V. 
37 a 
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For Jg4, we have 0 = —45? 
Eg, = (o, +p.) Jg4/2 


Ex = — (oi-0.) Ja4/2 


Vas 


V26 


For Jz3, we have 0 = —90° 
Er = p, J73/d 


Fis = 0 
For Jg», we have 0 = —135° 


Een = (o, F p.) Js4/2 


Es4 = (o)-p.) J62/2 


(e) The AMR effect is even in magnetization, i.e., it is invariant on magnetization 
reversal. It can be shown that 2e x (M Je)”. The AMR responses from the two 
sublatices add up. Using the AMR effect is hence appropriate for 
characterizing/detecting a collinear antiferromagnet [P. Wadley et al., Science 
351, 587 (2016)]. 

3 


2 4 
M, 
1 v —$ 5 
M, 
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6.2 Domain Wall Anisotropic Magnetoresistance 
(DWAMR) - chapter 3 


E Introduction/Reminder 
Several types of magnetic domain walls (DWs) in magnetic wires (83.5) are depicted 
in figure 66. An electrical current (I) flows along the wires. 


(1) Vortex wall (3) Bloch wall 


Hor MNM 


(2) Transverse wall Å (4) Néel wall 


Fic. 66 — Representation of 4 types of magnetic domain walls. Adapted with permission from 
IoP Science: R. L. Stamps et al., J. Phys. D: Appl. Phys. 47, 333001 (2014). Copyright 2014. 


E Questions 


(a) Which of these DWs display a signal of anisotropic magnetoresistance ĉe 
compared to the uniformly magnetized wire, and why? We recall that 
Ap — Pwithpw—Pwithoutpw — PDW — P0 (equation 43). 
Po 


P Pwithout DW 


(b) What is the sign of the effect in each case? 


E Solutions 


(a) Case (3), Z L M with and without DW => ^^ = 0. 


Cases (1), (2) and (4), Z || M without DW and 7 L M in the DW => 2e z 0. 
(b) The high resistivity state is obtained when I || M => py > ppw. 
Cases (1), (2) and (4) => ^? <0. 


6.3 Drift-Diffusion Equation for Spin Accumulation (u4) 
— chapters 3-5 


E Introduction/Reminder 
The distinct partial densities of current at the interface between materials of different 
types create a spin accumulation (or chemical potential imbalance) u, = —(u! — u!) 
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(83.2). Relaxation towards equilibrium conditions causes spins to diffuse near the 


interface, resulting in a diffusion-type equation for u,: V? u, — =z = 0 [equation (27)]. 
sf 


E Questions 
(a) Demonstrate equation (27), using the method proposed in $3.2. 


(b) Deduce how the law of total angular momentum conservation V -Js = 


e* (1-5?) 
m 


~(# — e) [equation (26)] takes the form V - Js; = u, [first term in 


5 sf 
equation (64), in §4.3]. 
(c) Show that equation (27) becomes V?u, — E = ( (V Ty! +S,V? T) — 0, with 
sf 
$,— -(s! — Bh if one considers the spin-dependent Seebeck effect (85.5), i.e. 
j(0 = MO yuQ) — gO gTOv T. 


E Solutions 


(a) Step 1: generalized Ohm's law 


vul? Z^ 


Step 2: charge and spin current 


J. —j +i! 
Js = -( - jl) 
Step 3: charge conservation 
V. Je =0 
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1 1 
Dig sco b ms gh 
=> Vu DELETE) 


Step 4: total angular momentum conservation (only spin-flip is considered in 


this calculation) 

ôn! dn! 
V Js = —e COT = so 
Ty Tg 

ôn! dnt 

— -(E-m 
Top Uy 
1 1 ôn! dnt 
Dy pes Dens 20 
=>V (ul —y') =e E + a) E =| (116) 


Step 5: obtaining the diffusion equation for spin accumulation 
Near equilibrium, we have 


" nD 
AR MW 
1 1 Niul Nut 
2(,1 — yt) = e | 
SU n-e( e) Tsp 


At equilibrium we have n! = N TQ). and equilibrium recovery after stochastic 
reversal imposes 
N! N! 
k= p= y 
Top Uy 


— Vy =0 (117) 
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(b) Combining equations (116) and (117) gives 


By afl 1\ fôn! dnt 
IX grec cse cl T gl qi. ql 
sf sf sf 
ls ,/(0 oV fôn! dn! 
ud [* a ola! Pal ql 
sf sf sf 


Using o?) = © (1 — (+ )$) the above expression transforms in 


Me Ae ôn! dnt 
S eH o 


Further introducting equation (26), one gets 


(c) Taking into account the spin-dependent Seebeck effect, the generalized Ohm's 
law writes: 


vil) =E Ot egtDy T 


al 


1 1 dS 
_ V2 — at E s 2 2 
=> —V?(u,) «(sv j V F) esa (VT)? — e8,V?T 


Based on point (a) we obtain: 


d$, 
> Vy, — = = (a (VT)? + SV? r) =0 
sf 
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A practical use of this type of thermoelectric effect and related equation is for 
example shown in A. Slachter et al., Nat. Phys. 6, 879 (2010). 


6.4 Spin Conductivity Mismatch (CPP-GMR) 
— chapter 3 


Bi Introduction/Reminder 

We consider a F/I/N/I/F multilayer ($3.3) consisting of two metallic ferromagnets 
(F) made of the same material, one metallic non-magnetic spacer (N) of thickness 
dy and spin diffusion length l; y, and two interfaces (I) (Inset of figure 67). 


AR/R' 


0.00 = c = 
1E-15 1E-14 1E-13 1LE-12 1E-11 1E-10 1E-9 1E-8 1E-7 1E-6 
* 2 
r, (Q.m) 


Fic. 67 — Calculated magnetoresistance (AR/ R?) vs. interfacial (I) spin resistance area 
product rj = rj(1— 52), for a F/I/N/I/F multilayer. Adapted with permission from the 
American Physical Society: A. Fert and H. Jaffrés, Phys. Rev. B 64, 184420 (2001). Copyright 
2001. 


For dy «& Us. ny» When an electric current is applied perpendicular to the sample 
plane, it can be shown that the difference AR between RP ) and RP is 


2(frr d yrs). 
2 
(rp ry) 4- 2 | (2) | n 


sf. N 
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and that 


2 | 2 L2 dx 
d (B — yY reri (Brr +y ri) aX rw 
N TN d 2(1 »)r, L2 PI 


(rr ri) + porn 


with 

RUP), the resistance of the stack when the magnetizations of the two F layers are 
antiparallel. 

R), the resistance of the stack when the magnetizations of the two F layers are 
parallel. 

ry, the spin resistance area product of the N layer (resistivity py times D, y). 

rp, the spin resistance area product of the F layers. 

rp, the spin resistance area product due to the interfaces, I (rp = Tinterface in $3.4). 
p, the spin polarization of the F layers. 

y, the spin polarization due to the interfaces, I 


E Questions 


(a) Justify why experiments on current perpendicular to plane magnetoresistance 
usually fulfill the following condition: dy < Uy N° 


(b) Express AR/R® (as a function of ff only) for the case with no interfacial spin 
resistance (ry = 0) and with only metals (ry ~ rr). Reminder: dy < [iy y. 


(c) Express AR/ RP, still for the case with no interfacial spin resistance (r, = 0) 


but with a semiconducting spacer (ry > rr). Consider 7X > Ti >L 

(d) From (c), what can you conclude on spin injection in a semiconductor? 

(e) Still for the case of a semiconducting spacer but with a non-zero interfacial spin 
resistance (r, Æ 0), we admit that it is possible to calculate the maximum value 
of magnetoresistance as AR/ R? = y?/(1 — y?). 
- From this expression, what is the result of adding a large interfacial spin 
resistance on spin-injection in a semiconductor? 
- What type of material that you know of gives a large spin-dependent 
resistance? 

(f) Calculate the three values of AR/R® from questions (b), (c) and (e), for: 
f —0.5, y 20.5, ry 24x10? O.m? (for the case of the semiconducting 
spacer), re = 4 x 10 P O.m?, I, y — 2 um, and dy = 20nm, 
and compare your findings with figure 67. 

(g) Illustrate the effect of the spin conductivity mismatch on the spin current. 
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B Solutions 


) Spin information is to be transmitted between the F layers, across the N layer. 
) ARRA EAR p) 
(P) WIN 
(c) AR/RU) — 4f? (cay, 
(d) 
) 


As is, spin injection in a semiconductor is not efficient because = uz * «I. 
- Adding interfacial spin resistances makes spin injection in a ciue 
possible. 
- Àn insulating tunnel barrier gives a large spin-dependent resistance. 
(f) Case (b), AR/ R? = 33% 
Case (c), AR/R™ = 1075 ~0 
Case (e), AR/ R? = 33% 
Agrees with data obtained in figure 67. 
(g) See figure 24. 


6.5 Intra-Band Scattering — Intrinsic Damping 
(xo) — chapters 2 and 4 


B Introduction/Reminder 

The damping parameter « illustrates how fast magnetization relaxes toward equi- 
librium, when initially brought out-of-equilibrium ($4.4). The intrinsic intra-band 
mechanism accounts for relaxation via electron scattering. 


N (eV 1) 


€ — Eg (eV) 


Fic. 68 — Calculated electron density of states (N) of Co,Feioox alloys for several 
Co-concentrations. All alloy compositions are aligned to a common Fermi energy, £p, at zero 
energy. Adapted with permission from Springer Nature: M. A.W. Schoen et al., Nat. Phys. 12, 
839 (2016). Copyright 2016. 
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E Questions 


(a) 


Give the typical equation describing magnetization dynamics, explain the dif- 
ferent terms and illustrate how they contribute. 

Which other mechanism that you know of also depends on electron scattering? 
Please detail how. 

From the calculations shown in figure 68, plot the dependence of intrinsic 
damping (o9) on the Co-concentration (x) in the Co,Fejg9_x ferromagnetic alloy. 
Please explain your plot. 

Which composition of the CoxFe100-x ferromagnetic alloy would you choose if 
you want to make a device with minimal fluctuations of magnetization after 
reversal, and why? 


B Solutions 


(a) 
(b) 


(c) 


See equation (65) and corresponding text, $4.4. 

For example, the electronic mean free path is governed by the Fermi golden rule 
and therefore by the DOS at the Fermi level (see 82.4). 

The intrinsic intra-band mechanism accounts for relaxation via electron scat- 
tering, meaning that, oj is proportional to the DOS at e = ep. From figure 68, 
we get: 


hasian tive 
0.003 " I" 
i r "omg 
a 
| = 
a F| 1.2 T 
0.002 - 4 — 
Si Kcu 
"x 309 < 
| 
0.001 F " 
—— Theory | 0.6 i 
Mankovsky et al. 
0.000 4 " 1 n 1 " i " 1 " 1i 
0.0 0.2 04 0.6 0.8 1.0 


Co — concentration 


Adapted with permission from Springer Nature: M. A.W. Schoen et al., Nat. Phys. 12, 839 (2016). 
Copyright 2016. 


(d) Large damping is suited for minimal fluctuations of magnetization, as it ensures 


fast relaxation. Based on (c), the largest intrinsic damping is obtained for the 
largest DOS at € = £r, i.e. near the CoggFe49 composition, from figure 68. 
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6.6 Spin Pumping (SP) and Inverse Spin Hall Effect 
(ISHE) —- chapters 4 and 5 


E Introduction/Reminder 

The spin pumping effect [Y. Tserkovnyak et al., Rev. Mod. Phys. 77, 1375 (2005)] 
refers to the ability of a magnetic material to generate a spin current J? when 
brought out-of-equilibrium (84.4). The technique usually involves inducing reso- 
nance in a ferromagnetic (F) spin injector — e.g., a NiFe layer — which is adjacent to a 
non-magnetic material (N) known as the spin sink — e.g. a Pt layer (figure 69). Spin 
pumping and spin transfer torque are reciprocal effect. An intuitive picture consists 
in comparing spin transfer torque to a water flow (spin current) moving the blades of 
a watermill (magnetization) and spin pumping to moving blades (magnetization) 
creating a water flow (spin current) (figure 41). 


Fic. 69 — (a) Illustration of the spin pumping effect due to sustained out-of-equilibrium 
magnetization dynamics. (b) Illustration of spin current diffusion across the non-magnetic layer 
(here, Pt). (c) Equivalent circuit considering spin current to charge current conversion due to 
the inverse spin Hall effect in the non-magnetic layer. Also figure 42. Adapted with permission 
from AIP Publishing: K. Ando et al., J. Appl. Phys. 109, 103913 (2011). Copyright 2011. 


E Questions 


(a) Because the system is out-of-equilibrium, spins accumulate at the F/N interface 
and diffuse across the N layer [figure 69b]. The time-varying spin density 
(nonequilibrium chemical-potential imbalance) can be written as follows: 
um ue. In this part, you will calculate the y-dependence of the spin current 


in the N layer: J,(y). 


Exercises and Solutions 115 


(al) Write the relation between J, and p,. 
This will be equation (al). 


(a2) Write the transport equation that regulates u,, and show that 


du, 1 
bou (a2) 
2 #2 s 
d qm 
with U = ELS and lsp = ,/ Dv, where D is the diffusion constant, 


and t*, is the average spin-flip scattering rate in the N layer. Here, Č, = 
sf g 8 y sf 


* # # 
y and I7, — I5 y. 
(a3) Is it realistic to consider that Lu ~ig? 


(a4) Write the boundary conditions at y = Oand y = dy. 


(a5) Use the result of question (a4) and the fact that the solution of equation 
(a2) takes the following form: u (y) = Aer + Be’, to explicit ,(y) 
us. y, dn, Usps e, and py. 


(a6) Use the result of question (a5) and equation (al) to explicit J,(y) vs. y, 
dy, Ip, and J}. 
This will be equation (a6). 
Due to the inverse spin Hall effect (ISHE) (chapter 5), the spin current is 
converted in a transverse charge current Je along x. The N layer then becomes a 


‘source’ of charge current. The spin current to charge current conversion is 
expressed as J.(y)= Osu J,(y), where Üsug is called the spin Hall angle. 


(b1) Use equation (a6) and the indications above to explicit the average charge 
current density Je = (Je(y)) =a 0" J (y)dy vs. y, dn, Uy, and d 


(b2) Plot the charge current I vs. dy. Comment the trend for dy < Uy and for 
dy > Uy. 


(b3) The equivalent circuit of the F/N bilayer is illustrated in figure 69c. 
Explicit the electromotive force (voltage Vigug) due to the inverse spin 
Hall effect in the N layer induced by spin pumping in a F/N bilayer. 
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(b4) Calculate the value of the ‘spin-charge conversion efficiency’: Osnelsp, 
considering that dyoy >> drop, and for Visae = 4 uV, dy = 10 nm, 5, = 
3 nm, ay = 4 x 109 Sm *, w = 0.5 mm, and JU = 4.8 x 109 Acm *, 


(b5) A contribution to the spin Hall angle is due to extrinsic sd scattering on 
defects. Given this specific contribution — intrinsic contributions are not 
considered here — how would you proceed to increase the value of the spin 
Hall angle in a material? Will the options you propose improve the 
‘spin-charge conversion efficiency’? 


It is possible to show that the angular dependence of the inverse spin Hall 
voltage is 


2 
wOsnel,tanh ( dn /(2 I) egy" (m hg) 


dyon + drop 8na c 


Visug = sin(0y)T, (c) 


i M s||Sim" (0r) + " (up MsySim (0). +40? 
(uy MsySin (8). +40? 

Üsug the spin-Hall angle, h,; the rf excitation field used to reach F 

resonance, and o the resonance angular frequency. Ow is defined in figure 70. 


with T = 2% 


, w the width of the layers, 


T is a parameter accounting for the trajectory of the magnetization pre- 
cession in the (zy) plane. It can be viewed as the elliptic- to 
circular-trajectory ratio. gl! is the real part of the spin mixing conductance 
per unit area per quantum conductance per spin channel, accounting for the 
ability of the F/N interface and N layer to absorb the spin component along 
M x (m x M). 


(cl) Show that T is a dimensionless parameter. 


(c2) What is the trajectory of the magnetization precession for 03; = 0? Why is 
it so? And should it always be the case? 


(c3) The angular(0g)-dependence of the magnetization's tilt 0j deduced from 
experimental data for the case of a 8 nm-thick NiFe film is given in 
figure 70. Oy is the angle between the applied magnetic field H and the 
normal to the sample's surface y. Comment this behaviour: what governs 
it? Use equation (c) and figure 70, to hand-sketch the 04-dependence of 
Visug. Comment the symmetry of Vijggg with H. 
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270 
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Oy C) 


Fic. 70 — Typical angular(0;;)-dependence of the magnetization's tilt 0y. From O. Gladii 
et al., Phys. Rev. B 100, 174409 (2019). 


E Solutions 


(al) The relation between J,and u, is 


(al) 


Note: take o! = o! = oy/2 


(a2) The time-variation of the spin density ji, results in spin diffusion, which is 
balanced by spin scattering. The transport equation becomes 


dü, jd, i, 


di dy? 7 , 
d? 
Lom = a m Pa 
dy Tsf 
du 1 
s 1,=0 a2 
d P Hs (a2) 


* 


i [ : . ; 
with ur = jae and Uy = Di, where D is the diffusion constant and Uy 


+ det 


is the average spin-flip scattering rate. 
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Ur 
a/1+ iot, 


frequency for a ferromagnet is f ~ 10 GHz (with uo9H ~ 0.1 T), and the typical 
value of spin-flip scattering time is týp~10fs, so ov ~ 10 x10 74 <1 


(a3) We have i = with œt}, = 2nf1;,. The typical value of the resonance 


=> Dr ~ D. Note that this condition is no longer true in the THz regime, for 
example for antiferromagnetic resonance. In the following, we consider that 
men. 


(a4) At the boundaries, we have: 
for y = 0, J; =“ (0) = J9, 


~ 2epy dy 


for y = dy, Js = 54-35 (dy) = 0. 


~ 2epy dy 


(a5) The solution of equation (a2) takes the form 


u(y) = Aer + Be "II, 


for y 20, 41-55 = J? => A— B = 2epy JT, 


> 2epy dy 
for y= dy, TR de =0 => B= Ae 
i. P s " 
u(y) = 2epy JO iz RENI. (e + ML evt) 
Holt) = 2p! m 1 T (efi + eit dtt) 


Hal) = Repu Ti p " 3 
N/ bef 


(a6) Using the result of question (a5) and equation (al), we obtain 
E sinh ((y — dx)/;) 
sinh (diy /t) 


Is(y) = (a6) 
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(b1) 1 pay 1 fd j sinh ((y = dy)/T,) 
Je= Ug) f wdd- f Osu d 
N Jo N Jo sinh(du/1.;) 


1 1 dy 
eS Osu J? 7 f sinh ((y — dy)/t) dy 
N sinh (0) Jo 


Je = ce L sosh(4u/1) 
e — SHE" s 


* dw sinh (dx/T;) 


AEN 2sinb? (ay / (21,)) +1 
Osun J? 
dy 2sinh (dy /(20,) )cosh (dy/ (21;)) 


(b2) I= Jeldy = Osypllt;tanh (#) J? 


dy 
ar, 


Ospe lly 7 


dy 
als 


For dy < D. I ~ Osypl ou d Spins are still polarized and get converted effi- 
ciently in the N layer. The thicker the layer, the more spins are converted. 
For dy > D I = Üsug Ue Js . The signal levels out. The part of the N layer in 
contact with the F layer (dy < 2/;,) converts spins. Above dy ~ 2l; spin-charge 
conversion becomes inefficient because the spins are depolarized. 


120 


(b2) 


(b4) 


(cl) 
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pru pyw 
V _ Ry Ry pate ly ig 
ISHE — R R — pru a ppu Y e&N 
F+ 4UN ldr ldy 
wdy 
VIsHE 


E dxN6y + dpop 


= Osme — tanh 0 
Visug SHE ae an | ) : 


The Hall angle is Osme ^ 696 and the ‘spin-charge conversion efficiency’ is 
Osun ls ~ 0.18 nm |J. C. Rojas-Sanchez et aL, Phys. Rev. Lett. 112, 106602 
(2014)]. 


Increasing the number of scattering centers will increase the extrinsic contri- 
bution to Ogg, because for this extrinsic contribution sgg ox i and increasing 
sf 


the number of scatterers will reduce Up This solution is however not appro- 
priate to increase the ‘spin-charge conversion efficiency’ because Osne bsy will 
remain unchanged. 

Increasing the scattering efficiency by using heavier materials, with a large 
atomic number Z, is another solution to increase Üggp. Because spin-orbit 
interactions roughly vary as Z*, this solution will result in a large increase in 
*spin-charge conversion efficiency'. Note however that other effects contribute 
to Üggg and the picture presented above is certainly more complicated. In 
particular, intrinsic contributions do matter, see chapter 5. 


+2 +2 2 , 
= uo M s|y|sin (Om) + "TEES (0)) +40? 
1220 


a =, 
(1M sySin (0r) +40 

y is in unit of Hz.T *; Lio M s is in T; and c is in Hz, so we can conclude that I is 
a dimensionless parameter. 


For 0y = 0, Fl — 1. By definition, the magnetization trajectory in the 
(x,y) plane is thus circular. In this case, the magnetization rotates about the 
out-of-plane direction. If there is no in-plane anisotropy, it is expected that 
the trajectory is circular. It would not be the case in the presence of an in-plane 
anisotropy. 


The magnetization is mostly in-plane, whatever the angle of the applied field. 
'This is because of the demagnetizing field. 
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Visue (norm.) 


Oy C) 


From O. Gladii et al., Phys. Rev. B 100, 174409 (2019). 


Visur is an odd function of H. Angular-dependent symmetries are useful to 
disentangle the inverse spin Hall effect and several others that may occur 
concurrently [M. Harder et al, Phys. Rep. 661, 1 (2016)]. When two effects 
share the same angular-dependent symmetries, frequency-, temperature-, 
stacking order- etc. dependences are used to unravel the contributions. 


6.7 Spin Pumping (SP) — Additional Damping (ap) 
— chapter 4 


E Introduction/Reminder 

'This exercise is related to the previous one. Its purpose is to demonstrate that spin 
pumping ($4.4) generates a damping term and to calculate this term. Magnetization 
dynamics (figure 69a) is described by the LLG equation. In the frame of spin 
pumping, one takes into account a term accounting for a loss of spin current due to 
the transfer of angular momentum at the F/N interface, at y = 0. The LLG equation 
takes the form [R. Igushi et al., J. Phys. Soc. Jap. 86, 011003 (2017)]. 


aM ao dM 
x 


| | SP 
di 1M x (ugH.er + uly) 4 M: di Tsrr 


y 


'The torque Tou can be obtained based on the conservation law of total angular 


y 


d h 
momentum, which local formulation is Aa |Y: (54) = —|)|T [equation (75) 
€ 


di — 
and corresponding text|. From figure 69a, we see that, when incoming electrons 
provide m to the magnetization M, this latter orients towards m, and M > 0, in 
agreement with the ideas discussed in chapter 4. 
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'The LLG equation then becomes 


dM 
di x yM x (HoH est ale Ligh) I 


ao dM I» h 0 
M 
Me ae donee! 


with J? the spin current at the interface [Y. Tserkovnyak et al., Phys. Rev. B 66, 
224403 (2002)] 


€ dM € tl dM 


o back __ TL 
J; = dee ^ = 2n M2, Ip eM x dt 2nMg Tit dt 


We recall that the definitions and uses of J; and y may vary, resulting in a factor i 


difference and +1 or cj, respectively depending on the articles and books con- 
sulted. In this textbook, the spin current density is J, = -( — j'), expressed in 


y 


A.m ? and the gyromagnetic ratio y<0 is taken as — 


E Questions 


(a) Inject J? in the LLG equation and show that the total damping takes the form 
& = a + Xp. Give the expression of «p. Note: we will make the approximation 
that the imaginary part of the effective spin mixing conductance hee averages 
out (FL term page 58). 

(b) For this question, we consider that: 

- the initial spin current pumped writes as follows: 


pump |. 
J =e 


- the backflow spin current at the interface, at y = 0, relates to the transfer of 
angular momentum from spin accumulation in the non-magnetic layer to the 
ferromagnet magnetization as 


0 0 
bak ___Ê M Asy € fl Hs 
Js aan or M * CE 3 p ter 


Note: we will also make the approximation that the imaginary part of the spin 


mixing conductance git averages out. 


(b1) Express gl! as a function of g!! when J>** = 0. 
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(b2) Based on the drift-diffusion equation for spin accumulation, give the expression 


of u?, then deduce the expression of J'*** 


| 1 20 ey b ~ + * 
lu € [- T how tanh (&)] , where Us = MIT 


and finally demonstrate that 


(b3) Comment this last formula when i; > dy. 


E Solutions 


(a) Based on the introduction/reminder part, and neglecting 255 the LLG 
equations becomes 


y A dM 
4 git )u x 


1 
TL = —|[y|M x (Het + uoh, 3 25 
ly| x (ut ff + Ho s) er Mg c Msdp An reff dt 


Due to spin pumping (sinking), the damping of the ferromagnetic layer 
therefore writes 
X= Xo + Ay 


with 
[y] A Ju 


Xp = Msdpr 4n Ir eft 


b1) If one neglects J"***, then J? = J?""P and 
g sS S S 


Tb rir 
Ir eff P gn 


(b2) The expression of u? can be obtained by solving the drift-diffusion equation for 
spin accumulation in the non-magnetic layer, for y — 0 
u, 1 


— 4, =0 
dy? m a 


At the boundaries, we have the following conditions: 


for y = 0, z4- 75 (0) = J?, and for y = dy, oa Te (dy) =0. 


The solution of the drift-diffusion equation takes the form. 


u(y) = Aer + Be "II, 
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=(0) = J} => A — B= 2epy 0L 


for y = dug. di (dy) 20 => B= Ae’ 


for y — 0, bns E 


We E " " 
u(y) = 2epy J? pale ds (e! y p e Mg ean v) 
— ECON sf 


e Only laf JE JU, 2d /t 
ij! y/ ts —y/ ley N/E, 
us) = 2epy 27 en anlly 1 eeanlly (e can - ) 


cosh((y — dy) /t) 
sinh (dx/L;) 


us y) E —2epy JS isy 


Wf = —2epy PE coth( dw/tp) 


2el 
p=- 5 coth( dy IE) 
zm 2M. 


Note: see also answers to question (a) — al to a5 — in the previous exercise. 


Neglecting the imaginary part, we have 


0 
back nu Hs 
J; "xai mx (su 


(at) 


yo, based on the expression of J}; Jem and Jor and 


e 
JE = " 7M3 gM x UE 


From J? = eum e 
S S 
neglecting the imaginary parts, one obtains 


+ M 
nM Pre dt mM O dp ^ 2g M Ir 
(m X Ge Leoth(dy/t,)J j 
M 
ERE gi EE TIE f gM 


e TA 
“eM x Mx 
2n MI |n "qp ag MZ) dt anM 
2e[, dM 
sf * etl 
th(d I) a ee 
joy PANES (imita i JJ 
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(b3) 


6.8 


dM dM 
| 
gh 4M x rum NM x s x g!M 
2el;, d aM 
x | M x —* coth A E 5 gl! gM x — 
2el* d e 
sf N 
Ir eff — gi gi! Z0 Om a a 
er, : 

grat = ott — gi "i cot ( du) ger 


2er 
gl stenb (dy /U) = gH tanb(dy/T;) = gu a ie 


ae) 4 gt 20 er] ot x 
‘eff tanh (42) +g, p =g} tanh (dy/T) 


gl tanh(ay/t;) 


Ir eff = 2 pk 
d 2e, 
tanh (2) + gil i 


e 
35 


D hon 


By comparing the results of questions (b1) and (b2), we conclude that taking 
into account the backflow spin current renormalizes the spin mixing 
conductance. 

We also see an expression that looks like resistors in series. More information 
on the spin-circuit representation of spin pumping can be found in K. Roy 
et al., Phys. Rev. Appl. 8, 011001 (2017). 


When Uy > dy, ee = 0. Indeed, no spin has been absorbed and all spins have 


been reflected: J?""P = J?***, and thus J? = 0. In other words, no net transfer 
of angular momentum took place. 


Spin Hall Magnetoresistance (SMR) — chapters 2-4 


Bi Introduction/Reminder 

'This exercise is based on a spintronic effect called the spin Hall magnetoresistance 
(SMR), published in H. Nakayama et al., Phys. Rev. Lett. 110, 206601 (2013). It 
illustrates knowledge gained in chapters 2-4. A bilayer is considered (figure 71). It is 
made of a metallic heavy *non-magnet' (HM), like Pt, and an insulating magnetic 
layer, like YIG (F) (YIG = Yttrium iron garnet, Y 3Fe5O,). 
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(a) (b) 


9, 


I J, 

* "y N 7 as 

s) d — Sy 
» sg AS 
o E. f i N 


J 
LJ 


Fic. 71 — Illustration of the spin Hall magnetoresistance effect. From S.-Y. Huang et al., Sci. 
Rep. 8, 108 (2018) - CC BY Licence. 


A charge current Je is passed along the X direction in the HM. The spin Hall effect 
(SHE) converts this charge current in a transverse spin current Js. The amplitude of 
the spin current is controlled by the orientation of the F magnetization M — M sM, 
due to transfer of angular momentum via spin transfer torque. We note a, the angle 
between M and y. 

In the illustrations of figure 71, Js = J” = J” (z)2, meaning that the spin current is 
made of spins polarized along y only (first part of the equation) and that those spins 
flow along Z only (second part), in a way that is z-dependent. 

For « = n/2, figure 71a, when M = y, the polarization of the spin current in the HM 
and the magnetization of the F are parallel. As a result, there is no transfer of 
angular momentum. Spin accumulation in the HM is p, = 1, (z)y, meaning that the 
spin accumulation vector has a component along y only. It creates a backflow spin 
current, which compensates for the triggering spin current. The resistivity of the HM 
is low. 

For « = 0, figure 71b, when M- X, the polarization of the spin current in the HM 
and the F magnetization are perpendicular. As a result, there is transfer of angular 
momentum and the spin current is partially absorbed by F. Spin accumulation in the 
HM still writes p, = “,,(z)y. It is smaller compared to the case « = n/2. The 
resistivity of the HM is high. 

The purpose of this exercise is to calculate the change in resistivity with o, i.e. with 
the orientation of M, which can be changed, in practice, by use of an external 
magnetic field H. 

'The general Ohm's law in HM is 


Je 1 ÜsugX x  ÜsugyX OsHEZX E 
Jg | o | osx 1 0 0 Vy.,/ (20) 
J! | ^F | suy x 0 1 0 Vu, / (20) 


J; ÜsugZ x 0 0 1 Vu,./(2e) 
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Here, E = E,X is the applied electric field, — e is the electron charge, oy is the HM 
conductivity, and Üggg is the spin Hall angle of the HM layer. 


E Questions 


(a) 


(b) 


We first start by treating the case x = nx/2 [until question (f)]. Give the 


expression of e as a function of [s ) 

S 
Use the diffusion equation for spin accumulation to find the expression of 
Hs, (z). For this part, consider (i) a semi-infinite layer and (ii) the following 
boundary condition at z = 0, i.e., at the interface between HM and F 


€ asl; YA 
J0) = - 5a M x (Mx =) 19 


What is the physical meaning of the above mentioned boundary condition? 
Using questions (a) and (b), give the expression of J. = J.(z)X 

Write Je(z) in 2 form J.(z)= J°+J°f(z) and use the approximation 
Je= JO i Ty o f(z)dz, to find the ‘average’ charge current. dy is the 
thickness of HM. 

Set Je = o,,E,, and expand to leading order in up to find the expression of 
Prz = Pxv,x—n/2- 

For «=0, go over questions (a)-(f) again and give the expression of 
Prz = Prva=0° 

The expression of the spin Hall magnetoresistance (SMR) is defined as 


Uy mG 
MP xx = Prra=0 T Prr=n/2 — Po T Prz Verify that AP on = m Pn Osun a dy ae 


The actual full angular dependence of the SMR signal is given as follows (see 
figure 72 for the definition of the angles): 


= 2 
Poe, smn (X) = por Ap... COS x 


«sun (P) = Po + Ap,,sin? ff 
Por sMR(Y) = po 


Explain why p,, sun (7?) = po is independent of y. 
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Fic. 72 — Illustration of the angles considered. Adapted with permission from the American 
Physical Society: H. Nakayama et al., Phys. Rev. Lett. 110, 206601 (2013). Copyright 2013. 


(j) We now consider that the HM layer is slightly magnetized due to proximity 
effects, meaning that it now consists of a metallic ferromagnet with magneti- 
zation Myy = M eM gw. 'This magnetization does not influence the SMR 


signal. However, it creates an additional signal due to anisotropic magne- 
toresistance (AMR). 


(j1) Explain what AMR is. 


(j2) Write how AMR depends on the orientation of the external magnetic field 
H, a, f, y 


(j3) Which angular dependence would you use to discriminate between SMR 
and AMR, and why? 


B Solutions 


(a) For « = 1/2, we have M =9, and Hy = ns (y: 


OSHE Obey 
E; 0 E 
(5)-» Tae " (3) 
Ji 0 eae Hey Z 


2 
y LE 0 
2 x2 sy 
dz ly 


For a semi-infinite layer: 
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Boundary condition: 


T € lar m mi A 
J(0) = — 5,9 M x Mx% y 


Here, with M = y, and p, = Hsy(2)9: 


Ji(0) =0 
Hence: 
on OH, 
—onOsup E, + —— (0) — 0 
o NÜspgg Es + De d (0) 
—ÜsggE, — — = 0 
SHE Ëz elk, 


A = —2el, Osun E; 


us (2) = —2el* Osun Eze "I 


(c) 


At z — 0, the spin current can be partially absorbed by the ferromagnet due to 


transfer of angular momentum via spin transfer torque. This is especially true 
for « = 0. 


(d) 


Je(2) = on Er + ypo Ere ^! v 


Jel) = J? + agp Le 


(e) Averaging over the film thickness: 


1 9g k 
Je = Jo | Gas Joes dz 
N JO 


* 


Č. 
Je = J? + d. n 
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Setting Je = os Es: 


; F 2 2 
Expanding to leading order in Ogg: 


Pax = Proxa=n/2 = PN (1 =< bs 


For x = 0, we have M = €, and p, = Moy (2)¥: 


E OSHE Os, 0 
i3 2 6 ^ 2e Oz 
s —Üsgg Er + ——— 
i SHE^'T 96 Oz 


Spin accumulation diffusion equation: 


2 

Duy 1 ines 
2 x2 sy 

dz Uy 


For a semi-infinite layer: 


Boundary condition 


1) - - | ott x (M xt) 


Here, with M =X, and Hy = ns): 


J4(0) = —eg!! u, (0)/h 


Hence: 


ON Ol, 
2e Oz 


—on Osun Er 4 


N) &) 


(0) = —eg}! u, (0)/h 


Exercises and Solutions 


onA 


—— = —eglA/h 
2e, °9r / 


—on Og E, — 


2et; JOSHE Es 
aek Na 


As a result: 


with 


e y 
Je(z) = oyvE,+ Osn E; 


14 2e%,g} /(hon) 


eis 


Jel2) = J? Out 
1+ Dg. Go/on 


with Go = 2@/h, the quantum conductance. 
Averaging over the film thickness: 


04 5 e s 
Jem Jat zl : dz 
: sind " + l a Go/oNn 


č 1 
Je= je ge 0 
(dy SE Egi Glan ° 


Setting Je = o4 Es: 


* 


Or = on{ 1+ ly 92 l 
A dy BEL +E gt! Go/on 


Expanding to leading order in 0gy,: 


Uy 2 1 


TT T TL, A=! = 1 0 
Por = Prr a=0 Zz du SU TE gl Giov 
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(h) Setting MP xx = Pura—0 T Parra=r/2 — Po — Pr/2' 


* 


D. 1 Le 
Ap, = py | 1— 2-02 (1 sp ) 
Pa; Z dy SHE] n Egl 1 G m DN dy SHE 


2 Uy ight Go 
SHE ; 
dy oy + Dg) Go 


(i) Whatever y, M is perpendicular to p. 
(j1) See §2.6. 


(j2) Here, the angular dependence of AMR can be written as follows: 


Prz AMR (2) —p,-— Ap , cos^a 
Pas, AMR(B) = Bi 
Pas AMR() = p, — Ap sin 


with Ap, — p, — P] 


(j3) y can be used to discriminate between AMR and SMR. (Reminder: 
Pox smR(?) = py) 


Additional note: Here, we used a semi-infinite layer to ease the calculations. 
The exercise can also be done if one considers a finite size layer of thickness dy. 
Basically, this change will modify the boundary condition in question (b), the upper 
bound of the integral in question (e), and the final results — one obtains: 


" I Tl 2 dy 
Ap, = pol Iy i9, Cotanh a% 02 zu P S 
li SHE dw on + g! Gocoth 42 SHE dy 2b, 
sf 
o E, yg, Gotanh? se 
~ po Osun > " TE 
dN oy + tg! Gocoth 2 


The physics at stake and the main conclusions remain unchanged, in particular the 
Goan dependence. For a full description, see H. Nakayama et al., Phys. Rev. Lett. 110, 
206601 (2013), and Y.-T.Chen et al., Phys. Rev. B 87, 144411 (2013). 
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6.9 Harmonic Analysis of the Anomalous Hall Voltage 


(v ve) — chapters 4 and 5 


E Introduction/Reminder 
This exercise illustrates the harmonic Hall voltage measurement technique, which 
has become a standard technique to determine the symmetry and magnitude of the 
effective fields Hp, and Hp, originating the spin-orbit torques Tfr, and Tp, (84.7). 
It is based on U. H. Pi et al., Appl. Phys. Lett. 97, 162507 (2010). 
A Hall-bar structure is patterned into a heavy metal non-magnet (N)/ferromagnet 
(F) bilayer and a periodic charge current I, = Isin(wt) is driven along the longi- 
tudinal arm (x) of the Hall-bar, therefore generating a periodic spin-orbit torque 
(figure 73). In this exercise, to simplify calculations, we consider that only a periodic 
SOT term (§4.7) originating from a periodic effective field along y, of the form 
Hsor,, = Hsorsin(cot) is generated. Under the application of a static transverse 
external field H,, the tilting of the magnetization out of the easy axis (here, Z) is then 
modulated by Hsor,, induced by I. The anomalous Hall resistance Rry, which 
relates to the z-component of magnetization (chapter 5) thus takes the form: 

AR», (H) 


Roy (H) = Fig Hy) + VUES Hsorsin(ot) 


Isin(wt) 


Fic. 73 — Illustration of the device used for the harmonic Hall voltage measurement 
technique. 


E Questions 
(a) Show that the Hall voltage takes the form 
V(H)— Vi, + V$ sin(@t) + Vie cos(2at) 


and explicit V? , V®, and V" as a function of Fa Hy); Sa Hsor, and I. 


zy? © vy? 
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(b) Express Hgor as a function of V9? 


©, or its derivative, and V5». 

Experimental data of V7 and m are given in figure 74. 
(c) Explain why Led is non-zero for the //Pt(3)/Co(0.6)/AlOx(1.8 nm) trilayer 
and zero for the //Pt(3)/Co(0.6)/Pt(3 nm) trilayer. 
(d) From figure 74 and point (b), estimate and plot Hgor vs. Je. Comment the data. 
) Compare the Oersted field Hoe created by I to Hsor and comment. 


J, (10 A/onf) 
80} —2— 2.1 
> ——42 > 
E! € 4-64 3 
a apl ——85 3a 
Sk - —— 106 = 
20} 
0r 
00 02 04 06 08 10 00 02 04 06 08 10 
HoH, (T) HoHy(T) 


Fic. 74 — First Vj, and second uv harmonic transverse voltages measured as a function of 
the external transverse field H, for a Ta(5)//Pt(3)/Co(0.6)/AlOx(1.8 nm) trilayer (main 
figures) and a Ta(5)//Pt(3)/Co(0.6)/Pt(3 nm) trilayer (inset). Adapted with permission 
from the American Physical Society: U. H. Pi et al., Appl. Phys. Lett. 97, 162507 (2010). 
Copyright 2010. 


B Solutions 


(a) V4(H) = Roy(H) In 


_ day (H) 


1H Iorsintot) Isin(ct) 


with 
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V? = IRA(H) 


dRa(H) 


20 — | Hao /2 
Va jp or / 
20 
b Hsor = 2——3— 
(b) SOT "dVo/dH 


(c) V2 is non-zero for the / /Pt(3)/Co(0.6) / AIOx(1.8 nm) trilayer and zero for the 
//Pt(3)/Co(0.6)/Pt(3 nm) trilayer because the SOT considered and its related 
effective field requires symmetry breaking, here structural inversion asymmetry. 

(d) Estimating and plotting Hgor vs. Je gives something similar to the actual data 
plotted below. We note that Hsgor is linear with Je, as expected from its origin. 
Above a threshold current density J,~8 x 109 A.cm ?, Joule heating makes 
Aggy depart from its linear relation with Je. From the slope of Hsor vs. Je, we 
obtain ji Hsor[mT] ^ 3J. [109 A.cm ?]. 

We note that V2 is not used for the determination of Hsor because the dc 
voltage is usually contaminated by spurious contributions for example due to 
misalignments of the different branches of the Hall bar. 


M~ HoHy- 
— | -—oe9T 
= —— 08T 
E W LoT 
S 20 
S 20 
X 
© L 
S10 

Or 


024 67 9 t 
Je (10° A.cm?) 


Adapted with permission from the American Physical Society. 
U. H. Pi et al., Appl. Phys. Lett. 97, 162507 (2010). Copyright 2010. 


(e) From pp Ho. * 2 * width uy Je * surface, we estimate uoHoeļmT] ~ 0.5Je 
[109 A.cm ?] . The contribution from the Oersted field is thus smaller compared 
to the findings. 
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Note: in this simplified exercise, we considered a magnetization pointing 
out-of-plane and a damping-like torque only, based on U. H. Pi et al., Appl. Phys. 
Lett. 97, 162507 (2010). A general scheme to measure the amplitude and the 
direction of the effective fields originating SOTs, when the magnetization direction is 
random and when field- and damping-like components are at stake, is provided in 
K. Garello et al., Nat. Nanotech. 8, 587 (2013). Symmetry arguments also allow 
unravelling the different contributions to the transverse voltage, including, for 
example, the one due to anisotropic magnetoresistance (see exercise 6.1). 


6.10 Intrinsic Anomalous Hall and Nernst Effects 
(AHE, ANE) - chapter 5 


E Introduction/Reminder 

A prototypical example illustrating Berry curvature physics is the case of a 
spin-polarized 2D electron gas with Rashba spin-orbit coupling (§5.4). The corre- 
sponding Hamiltonian (figure 59, bottom right) is 


27,2 


h 
H = 


M +op(k x o) -Z — ^o; 


From this Hamiltonian, the aim of this exercise is to demonstrate the expressions of the 
eigenvalues e4, eigenvectors |+), and Berry curvature Q4 [equations (107)-(109)]. 
The intrinsic anomalous Hall conductivity (figure 60) is then obtained from 
[equation (100)] 


on =< [ 3505001609) 5 


and its complementary anomalous Nernst conductivity from [equation (112)] 


E Oc. (&r) 
Go Oe 


Ory = —2e 


E Questions 


(a) Show that the above Hamiltonian can take the form 


+b(k)-6, with b(k) = (—arky, ok, —A). 


Reminder of the Pauli matrices: or = t 3 oy = ( n and 
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(b) Show that the eigenvalues can be expressed as 


+(k) = eg(K) + |b(K)]. 


[sz] 


(c) With these eigenvalues, demonstrate that the eigenvectors of H are of the form 


dE 
Tr 
~~ 
I 
j= 
| 
= 
=> 
~~ 
+ 
> 
ES 
[aun 


A 
c ibo e 


where |J), = 5 pitt ( J and |y), = 25 e** (2) are the eigenvectors for free 


electrons. 
(d) With these eigenvectors, from equation (87) 2,(k) = Vy x (W,(k)|¢Velw,(k)), 
show that the Berry curvature takes the form 


2 
ARA 


(=F 
i 2 (002, k? + A’) 


3/2 


Use the fact that |+) = e*rlus), and apply equation (87) to |u.). 


e 


E Solutions 


(a) The Hamiltonian can be transformed as follows: 


Wk? x 
H = —— +arlk x o) -Z — Ao, 
2me 


2,2 
diet +ar(zxk)-o- ^o; 
2me 
Ve 
H= +b(k) -o 
2Me 


with b(k) = (—arky, pk, —A). 
(b) In this type of formulation, the eigenvalues can be expressed as 


+(k) = eo(k) = |b(K)] 


hk? 
e (k) = 2m V OE + A? 


e 
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More precisely, in the formulation of the Hamiltonian, b(k) represents a 
point on the Bloch sphere, which has the consequence that the solutions of 
the eigenvectors of the total Hamiltonian are in the (V4); y/,)) basis of the 
original wavefunctions in the absence of the b(k) - o perturbation. Thus, the 
Schródinger equation 


H|+) = e|) 


has eigenvectors of the form 


|) = Ax |i) + Bx Iv) 


with |y) = 25 ent ( J and |yj) = Z prt ( ; ) , the eigenvectors for free 
electrons. 

After insertion into the Schrödinger equation, and taking ks+ iky = 
k(cos $ + isin $) = ke'?, and thus k, — ik, = k(cos — isin $) = ke“, we 
get 


hèk? ; 
A-é itg ke? 
2m, Ax 
PII ej 0 
iapke'® Jo es x 
Me 


Zeroing the determinant of the square matrix on the left returns the energy 


eigenvalues 
n 
& (k) = ES 1/03, k2 + A? = colk) + |b(k)| 


We now detail how to find the coefficients A, and B+}, which govern the 
eigenvector |+) = A, |V4) + B+ |Wy) associated to the eigenvalue £4 
Substituting € in the above matrix, we get 


—A — |b| —iopke*? Ax zT 
ike?  A—|b| )\ By) ~ 


This leads to a system of 2 equations and 2 unknowns 


(—A — |b|)A, — ixgke'? B, =0 


ixgke 1? A, +(A— |b|)B, —0 
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However, the system is degenerate, so we need to choose one variable and 
find the other. Starting from zero, we would choose A, = 1, then find B4, 
then use A*, A, + B*. B, = 1 for normalization. Because the solution is 


given to us, we choose A, = ,/1— a and find B, from 


A+ |b| A+ |b| A 
D ixgkei ^ * ig keió Ib| 
ix ke 1o iake’? A 
B = A — 
P A-|| ^ — |p| lb] 
S A+ [b A ae A 
B c 2i eee 2 —2ip 1 P 
eee smi i) idi Qi" 
A 


The coefficients A- and B_, which govern the eigenvector |—) = 
A. |V4) + B_|y) associated to the eigenvalue &.. can be found by using the 
same methodology. 

From 


|) = Ax |i) Beli) 


with Yp — 45 s(a) and |y) =. se(1) and having determined 


A,, B}, x and B_, we get the expressions of the eigenvectors 


I 
= 


wr uU 


i "ii k)| i k)| 


(d) The expressions of the eigenvectors can be transformed into 


) = eT lus) 
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with 


dis A i can) tay taai) 


and $ = arctan(k,/k,) and k? = k? + ky. 
We now detail the steps necessary to calculate Q4 (k), which is the Berry 
curvature along z corresponding to the eigenvector |u.,). 


Take 
= fi kr, ky 
Jus) = É kinky 


We then have 


io Of 
Lok, Ok, Q: 
(u4|Velu+) = qu Of, Ef Oh, = 4) 
1 ðk, ? ðky 0 
0 


Thus, we get 
(k) = Q,() -2 


Q5 (k) = (Vx x (us |iVa]us.)) 2 


Pine eel (On ON. s 
aid) = sw ( X 23] ; 


aA 
Q? (k) = R 
H 2 (922 + a)?” 


Using the same methodology for Q? (k), we obtain 
oA 


(ak? + AP? 


Q (k) =F 
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Note 1: the eigenvalues, eigenvectors, and (zero) Berry curvature for a Rashba 
Hamiltonian (non spin-split) can be straightforwardly obtained from the above 
expressions by setting A — 0 (85.4, figure 59 top right). 


Note 2: to go further, the general expressions of eigenvalues, and Berry curvature 
can be found in D. Culcer et al, Phys. Rev. B 68, 045327 (2003) for two-level 
Kk 
Me 
b(k) = (—arf(k), arf(k),—A). In particular, it is possible to show that Q X am 


degenerate systems with a Hamiltonian of the generic form H = +b(k) -o with 


which highlights the expression of a ‘magnetic monopole’ in k-space, sitting at b = 0 
[D. Xiao et al., Rev. Mod. Phys. 82, 1959 (2010)]. 


Chapter 7 
Conclusion 


By now, the objectives of giving the reader the tools to go further and understand 
other related or more elaborate phenomena of spintronics should have been 
achieved. Of course, “Everybody calls ‘clear’ those ideas which have the same degree 
of confusion as his own" — M. Proust. The aim of this conclusion is to open up the 
field of vision by presenting in a non-exhaustive way some of the current topics to 
which the interested reader may wish to turn. To this end, the selected topics are 
briefly mentioned and accompanied by references consisting of comprehensive 
review articles that provide a complete overview of the topic. 

It is important to mention at this point that the experimental measurement of 
the effects presented in this book is extremely delicate as it involves the fabrication 
of nano-sized structures, the detection of signals that can be of a few nanovolts and 
that can be obtained in some cases by means of pulses of a few nanoseconds or even a 
few picoseconds, as well as the use of advanced electrical, optical, and thermal 
techniques, not to mention the use of large instruments such as synchrotrons. 

In general, spintronics and the disciplines that interface with it, such as micro- 
electronics, quantum computing, biology, and medicine, call for teams of researchers 
and engineers with a wide range of skills in, for example: 


Bi solid state physics 

Bl materials science and engineering 

Bl instrumentation 

B theory 

Bi analytical and numerical computation 
Bi signal processing 

Bl integrated circuit design. 


'This wealth of skills needed by the discipline and of interest to society promises a 
bright future for those who wish to join it, whether they approach the discipline from 
the point of view of fundamental science or from an applicative perspective, in 
physics, mathematics, biology, or chemistry. 
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Moreover, the great diversity in which magnetic materials appear in nature — 
from metals to semimetals to insulators, not to mention the great wealth of magnetic 
textures and topological structures described by group theory, provides a vast and 
fascinating playground. Among the materials that are currently the subject of many 
studies are, for example: 


Bl oxides with their unique versatility and multiple functionalities — M. Coll et al., 
Appl. Surf. Sci. 482, 1 (2019) 


BB 2D materials because of their unusual nanoscale-related characteristics — 
S. Z. Butler et al., ACS Nano 7, 2898 (2013), Q. H. Wang et al., ACS Nano 16, 
6960 (2022) 


E other topological materials and textures, including: 


© topological insulators exhibiting topologically protected spin-polarized surface 
states — B. A. Bernevig et al., Nature 603, 41 (2022) 

o Weyl semimetals displaying Dirac nodes of opposite chirality - N. Kumar 
et al., Chem. Rev. 121, 2780 (2021) 

O the newly discovered unconventional ‘alter’ magnets with anisotropic 
higher-order partial-wave forms of magnetically ordered phases — L. Smejkal 
et al., Nat. Rev. Mater. 7, 482 (2022) 

© topological magnetic textures like skyrmions and other non-collinear struc- 
tures, exhibiting locally topologically protected states as opposed to the sur- 
rounding magnetic structure — A. Soumyanarayanan et al., Nature 539, 509 
(2016), B. Góbel et al., Phys. Rep. 895, 1 (2021) 


Bi antiferromagnets for which the presence of two sub-lattices of magnetic 
moments gives them unique properties in real and reciprocal space — 
T. Jungwirth et al., Nat. Nano. 11, 231 (2016), V. Baltz et al., Rev. Mod. 
Phys. 90, 015005 (2018) 


Bi ferrimagnets, which combine at compensation some favorable properties of an 
antiferromagnet with the net spin polarization of a ferromagnet — S. K. Kim et al., 
Nat. Mater. 21, 24 (2022) 


Bi heusler alloys for their specific density of states and corresponding large spin 
polarization and low damping — T. Graf et al., Prog. Solid State Chem. 39, 1 
(2011), K. Elphick et al., Sci. Tech. Adv. Mater. 22, 235 (2021) 


Bl organic materials, which are more abundant and cheaper than other materials and 
could be used for molecular spintronics — A. Cornia et al., Nat. Mater. 16, 505 
(2017) 


Bl artificial spin structures, such as artificial spin ices, which give access to some 
non-conventional behaviors such as emergent magnetic monopoles and their 
dynamics — S. H. Skjerve et al., Nat. Rev. Phys. 2, 13 (2020) 


E heterostructures combining different spin orders or ways to propagate angular 
momentum, thus being a source of new physics and a promise of new function- 
alities that arise from proximity effects and competing interactions. Striking and 
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exciting examples of how several areas of condensed matter physics can look 
forward to a common future are provided by the combination of: 


o magnetic and superconducting order, the latter bringing quantum coherence 
and dissipationless transport — J. Linder et al., Nat. Phys. 11, 307 (2015) 

o strongly correlated antiferromagnetism and high-T, superconductivity, where 
the intermingling between spin transport and antiferromagnetic order proves 
to be decisive — F. Giustino et al., J. Phys. Mater. 3, 042006 (2020) 


The wealth of materials available, both crystalline and amorphous, allows the 
selection of specific materials and combinations of magnetic (e.g., anisotropy, 
exchange, DMI, damping) and transport (e.g., polarization, conductivity, Berry 
curvature) properties for the study of specific current research areas, whether for 
example: 


E symmetries and topology, towards quantum (band-topology) spintronics and 
non-dissipative spin-conservative electronics — C. K. Chiu et al., Rev. Mod. Phys. 
88, 035005 (2016) 


H angular momentum transport and transfer, where several carriers of angular 
momentum are involved — electron, phonon, and photon — and wherein the 
amplitude or the phase is used for information encoding and transfer. This 
research includes several fields: 


© spin currents/spin orbitronics — S. Maekawa et al. (eds), Spin Currents, 
Oxford University Press (2012) 

o orbital spintronics — D. Go et al., Europhys. Lett. 135, 37001 (2021) 

© spin mechatronics — M. Matsuo et al., J. Phys. Soc. Jpn 86, 011011 (2017), 
J. Puebla et al., Appl. Phys. Lett. 120, 220502 (2022) 

© opto spintronics — A. Kirilyuk et al., Rev. Mod. Phys. 82, 2731 (2010) 

o magnonics — A. V. Chumak et al., Nat. Phys. 11, 453 (2015), IEEE Trans. 
Mag. 58, 0800172 (2022) 


Bi spin caloritronics, which considers the interaction between heat and spin, for 
example, for thermal energy conversion - S. R. Boona et al., Energy Environ. Sci. 
7, 885 (2014), K.-i. Uchida et al., J. Phys. Soc. Jpn 90, 122001 (2021) 


Bl magneto-ionic and electrostatic gating, in which an excess/depletion of charges is 
used to modify magnetism — M. Nichterwitz, Appl. Phys. Lett. Mater. 9, 030903 
(2021) 


Bl ultimate time scales, at which specific ultrafast physical processes occur like heat 
to spin conversion through lattice excitation, with potential for T'Hz applications 
— J. Walowski et al., J. Appl. Phys. 120, 140901 (2016) 


Bl magnetic nanoparticles, which are highly-functionalizable objects for potential 
use in the health sector - K. Mahmoudi e£ al., Int. J. Hyperthermia 34, 1316 
(2018), C. Naud et al., Nanoscale Adv. 2, 3632 (2020) 
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Bi 3D spintronics, wherein complex magnetic configurations, like Bloch points, with 
specific properties become possible — A. Fernández-Pacheco, Nat. Commun. 8, 
15756 (2017) 


Bl and more mature fields at present, like: 


o MRAM/design for IT applications like memory, processors, data security — 
B. Dieny et al. (eds), Introduction to Magnetic Random-Access Memory, 
Wiley-IEEE Press (2016), V. Krizakova et al., J. Magn. Magn. Mater. 562, 
169692 (2022) 

© sensors for IT, automotive, space, and biological applications — B. Dieny et al., 
Nat. Elec. 3, 446 (2020) 

o oscillators for use as transceivers for telecommunication applications or 
neurons and axons for neuromorphic computing and related AI applications — 
J. Grollier, Nat. Elec. 3, 360 (2020) 
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1. Abbreviations 


AHE anomalous Hall effect 

AMR anisotropic magnetoresistance 
AP antiparallel 

CIP current in plane 

CPP current perpendicular-to-plane 
crit critical 

DL damping-like 

DOS density of states 

DW domain wall 

DWMR domain wall magnetoresistance 
eff effective 

e.g. ezempl gratia (for example) 
et al. et alii (and others) 

F ferromagnet 

FL field-like 

GMR giant magnetoresistance 

ibid. ibidem (in the same place) 

i.e. id est (that is) 

IEC interlayer exchange coupling 
ISHE inverse spin Hall effect 

LLG Landau-Lifshitz-Gilbert 

N non-magnet 

OHE ordinary Hall effect 

P parallel 

QHE quantum Hall effect 

QWS quantum well state 

REE Rashba-Edelstein effect 
RKKY Ruderman-Kittel-Kasuya-Yosida 
SAHE spontaneous anomalous Hall effect 
SHE spin Hall effect 

SML spin memory loss 

SO spin-orbit 

SOT spin-orbit torque 

SP spin pumping 

STT spin transfer torque 

THE topological Hall effect 


TMR tunnel magnetoresistance 
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2. Symbols and Units 


Bold is used for vectors, e.g. r= T, an overcaret is used for unit vectors, e.g. 
T =r/r, an overdot indicates a time-derivative, e.g. r = Or/Ot, a double bar is used 
E id J^ JA 
E S S S 
for second-order tensor, e.g. J, = [JO J(? JO) =| JZ J" J? 
Jo Xue. Je 
Many of the parameters used are spin-dependent. In the text, they are accom- 
panied by arrows as superscripts, where necessary. Additionally, some of the 
abbreviations detailed above are combined as subscripts or superscripts with the 
symbols, where necessary to specify that a given effect relates to a given phe- 
nomenon, which is then detailed in the text. 
Some symbols may play a dual role, in which case their meaning is obvious from 
the context. 


Symbol |Formula | Quantity Unit Value 
Physical constants 
c velocity of light ms! 2.998 x 108 
magnetic permeability " E 
Mo of free space TmA™ 4r X 1077 
kp Boltzmann constant Jk! 1.381 x 10°73 
e elementary charge C 1.602 x 10 7? 
me electron mass kg 9.109 x 10?! 
h Planck constant Js 6.626 x 10 
h , 
h Dx Planck constant /(27) Js 1.055 x 10 %4 
eh : 
Lp 2m, Bohr magneton An? 9.274 x 10 ?* 
2€ electrical conductance E = 
Go E Q™ (S) 7.748 X 107 
quantum 
2 
Rx Go von Klitzing constant Q 2.581 x 104 
n2 A T thermal conductance l 
9o 3h quantum WK? 9.464 x 10 ^ x T 
Poy Bu magnetic flux quantum T m? (Wb) 2.068 x 10 7? 
1 T! m? 14 
m. S S 3 zs ; x 
Ky d, Josephson constant (Wb m 4.836 X 10 
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Symbol} Formula Quantity Unit 


Arrows and relational operators 


(spin) magnetic moment m 
orientation? 


(spin) magnetic moment m 
orientation 


spin angular momentum 8 
orientation? 


spin angular momentum 8 
orientation? 


n spin mixing — moment mixing 
| parallel 
L perpendicular 


Operators, distribution functions and related quantities 


0 0 8 gradient 
v (2:5; , =) formula in cartesian coordinates 
fa) fa) ð divergence 
y. aota Tta . ‘ E i 
Ox Oy Oz formula in cartesian coordinates 
T transpose 
H Hamiltonian 
FECE 
EE TUE 
ip Or 
o dip E) Pauli operator 
= 1 0 
(O40 -1 
S h A spin angular momentum operator 
2 
ĝa aað orbital angular momentum 
ih| — — — @— t 
L sind 0$ n 00 opera i) i , 
formula in spherical coordinates 
gut spin angular momentum ladder 
S operator 
RE orbital angular momentum ladder 
L operator 
nh 
—i— V : ; 
v Me velocity operator 
y wavefunction m ?? 
parity or space inversion 
P symmetry 
T time reversal symmetry 
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(continued). 
Symbol | Formula Quantity Unit 
1 
5- Fermi-Dirac distribution function 
f exp (4) +1 
g f+ôg distribution function 
ôg out-of-equilibrium part of g 
Quantum numbers 
n principal quantum number 
E spin quantum number 
m; spin projection quantum number 
l orbital quantum number 
orbital projection quantum 
mı number 
i total angular momentum 
J quantum number 
total angular momentum 
my projection quantum number 


Energies, forces and related quantities 


E€ energy J 

En energy of the nth band J 

Ep Fermi energy J 

A Zeeman energy J 
Eso spin-orbit coupling energy J 

J sd exchange constant J 
Me —edo electrostatic potential J 
ôn ' ; 

Un HUn + Ne) chemical potential J 
lm £,at T —0 equilibrium chemical potential J 

H Ue + Ln electrochemical potential J 

Hy —(u! — u^) spin s accumulation? J 

D electric potential V 
E —Vo electric field Vm! 
F force N 
aR Rashba parameter Jm 
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(continued). 
Symbol | Formula Quantity Unit 
Times and related quantities 
t time S 
Ts electronic relaxation time S 
sd Rl J sa exchange relaxation time S 
Taf spin diffusion relaxation time S 
Tsd 
é Tef spin mistracking parameter 
f frequency s! 
o 2nf angular frequency rad s! 
Lengths, areas, volumes, angles and related quantities 
d thickness m 
A area m? 
V volume m? 
L,Y, 2 cartesian coordinates m,m,m 
. , m,rad, 
r,0, Q spherical coordinates rad 
k wavevector m! 
kp Fermi wavevector m! 
extremal spanning wavevector of z 
kr the Fermi surface m 
As Upte electronic mean free path m 
lsd V Dt exchange length m 
Asf UPTsf spin-flip length m 
Ley Diss spin-diffusion length m 
-1/2 
sf —— du» average spin-diffusion length m 
lr De 
WDW domain wall width m 
Q, Vy x Qi) [iV [v (K)) Berry curvature of the nth band m? 
A, (v, (Kk) | Vx |v.,, (K)) Berry connection of the nth band | m 
Lo 21, 
Yn $ An(k) dk = [f Qn(k) dk Berry phase of the nth band rad 
C L p Q,,(k) ak Chern number of the nth band, 
4 27 integer 
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(continued). 

Symbol | Formula Quantity Unit 
Particle density and density of states 

n no + ôn electron density m” 

m l 

no n f(s) N(s) de equilibrium electron density m? 

ôn excess electron density m? 

N density of states Jim? 

p —e(y|v) charge density Cm? 


v electron velocity ms 

Up Fermi velocity ms! 
A "DR A 2-1 

D 3 UF Ae electronic diffusion coefficient mî s 


Conductivities, resistivities and related quantities 


c N(e,)e? D electric conductivity Sm! 
ot o(1 = B’) effective electric conductivity Sm! 
1 ES 
p P electric resistivity Om 
1 
p* pm effective electric resistivity Qm 
T spin ‘resistance’ (resistance area $ 
ply product) Om 
spin ‘conductance’ (conductance E 
c/T; : (4) Sm? 
sf per unit area) 
—A spin-coupled interface ‘resistance’ 
à He : Q 2 
Ts zJ (resistance area product) m 
č 
p! 
a of bulk spin asymmetry parameter 
f interface spin asymmetry 
: parameter 
ô interface spin-flip parameter 
a—l 
p a+1 spin polarization 


g nm nm i 4 
Gi DU ios > " (ri ) spin mixing conductance’ S 


meN 
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(continued). 
Symbol| Formula Quantity Unit 
Gi spin mixing conductance per unit 
g'! AG) /2 area A per quantum conductance m? 
per spin channel Go/2(? 
Gt Re(G"') real part of G"! 
a! 3Im(G") imaginary part of Gl 
_, | effective spin mixing conductance 
1 - d - unit : -2 
is E f tanh-! EN per unit area per quantum m 
n Goon Uy conductance per spin channel 
g! 1— G try rt) spin mixing coefficient 
=I 
h (Czy, Ortz, Tyr) Hall vector (S m ) 
Charge and spin current densities and related quantities 
, Vu ? 1 2 
J PA partial charge current density Am 
Je "ES total charge current density Am? 
J -( -j') total spin 8 current density?) Am? 
h total spin s current density?) Jm? 
o ad = Re(IS viv) 
Ju hxE total charge Hall current density Am? 
Je 
OSHE J. spin Hall angle 
Magnetizations and flux densities 
M Lig (n! — n!) magnetization Am! 
M, saturation magnetization Amt 
m —|y|8 magnetic moment Am 
m m/V magnetic moment density ?? Amt 
magnetic flux density or magnetic 
HoH B-ficld t 
HoH ete effective magnetic flux density T 
radiofrequency magnetic flux 
Ho hing densit * 
sity 
(v x E) 
Hol so men spin-orbit magnetic flux density T 
TDL damping-like magnetic flux 
TDL xi = m” x (m x M) density amplitude T 
m Ey s ULAXxM field-like magnetic flux density T 
Ms amplitude 
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(continued). 
Symbol| Formula Quantity Unit 
Angular momenta, torques and related quantities 
1 
8 z p” spin angular momentum Js 
s (W|S|p) = 8/V spin angular momentum density Jsm? 
Jt M (macro)spin angular momentum EM 
Su ly| density of a ferromagnet Jsm 
4 orbital angular momentum Js 
torque from a spin density on a 
Ty VT magnetic moment J 
1 dM torque from a spin density on a E 
T X ———— as (10) Jm" 
|p| dt magnetization 
t transmission coefficient 
r reflection coefficient 
g Landé ¢ factor 
e 
y NS gyromagnetic ratio) C kg” 
dM „Mx dM i r 
[ri di ed Ms di damping parameter 
aX A = Xo + Ap intrinsic damping parameter 
Lp extrinsic damping parameter 
Heat transport 
T temperature K 
m : AS Sm! 
a J.=-«:-VT thermoelectric conductivity VK! 
Sj bi = OHS i Seebeck coefficients VK! 
Skyrmion transport 
1 aT VW VW 5 é 
Q In M. (oM x Q,M ) dxdy| topological index 
Msd E 
G — " = An Qz gyromagnetic coupling vector kgs! 
) 
Msdr Be DL MT 
Dij bi (aim :0;M ) didj dyadic coefficients kgs! 
i ) 
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In the literature, the terms moment m and spin 8 are sometimes used inter- 
changeably, and the 1 and | symbols sometimes refer to m or 8, which can lead to 
confusion as m and 8 are of opposite sign from 8 = —4n/|y|. In this book, the term 
moment refers exclusively to m and spin to 8; the 1 and | symbols refer exclusively 
to the orientation of m, and the f; and |} symbols refer exclusively to the orientation 
of 8. Note that, in the literature, the + and — symbols are sometimes used for the 
orientation of 8. 

(I the literature, the spin accumulation is defined either as u, = uil — pl, or 


u, = — (ui — pi), or u, = x (u! — u!) /2, resulting in a factor +1/2 difference in some 
of the equations depending on the articles and books consulted. Throughout this 
textbook, the term spin refers to 6 and the spin 8 accumulation is thus defined as 
p, = ab — u = - (uw = ut), because the f and | symbols refer to the orientation of 
the moment m, whereas the ¢} and |J symbols refer to the orientation of the spin 8, 
which points in the opposite direction to m. 

(In the literature, Uy is most often written l. The notation l may implicitly infer 
that the spin-diffusion length is spin-independent, which is not always the case. 
Hence, we prefer to use a superscript as a warning. 

Tt is customary to incorrectly use the term spin conductance for o/ Us although its 


unit is S.m^?. Note that this is the same unit as GT JA. In practice, the spin 
; Dn 
conductance is often used per quantum conductance A It can be compared 


directly to g! = n 


(9Tn the literature, the spin mixing conductance is defined either as G^, in units of S 
as here and throughout this textbook, or as G / A, in which case it is in units of 
S.m~?, resulting in a factor A difference in some of the equations depending on the 
articles and books consulted. 

The coefficients tj) and ri; account for the transmission and reflection of a 


moment m by a ferromagnet of magnetization M. Equivalently, the coefficients ty) 
UN 

hl 
magnet of macrospin density Sy = -gM. Therefore, thy) = tj) and Tay = rq), 


2 * i ' 
e = FE net Eevr(m)]eimaq-mÍeer(m)] an 
gt+=1- («t + nri) =1- (tty + ni). 


(Tt is customary to incorrectly use the term spin mixing conductance for g!! = 


although its unit is m~?. 


and ry) account for transmission and reflection of a spin 8 = m by a ferro- 


qt 
AGo/2? 


(tm the literature, the spin current density Js = -( — J!) is expressed either in 
A.m^?, units of charge current density, as here and throughout this textbook, or in 


J.m"?, units of the density of spin angular momentum ^, in which case 


2: 
J.—-— i Gg = 3*5. resulting in a factor i difference is some of the equations 
depending on the articles and books consulted. In this textbook, we use Q — £J. In 
the literature, the spin current is also defined either as J, —j! —j', or as 


Js = -( = j^). resulting in a factor —1 difference in some of the equations 
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depending on the articles and books consulted. Throughout this textbook, the term 
spin refers to 8 and the spin 8 current is thus defined as J, =j" — j^ = -( — J*). 
because the f and | symbols refer to the orientation of the moment m, whereas the f| 
and |J symbols refer to the orientation of the spin 8, which points in the opposite 
direction to m. 

(In the literature, m is used either to define a property, in A.m^!, units of mag- 
netization, as throughout this textbook, or for the dimensionless magnetization unit 
vector, in which case m — A resulting in some possible misunderstanding and a 
factor Ms difference is some of the equations depending on the articles and books 


consulted. In this textbook, we use an overcaret for unit vectors, i.e., M = x 
(Tm the literature, the torque T = — Ba is expressed either in J.m^?, as here and 
throughout this textbook, or in A.m !.s^!, in which case T = — M, resulting in a 


factor |y| difference in some of the equations depending on the articles and books 
consulted. 

(UT the literature, the gyromagnetic ratio y = — >= , which is negative for electron 

, gy B y Im. I 8 

spin is used either as is in the equations, e.g. 4 = yM x (uH). .., or in absolute 
value ||, as throughout this textbook, e.g. d = —|y|M x (uoHerr)..., or replaced by 
Yo = —yuo, to make it positive and include wy to fit the cgs unit system: 
dM — yoM x Hef. .., resulting in a factor —1 or cj difference is some of the 


dM — — 
equations depending on the articles and books consulted. 


